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ABSTRACT 

In this paper Riemann-Lioville fractional integral and derivative formulas 
for Fox-Wright ,,¢,(z) generalized hypergeometric functions are obtained. Some 
special cases of the established formulas are also discussed. 
2000 Mathematics Subject Classification : 26A33, 45D05. 
Keywords : Fox-Wright generalized hypergeometric Function, Riemann-Liuville 
fractional Integral and differential. 

1. Introduction. The Fox-Wright generalized hypergeometric ,9,(z) for 
zeC is defined in series form as [4]. 

p 
(a,,A,), . Tr(a. + A;4)z 


9, (2) =, €. HESSE s 
p Yq p Ya (b B;), 120 [Irho 2,5): ALADO) 
11 


where a,,b, eC,A, >0,B; > 0 
q p A ; 
1+ DB, — XA > 0; A,,B, E R(A,B, TA O)(i 21. PJ 5 15-0) for suitably 
J=1 


jal 

bounded value of |z|. 

As special case, for Ay=...=A,=1, B,=...=B,=1, (1.1) reduces to generalized 
hypergeometric function [3] 


$ (2:3), (a1). © T(a,),....P(2,) F [encia 
NER). AE) r(b),..r(b) ^ 'tbu-»5: 
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(1.2) 


q 


I — 


For p=1, q=1, a, 26,4, 2 1,0, =P and B, =o, (1.1) reduces to 


(5,1) e T(5+k) 2" 6 
=) ———— — -I(6)E 
| D AA Lemon) ay | Pan) (0.3) 
called generalized Mittag-Leffler function defined by Prabhakar [2]. 
If §=1 then (1.3) reduces to generalized Mittag-Leffler function 
(1,1) = F(1+k) z^ 
= — —— = E 
iro z Apra) k! me) 
a If §=1 and B-1 then (1.3) reduces to 
IN) 2! ee 
oa era P0) ...(1.8) 


called Mittag-Leffler function [1]. 7 
The object of this paper is to derive the relations which exist between the 
Fox-Wright generalized hypergeometric function no (z) and left and right side 


Riemann-Liouville fractional integral and differential operators. 


The fractional integral and differential operators defined by Samko, Kilbas - 
and Marichev [7] for q. >0 are given by | 


| a 1 x a= 1 

| Were) eT DES "e E 

C 
(CA e Far m E 

r(o) x , "our ( 

ha aa J (c£) Agar (1.6) 

| " P(t-{o})\de) ^ , Be 

| [el í 

B - = 1 d al = ta} 

à Ba aa) (ey r()at SORD 

| where [a] means the maximal integer not exceeding « and {a} is the fractional | 


di part of a. 
| 


2-Properties of generalized Fox Wright hypergeometric function. I 


In this section we derive several interesting properti 


1.3) 


1.3) 


the 
ide 


bas 


1.4) 


1.5) 


1.6) 


Wright hypergeometric function 0, (2) defined by (1.1) with the help of Riemann 


-Liouville fractional integral and derivative formula [4]. 


Theorem 1. Let o>0,8>0,y>0 and acR,let IP be left sided operator of 


Riemann-Liouville fractional integral (1.4). Then there holds the formula 


(2,4), 1 E (A) yen) | 
5 at" (ea) = x? pom a 


"et ax 3m 
pra (b (8, ), (eB) eA 2711), 


J' PAS 


Proof. By virtue of (1.1) and (1.4) we have 


à p 1 “e Ir (a, as Ak) (x -ty d a, thr 

at (J^ rey 2 k! 

th Tra 
j=l 


Interchanging the order of integration and summation; and evaluating the inner 
integral with the help of beta functions, by setting (=xy, we get 


Hr (a, + Ak) ep : 
LAS. ya ———_ 7 “(1-97 y" dy 
^^ ro, + Bk) I (o)! 
J=1 


(a, A), , (1,8) 


aorta 1 : B 


nn (8,8), (cB) 


Interchanging the order of integration and summation is permissible under the 
conditions stated with the theorem, due to convergence of the integral involved in 
the process. This completes the proof of Theorem 1. 


Corollary 1.1 .For a > 0,8 » 0,y,A» 0, there holds the formula 


T a 
Te | 10; bas 


Corollary 1.2. By setting § — 1 in (2.2) there holds the formula 


ar o =1(5)2 Ej... (ax") (2.2) 


Ly, (E78, (at?) (x) = ° Epa (az) 3) 
establised by Saxena and Saigo [6]. 
Theorem 2. Let o.» 0,B>0,y>0 and a e R,let I" be right operator of Riemann 


Liouville fractionalcfmeegvattul Kayi Pestiothseekesl ds re aforinitalse 


-ammm 


4 
\ 


| (a A). s l ) E ta, A), p(B) | E 
| JU ep) at ge) Orr ax” 
i p®, (b,,B,), | pl t gel (b,,B,), (at yB) (2.4) j 
Y 5 d 4 4 J 
Proof. By virtue of (1.1) and (1.5), we have 
4 p , 
Pla, + Ak uc utes 
rte. at? | (x) ; lis H ai (tx) E - dt 
EM pYy ~ 1 R 4 277 3 
(b,.B, m i r(x) h=0 r(o, ES Bk) k! ; 


Interchanging the order of integration and summation and then evaluating the 
inner integral by beta function formula, we get 


f (Ken ie de 


a . J^ yal (x at, as dt (by putting t-x=y) 


(putting y=wx). 


Sy ieee 
=x pua ( 
j 


Interchanging the order of integration and summation is permissible under the 
conditions stated with the theorem, due: to convergence of the integrals involved 
in the process. This completes the proof of Theorem 2. 


| ` Corollary 2.1. For œ> 0,8» 0,y » 0,ae R, there holds the formula 


1 I = 191 ie 


Corollary 2.2. By setting § =1 in (2.5), there holds the formula 


|o SS) En (a a — (2:5) 


j tio Heid war. otri Initiative 
(uote X Baty (arty 


estabilished by Saxena and Saigo [6]. 


Theorem 3. Let o50,8> 0,70 >0s0aem and D; be left sided operator of 


(2.4 : : : 2 2 
! | Riemann Liouville fractional integral (1.6). Then there holds the formula 
(a A.) | | (a A | A 
Y 3 exister uo ; - pA) UG B) 
e OL » 9, ai at" (y= PEG ae lax" 
a d (6,,B, A HIR (b,.B, ), (eB) TA 
Proof. By virtue of (1.1) and (1.6), we have 
,44,), fajas (aA) | 
4- | d : 1? Lp 
DILO O, “at («)=(4] Tes p lat! (x 
the ie (bB) Sid ai BE ie 
| 
$ Pla, + Ak > 2 
Rh : 1 * zl £l "pegha dt. 
^ T [r (o, + B,2) F(1- fe) eA dx 
j=1 
Interchanging the order of integration and summation and evaluating the 
inner integral by beta function formula (by setting t=xy), we obtain 
Ir a, + Ajk) gi qui 
à 3 1 3k ‘ "ETT 
L-H.S.= ye : a r(y*p ) (=) NS 
E [[r(o, +Bjk) k! T(y+Pr+1-[0)) 
j=l 
p 
r the d JIT (a, n NES eat q 


lved = i=l 
lv Daan z 


(a; A; js (Y, B) 


SEL o 
(2.5) me" (b, B) (01,8) 
Interchaning the order of integration and summation is permissible under the 
condition stated with the theorem, due to convergence of the integral involved in 
th 3. 
| € process. This completes the proof of Theorem 
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EE Aa “—- 
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Corollary 3.1. For œ > 0,82 0, y, » 0, there holds the formula 


a (OY) 
|] Spin 


Corollary 3.2. By setting 5-1 in (2.8) then there holds the formula 


11 ) youl 3 
Dy. le 19) jea at! Je m EN (ax! ) (2.9) 


established by Saxena and Saigo [5]. 


Y D; 


04 


at lis x)=. (S)x yes A ¿(ax ) ...(2.8) | 


Theorem 4. Let o>0,8>0,y,A>0, ae R and D" be right sided operator of 


Riemann Liouville fractional integral (1.7). Then there holds the formula 


(a;, A), p (Y, B) -B 
D+ P, + ax 
priua (b. B,), (c. Y B) 


(a;,A;) 
Mat Mo) at 


(o, B;), 


Proof. By virtue of (1.1) and (1.7), we have 


Dry py * (2.10) 


aka (a;,A,) T 
"DT m at" a E O EA (x) 


dx (b, B; ) q 


SS e, 


= $E IC NAE (E) Te = x) p Meg 


| p 

li Unta, tAk) a! | di le re M ud 

n = re oe => x+y rs y Ik dy 
Ih aa) dx 0 l y 

D, j=l 


(putting t-x«=y) 


- Tr (a+ Ak) Ng 
y-[0]-1) 
O ARE, TON Ab 2 


d [edet uec 
-ph-H od 
ee Mitiative (putting y=WwX ) 


(2.9) 


or of 


2.10) 


yx ) 


Pla, + Ak 
S E (nse) a'T (Bk 4- y) 


per) II rb, E B k) kIT (Bk y — a) 


j=1 


- Bh 


7 (a,,A, (58) 
=x 19, ; 
i (6,,B,), (.—v.B) 


Interchanging the order of integration and summation is permissible under the 
condition stated with the Theorem, due to convergence of the integral involved in 
the process. This completes the proof of Theorem 4. 


Corollary 4.1. For 9» 0,8» 0,y,1 5» 0 , there holds the formula 


8,1 
llt 


Corollary 4.2. For 8-1, (2.11) takes the form 


at “feo = TS) "Ej... (ax) Ee 9 


E Li = E 
oe elt j i Jlez ‘Epu (a?) (2.12) 
established by Saxena and Saigo [6]. 3 
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ABSTRACT 

In the Present paper, an attempt has been made to bring basic 
hypergeometric functions within the purview of Lie theory by constructing a 
dynamical symmetry algebra of basic hypergeometric function 201- Multiplier 
representation theory is then used to obtain generating function for basic analogue 
of Gauss polynomial. The results obtained in this paper are extensions of the 
results derived earlier by Miller [3] and Sarkar-Chatterjea [4]. 
2000 Mathematics Subject Classification : Primary 60E07, 20G05; Secondary 
33C55, 14A17. 
Keywords : Lie algebra, Generating function, Basic Gauss polynomials. 

1. Introduction. The q-analogue of the Gauss functions or Heine's series 
[1] may be written as 

»0, (a,b; ¢;q; x) = Sila: q,n|[b;q,n]/ [c;q,n][n; al! ( (es 0,—1,—2,...) 

u=0 

where |g| «1 and |x|<1. 
Here [a;q,n] and [n;q]! are respectively the basic Pochammer's symbol and basic 
factorial function defined as [a;q,n]-[a;q]la*l;q]...[a-*n-1;q] and 


[n;q]!= [1;¢][2;q]...[zsq]. 
The basic differential Oper oa sx is defined by [1] through the relation 
2011) 


B, 0x) = 19 (ae) - 6(x))/x(a -1). 
2. The A Rf.edng o hdigamical symmetry 


R 
m 
i 


algebra of the hypergeometric function has been defined by Miller [2]. We use the 


same technique to define the dynamical symmetry algebra of ,6,. Let 


Qua = Pa (Y pa o)r, (a)/T, (y)-2 [o,8;v;q;x]s"u^t* s: C210) f 


be the basis elements of a subspace of analytical functions of four variables x,s,u 
and t, associated with Heine's basic hypergeometric function of Heine's series ,p,. 
Introduction of variables s,u and ¢ renders differential operators independent of 
parameters o,p and Y and thus facilitates their repeated operation. 

The dynamical symmetry algebra of ,0, is a 15-dimensional complex Lie 
algebra isomorphic to sl (4), generated by twelve ps -operators termed as raising 
or lowering operators in view of their effect of raising or lowering the correspond- 


ing suffix in dus - The ps -operators are 


(i) Ef =s (x(1- x)B; , +tB; —- sB” - xuB? hh 


AN qe qus qu 
(NEN. ut (x(1 2x)B,.—xsB, tB = 1), 29) 
The action of these operators on Yus is given by 
dle Bang = [e “4 Lalor, Bang 


E c aa Pahang = [Y -a- ASCH eem E09.) 


The upper factor in each bracket is to be associated with plus sign and 
lower with minus sign. Twelve E-operators together with three maintenance 


operators J,,J,,J., and Identity operator I form a basic for g1(4) =sl(4)U), where 


(D is the 1-dimensional Lie algebra generated by 1. 
Here 


Jr. = SB, Jo. ims UB? AS), = tB, and I^ “1 (2.4) 
with the results 


JA Ona = CA D 
Jo Dag = [Bid] Ouro ? 


Paba = [v:2]0,5,,., D 
and 


ey x: Danya a (2.5) o 


3. The Generating Functions for Basic Analogues of Gauss 
Polynomials. On comparing the results obtained by the action one parameter 
subgroup (exp, aE^, +4) generated by the operator E”, defined in (2.2) on Qapi 


-4,q 


defined in (2.15 AMPARA PERA 


tt 3 . “TI Y- z 
he (st / (ax 4 st) [al (st s JJ 1 O [0,B;¥;q;x(a + st )/(ax st)| 


y Li mia S > s * pe 
(2.1) ! = De [ ues ql, /[m;q]!,®, [oq "rale eng pio cts líst) (3.1) 
X,5,U m0 3 
0. Taking, œ > 0,B ~A+p+m-1,y 2 q"",st 3 1,a 31, we get 
nt of qe. 


(1+ x) = [2:q] = Val, /[m;q]/., b, [-m, ù +u + m-i; q; x (3.2) 


m0 


z Lie By definition of basic Gauss polynomial [1] 


^M ^ || 
ising | G"(q;x)- ,,|-m,.X * u* m - 1;3;q;x], | 
ond- | where 4 #0,-1,-2,3,... - 13.3) | 
Using (3.3) in (3.2), we get the generating function | 
1 | 
` l-i- re A-1 3H 
(1x) ""[2;g] = Y [v:a], [m4] C)" (a; x) | 
m=0 j 
| 
(22) for basic Gauss polynomials. | 
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ABSTRACT 

The aim of the present paper is to establish two new recurrence relations 
for the H -function. Some results for the generalized Wright hypergeometric 
function and gen Riemann Zeta function are special cases of our main 
findings. 
2000 Mathematics Subject Classification : 33C60 
Keywords and Phrases : Generalized Wright Hypergeometric Function, 
Generalized Riemann Zeta Function and H -function. 

1. Introduction. The j -function was introduced by Inayat Hussain [4] 
and studied by Bushman and Srivastava [1]. 


The -function is defined and represented in the following manner: 
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| 1 The behaviour of jj -function for small valkues of |z| follows easily from a result 
b recently given by Rathie {[6],p.306, eq.(6.9)7. 
| 2. Main Results. In this section, we establish two recurrence relations, ! 


First Recurrence Relation 
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Second Recurrence Relation 
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Recurrence relations (2.1) and (2.2) are valid under the conditions of (1.1) and 
(1.2). > 
Proof. The recurrence relations (2.1) and (2.2) can be easily established by app 


ing to the definition [4], [1] and use the well known relation (Rainville [5], poh 
Th.8) 


T 2 
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involving l'efupetion and reintespuetidhe mesubigini terms of the jj -function E 
arrive at the right hand side of the desire results (2.1) and ' 


3. Particulars Cases. 


res : : Neha PA AT oe 
ult (1) If in the recurrence relation (2.1) we reduce Hj) to generalized Wright 


hypergeometric function (3, p.271, eq. (7)] 
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(ii) Ifin the recurrence relation (2.2), we reduce Hy," to generalized Riemann 
..(2.1) Zeta function (3,p.27, eq.(1); 7, p.314-315, eq. (1.6) and (1.7)] 


e 


21.0 (0,1),(—n,1;P) 


lz, p,n) = 1 jem ae | 
(0,1,1),(1- n,1; P) 


i o7 axi ET V2 —gg "xu 
; à * — (o1, (-m;P) 


ii) If A;=B;=1 in (2.1) and (2.2) then [7 reduces to Fox H-function and we get 
known recurrence relations given earlier by Chaurasia [2]. 
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ABSTRACT 

A large body of literature is devoted to the estimation of covariance matrices 
in a large-scale setting. Recent work in this area includes the shrinkage approach 
proposed by Schäfer and Strimmer [16], where the authors analystically calculate 
the optimal shrinkage intensity, yielding a good, computationally inexpensive 
estimate. Our focus is an estimate with the property that the corresponding inverse 
covariance matrix is sparse. 

Dempster [4] introduced the concept of covariance selection, where the 
number of parameters to be estimated is reduced by setting to zero some elements 
of the inverse covariance matix. Covariance selection can lead to a more robust 
estimate of S if enough entries of its inverse are set to zero. Traditionally, a greedy 
forward/backward search algorithm is employed to determine the zero patterns 
Lauritzen [9]. However, this method quickly becomes computationally infeasible 
as p grows. 

2000 Mathematics Subject Classification : Primary 90C25; Secondary 90C47, 
49K30 

Keywords and Phrases : Applications of Mathematical Optimization, 
Robustness, Duality and Bounds, Convergence and Property of Solut:.:: Convex 
Optimization, Nesterov's Method. 

1. Introduction. In this paper we investigate the following related idea. 
Beginning with a dense empirical covariance matrix S, we compute a numum 
likelihood estimate of S with an II-norm penalty added to encourage sparsity in 
the inverse. The authors Li and Gui [10] introduce a gradient descent algorithm 
in which they account for the sparsity of the inverse covariance matrix by defining 
a loss function that is the negative of the log likelihood function. Recently, Huang, 
Liu and Pourahmodi [7], and Dahl [2] considered penalized maximum liklihood 
estimation, and Dahl [2] in particular, proposed a set of large scale matiada to 
solve problems where a sparse structure of S-1 is known a priori. Qua contribution 
is threefold: we presentsaquaeablycenamer gendualgaridangothakidmefficient for large- 


a sparse, invertible estimate of S-1, even for n<p; we 
exity estimates for the problem; and finally we test our 


scale instances, yielding 
obtain some basic compl 


algorithm on synthetic data as well as gene expression data from two datasets, 


Notations. For a pxp matrix X,X20 means X is symmetric and positive semi- 
definite; |x| denotes the largest singular value norm, |x]; the sum of the absolute 


values of its elements, and |x|], their largest magnitude. 

2. Preliminaries. In this section we set up the problem and -discuss some 
of its properties. 
2.1. Problem Setup. Let S20 be a given empirical matrix, for data drawn froma 
multivariate Gaussian distribution. Let the variable X be our estimate of the 
inverse covariance matrix. We consider the penalized maximum- -liklihood problem 


max log det X -(X,X)- pls, GU) 
where (S,X = trace( SX ) denotes the scalar product 

between two symmetric matrices S and X, and the 

term lle, = Y, A penalizes nonzero elements of X. 

Here, the scalar parameter p>0 controls the size of the penalty, hence the sparsity 
“of the solution. The penalty term involving the sum of absolute values of the entries 
of X is a proxy for the number of its nonzero elements, and is often used-albeit 
with vector, not matrix, variables-in regression techniques, such as LASSO in 
Tibshirani [17], when sparsity of the solution is a concern. 

The classical maximum likelihood estimate of X is recovered for p>0, and 
is simply S, the empirical covariance matrix. Furthermore, as noted above, for 
p>>n, the matrix S is likely to be singular. It is desirable for our estimate of 5 t0 
be invertible. We shall show that our proposed estimator performs some 
regularization, so that our estimate is invertible for every p» 0. 

2.2. Robustness, Duality and Bounds. By introducing a dual variable U, W° 
can write (1) as max min log det X + (X,S+U). 


Here |U||. denotes the maximal absolute value of the entries of U. This correspond? 


to seeking an estimate with maximal worst-case likelihood, over all component 
wise bounded additive perturbations S+U of the empirical covariance matrix s. 
Such a "robust optimization" interpretation can be given to a number of estimation 
problems, most notably support vector machines for classification. 

We can obtain the dual problem by exchanging the max and the min: 


min E log det {S$ UA LaS. tl aud? bv eGangotri Initiative 0 
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corresponding covariance matrix estimate is -$ :=S+U. Since the above dual 


problem has a compact feasible set, the primal and dual problems are D ub 
The optimality conditions relate the primal and dual solutions by Y MEN 


The following theorem shows that adding the /,-norm penalty regularizes 
the solution. 


Theorem-1 For every p>0, the optimal solution to the penalized ML problem (1) 
is unique, and bounded as follows : 


1 
a(p) € X <B(p)l, where 0(p)===—B(p)= i 
IS] pp p 
Proof. An optimal X satisfies X-(S--U)-!, where |U| <p. Thus, we can without 


loss of generality impose that X 2 a(p)I, where o(p) is defined in the theorem. 
Likewise, we can show that X is bounded above. Indeed, at optimum, the primal- 
dual gap is zero: 
0- -logdet(S 4 P)- p -logdet X * (S, X) * p|X]| 

--p*(S,X)*p|X|,, 
where we have used (S-- U)X— 1. Since S,X are both positive semi-definite, we obtain 


|X| SIX Ie II, <B(p)7 as claimed. Problem (2) is smooth and convex. When 


p(p+1)/2 is in the low hundreds, the problem can be solved by existing software 


that uses an interior point method Vandenberghe [18], The complexity to compute 
an e-suboptimal solution using such-second-order methods, however, is 


Olps log(1/e)) making them infeasible for even moderately large p. 


The authors Dahl et al. [2] developed a set of algorithms to estimate the nonzero 
entries of x-! when the sparsity pattern is known a priori and corresponds to an 


undirected graphical model that is not chordal. Here our focus is on relatively 
large, dense problems, for which the sparsity pattern is not known a priori. Note 
that we cannot expect to do better than O(p?), which is the cost of solving the non- 
penalized problem p=0 for a dense sample covariance matrix S. 

2.3 Choice of Regularization Parameter p. In this section we provide a simple 
heuristic for choosing the penality parameter p, based on hypothesis testing, We ' 
emphasize that while the choice of p is an important issue that deserves a thorough 
investigation, It is not the focus of this paper. 


the observation that if p< Is, E then there cannot be zero 
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The heuristic is based on. 


an. 


in that element of our estimate of the covariance matrix Xj #0 suppose we choose 


p according to 
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where ¢,_,(y) denotes the two-tailed 100y% point of the ¿-distribution, for ng 4 


n-2 


degrees of freedom. With this choice, and using the fact that S20, it can be shown! € 
that p< |S, | implies the condition for rejecting the null hypothesis that variables 


i and j are independent in the underlying distribution, under a likelihood ratio 
test of size y Muirhead [13]. We note that this choice yields an asymtotically | 
consistent estimator. As n > eo, we recover the sample covariance S as our estimate 
of the covariance matrix, and S converges to the true covariance 2. 

3. Block Coordinate Descent Methd. In this section we present an 
efficient algorithm for solving the dual problem (2) based on block coordinate ( 
descent. : | 
3.1 Algorithm. We first describe a method for solving (2) by optimizing over one < 
column and row of S+U at a time. Let W:=S+U be our estimate of the true ! 
covariance. The algorithm begins by initializing WW=S+pl. The diagonal elements | 
of W? are set to their optimal values, and are left unchanged in what follows. | 

We can permute rows and columns of W, so that we are optimizing over the 


last column and row. Partition W and S'as 
W- E Wiz | S= Ea Sy» | 
Win o> Sey Shy 


where wS e R’ the update rule is found by solving the dual problem (2), with 
U fixed except for its last column and row. This leads to a box-constrained quadrati? 
program (QP): 

Wp arg min y Way :Iy- ss]. s p} d 
We cycle through the columns in order, solving a QP at each step. After each sweep 
through all columns, we check to see if the primal-dual gap is less than €, 4 give 
tolerance. The primal variable is related to W by X-W-. The duality gap conditi | 
is then 


(S, X) p|X], < pre. | 


3.2. Convergence and Property of Solution. Iterates produced by oh 


: C0. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative Y 
coordinate déscent dorem are strictly positive definite. Indeed, since S20, 


(3) 


for n- 
shown 
riables 


d ratio 
tically 
timate 


ent an 
~dinate 


ver one 
e true 
ements 
WS. 

ver the 


that the following Schur complement is positive : w,, —w,W,'w,, >0- By the update 


rule (4), we have 
A er pure 
Way — Ú, War Wio > Wy — Wio Wii Uy, >0 


, 


which, using Schur complements again, implies that the new iterate satisfies W »0- 


Note that since the method generates a sequence of feasible primal and dual points, 
the stopping criterion is nonheuristic. As a consequence, the QP (4) to be solved at 
each iteration has a unique solution. This implies that the method converges to 
the true solution of (2), by virtue of general results on block-coordinate descent ' 


algorithms Bertsekas [1]. 


The above results shed some interesting light on the solution to problem 


(2). Suppose that the column sj, of the sample covariance satisfies |s,, € p, where 


the inequalities hold component wise. Then the corresponding column of the 


solution is zero X,, = 0. Indeed, if the zero vector is in the constraint set of the QP 


(4), then it must be the solution to that QP. As the constraint set will not change 
no matter how many times we return to that column, the corresponding column of 
all iterates will be zero. Since the iterates converge to the solution, the solution 
must have zero for that column. This property can be used to reduce the size of the 
problem in advance, by setting to zero columns of W that correspond to columns 
in the sample covariance S that meet the above condition. 

Using the work of Luo and Tseng [11], it is possible to show that the local 
convergence rate of this method is at least linear. In practice we have found that a 
small number of sweeps through all.columns, independent of problem size p, is 


sufficient to achieve convergence. For a fixed number of val ST e 


method is O(Kp*), since each iteration costs OD), 7 


3.3 Connection to LASSO. The dual of (4) is 
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denote the square root of Wib := aay" , then we can write(5) as min|Qx -bj * olx], 
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The above is a penalized least-squares problem, often referred to as LASSO. 


If W,, were a principal minor of the sample co 


equivalent to a penalized regression 


variance S, then the above would be 


of one variable against all others. Thus, the 


approach is reminiscent of the approach explored by Meinshausen and Buhlmann 


[12]. but there are two major differences. First, we begin with some regularization, 
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in particular, W}; is never a minor of S. In a sense, the coordinate descent methog 
can be interpreted as a recursive LASSO method. 

\ 4. Nesterov's Method. In this section, we apply the recent results due to 

Nesterov [15] to obtain a first-order method for solving (1). Our main goal is not | 

to obtain another algorithm, as we have found that the coordinate descent is already 

quite efficient; rather, we seek to use Nesterov's formalism to derive a rigorous | 

complexity estimate for the problem, improved over that delivered by interior point 
methods. | 
As we shall see, Nesterov's framework allows us to obtain an algorithm 
that has a complexity of O(p" le), where s>0 is the desired accuracy on the ' 
objective of problem (1). This is to be contrasted with the complexity of interior- 
point methods, O(p* log(1/e)). Thus, Nesterov's method provides a much better | 

dependence on problem size, at the expense of a degraded dependence on accuracy. — , 

In our opinion, obtaining an estimate that is accurate numerically up to dozens of 

digits has little practical value, as it is much more important to be able to solve 

larger problems with less accuracy. Note also that the memory requirements for 
Nesterov's methods are much better than those of interior-point methods. 
4.1 Idea of Nesterov's Method. Nesterov's method [15] applies to a class of 


non-smooth, convex optimization problems, of the form ] 
mind (x) :xeQ,} (6) 

where the objective function is described as | 
f(x)= f (x)+ max{(Ax,u), :ue Q,}. 


Here @, and Q, are bounded, closed, convex sets, f(x) is differentiable | 
(with Lipschitz continuous gradient) and convex on Q;, and A is a linear operator , 


Observe that we can write (1) in this form if we impose bounds on the eigenvalues 
of the solution, X. To this end, we let 


4 Q ={X sol « X «pr), 

| Q, - (U :|U|. <p), 

where o,B(0 « o « B) are given. We also define 
f(x) - -logdet X - (S, X), and A=pl. 

To Q, and Q, we associate norms and continuous, strongly convex function? « 


called prox-functions, d, (X) and d,(U). For Q, we choose the Frobenius norm: 


an a prox-function d, (X js —log det X + logB: For Q,, we choose the Frobenit® : 
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Soros. The above function turns out to be a smooth uniform 

"pont differentiable, convex on Q, and has a Lipschitz-continuous gradient, with a 
constant L that can be computed as detailed below. A specific gradient scheme is 

rithm ^ then applied to this smooth approximation, with convergence Wi Ot Lc). 

on the 4-2 Algorithm and Complexity Estimate. To detail the algorithm and compute 

terio. the complexity, we first calculate some parameters corresponding to our definitions 
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approximation. It is 


better above. First, the strong convexity parameter for di(X), on Q, is c, =1/(%, in the 
curacy. | sense that V'd, (X)[H,H]= trace(X HX" H)z p^|H|;. for every symmetric £7. 


Zens of 


O Furthermore, the center of the set Q; is X, =angminy.,, d,(X)= fM , and satisfies 


nis for d, (X,) » 0. Without choice, we have D, = max... d,(X) - plogp/a. 

lass of Similary, the strong convexity parameter for d5(U) on Q; is c, =1 and we 
have D, = max,;,o,d,(U)= p^ /2. With this choice, the center of the set Q; is 

(0. U, = arg ming., d;(U) - 0. 

For a desired accuracy e, we set the smoothness parameter u =¢/2D,, and start 
with the initial point X, - I. The algorithm proceeds as follows. 
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A M =1/o° is the Lipschitz constant for the gradient of fi and the norm ¡Al is 


| induced by the Frobenius norm, and is equal to p. The algorithm is guaranteed to 
produce an ¢-suboptimal solution after a number of steps not exceeding | 


|D,D, 1, ¡MD 
NG) AAN oto IR 8 
Voo, € Y 0,2 


z kd p(log k) (4 Kari Ye) 
' € i 


where k=(/c is bound on the condition number of the solution. 

Now we are ready to estimate the complexity of the algorithm. For step 1, 
the gradient of the smooth approximation is readily computed in closed form, via 
the computation of the inverse of X. Stcp 2 essentially amounts to projecting on | 
i Q4, and requires an eigen value problem to be solved: likewise for step 3. In fact, 
each iteration costs Op). The number of iterations necessary to achieve an 
objective with absolute accuracy less than e is given in (8) by, N(e)= Op" le), if 
the condition number k is fixed a priori. Thus, the Complexity of the algorithm is 


O(p** le). 


5. Numerical Results. In this section we present some numerical results. 
We begin with a small synthetic example to test the ability of the method to recover: 
a sparse structure from a noisy matrix. Starting with a sparse matrix A, we obtain ' 


S by adding a uniform noise of magnitude g=0.1 to 4”. 


E In figure 1 we plot the sparsity patterns of A, S-!, and the solution X to (1) 
h using S and p=0. 


| p 


A 


A S X 
Figure 1. Recovering the Sparsity Pattern | 
We plot the underlying sparse matrix A, the inverse of the noisy version of A^ 
and the solution to problem (1) for p equal to the noise level. | 
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Figure 2. Recovering Structure 
Gene Expression Properties 


The 300 experiment compendium dataset contains n=253 samples with p=6136 
variables. With a view towards obtaining a very sparse graph, we set y=0.1 in the 
heuristic formula (3) of section (2.3) to obtain p=0.0313. 

Applying the property of the solution discussed in section (3.2), the size of the 
problem was reduced to P =537. Three sweeps through all columns were required 
to achieve a duality gap of = 0.146 , with a total computing time of 18-minutes, 
34-seconds. The resulting estimate of the inverse covariance matrix $-!is 99% 
Sparse and has a condition number of 21.84. Figure (4) shows a sample subgraph 
obtained from $-!, generated using the Graph Explore program developed by Dobra 
and West [6]. The method has picked out a cluster of genes associated with amino 
acidmetabolism, as described by Hughes et al. [7]. 


Figure 3. Application to Hughes Dataset Using p-0.0313. 
As we have seen, the penalized maximum likelihood problem formulated here is 
useful for recovering a sparse underlying precision matrix y-! from a dense sample 
covariance matrix S, even when the number of samples n is small relative to the 
number of variables p. In preliminary tests, the method appears to be a potentially 


valuable tool for Sr BARS eiplresstotlvaata sGdrdestiiietfarther testing is 


E 


required. 
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ABSTRACT 

Takano [13] studied decomposition of curvature tensor in a recurrent space. 
Sinha and Singh [12] defined and studied defined decomposition of recurrent 
curvature tensor fiela in a Finsler space. Negi and Rawat ([2][3],[4]) studied 
decomposition of recurrent curvature tensor field in a Kaehlerian space. Rawat 
[5], Rawat and Silswal ([6],[7]) defined and studied decomposition of recurrent 
curvature tensor fields in a Tachibana space. Further, Rawat and Dobhal ((8],[9]) 
studied decomposition of recurrent corvature tensor field in a Kaehlerian recurrent 
space. Rawat and Singh ([10],[11]) studied the decomposition of curvature tensor 
fields in a Kaehlerian recurent space of first order. 

In the present paper, we consider the decomposition of curvature tensor 
field RA! jp in terms of two non-zero vectors anda tensor field. Also several theorems 
are established and proved therein. 

2000 Mathematics Subject Classification : 53A45 
Keywords and Phrases : Decomposition, Kaehlerian, Recurrent, Projective 
curvature tensor. 

1. Introduction. When in a 2n-dimensional real space Xy, of class 
C'(r»2), there is a mixed tensor field F; Ri; =1,2,3,...,2n satisfying 

F'F} =A}, : (E1) 
we say that the space admits an almost complex structure and we call such a space 
an almost complex space. 

If an almost complex space has a posit 
ds? = Ejd5/d& which setisfies 

FE gy, = Bj 


t 


Then the space is called an almost-Hermitian space. 


ive definite Riemannian metric 


...(1.2) 


etric(or skew-symmetric) in i and h. 


1 In * def l 5 D 
this case the tengo &urakdf] él CRS , Haridwar. An eGangotri Initiative 


If an almost-Hermitian space satisfies 
Verte Wy Da V,E, =0, 


«+» (1,3) 


where V, denotes the operator of covariant differentiation with respect to the 


3 : eae M] 
metric tensor g, of the Riemannian space then it is called an almost-Kaehlerian 


space and if it satisties 


VF, +V, F, =0, 


roth 
then it is called a K-space. 
In an Almost-Hermitian space, if 


V, E, =0, or F, =0, 


ihj 


then it is called a Kaehlerian space. 
The Riemannian curvature tensor field is defined by 


Ra 2354-9, 0003 ai aah 
where 9, =0/dx' and Ix ‘| denotes the real local coordinates. 
The Ricci tensor and the scalar curvature are given by 
R, =R; and R=R,g" respectively. 
It is well known that these tensors satisfy the following identities 


Rira =R- R (NI 


R, -2R; 


ia 


F*R,, - -R,F" 


ia^ j 
and NA = TN. 
The Berres projective curvature tensor is defined by 


PS 


h 
uk 5 = Rir + 


CS may et - R18! +S,,F} SE" +28, F/), 


where S; =F Rj. 
The Bianchi identities in K,, are given by 


h h h 
Ri, + Ry; Ry 50, 
h h h 
and Rin, a + Rira J + Roja = i 0 . 


(1.4) 


(1.5) 


A) 
...(1.8) 
20159) 


(1.10) 


(1 


(1.12) 


(1.18) 


The commutative formulae for the curvature tensor field are given as follows 


Wir -Ty =W R 


(1,14) 
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(1.5) 


(1.6) 


A) 
18) 
O) 


(1.10) 


„(1.10 


(1.12) 


(1,19) 


OWS 


(34) 


h hk _ ma ph h 
and Tru rp tiie vs 1, Rani =I TR 


umi? (1.15) 


A Kaehlerian space K, is said to be Kaehlerian recurrent space of first order, if its 
curvature 


h 3 ph 
V, RA =R 


ati» 


Le. R; 


UTR 


4 h 
= An Ris, (1.16) 
where 4, is a non-zero vector and is known as recurrent vector field. The space is 
said to be Ricci-recurrent space of first order, if it satisfies the condition 
Ria = ùR; La.) 
Multiplying the above equation by g", we have 
R, -À,R. Us) 
2. Decomposition of Curvature Tensor Field mn . We consider the 


ope d .. h . . 
decomposition of recurrent curvature tensor field Kj, in the following form 


Rin =U" OM 4 (2.0) 
where two vectors u^, d,and tensor field W,, are such that 
dl 2:913) 


Theorem 2.1. Under the decomposition (2.1), the Bianchi identities for R, take 


the forms 
TA IN E (2:3) 
and A.V RA EAS =0. (2.4) 
Proof. From equations (1.12) and (2.1), we have 
OY jp O Va 01V; =0. Xs 
Since v" z0- 
From equations (1.13), (1.16) and (2.1), we get 
v" à; [A ws +h Wig t AM =0. (2.6) 
| Multiplying (2.6) by X, and using (2.2), we obtain 
$; La + Wis + Way |= 9: (2.7) 
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x refore we derive 


www m O 


we 


AWW, FAW, FA Ya =O. 


| This completes the proof of the theorem. 
D] 


Theorem 2.2. Under the decomposition (2.1), The tensor field JUSTE and v.i 
satisfy the relations 


AVR, = XR, i AR, =O, jh ...(2.8) 


a yk 


Proof. With the help of eguations (1.7), (1.16) and (1.17), we have 
A,B, = AR, T XR (2.9) 


a^ uk 


Multiplying (2.1) by A, and using relation (2.2), we get 


Rin = OW jy ...(2.10) 
From equations (2.9) and (2.10, we derive the required relation (2.8). 
T Theorem 2.3. Under the decomposition (2.1), the quantities A, and y" behave 
the recurrent vectors. The recurrent form of these quantities are given by 


SMA AON 


am ma 


and v^ 5 Lu (2,12) 


EI 


Proof. Differentiating (2.8) covariantly w.r.t. y" and using (2.1) and (2.8), we obtain 


Mor OW ; = Ai tt = Àj Li D G3 
| Multiplying (2.13) by A, and using (2.1) and (2.9), we have 
^ kam (AR p -AR,)- he (X, Ro OR QM 


Now, multiplying equation (2.14) by X, we get | 
Na (XR; e AR) An = VA, (rin Ry us X Fa) . (2.15) 
Since the expression on the right hand side of above equation is symmetric in 4 


and h, therefore 

Mam Ba = Midas (2.16) 
provided that 

AR. -À R, #0. 


i jh Jo sik 


The vector field à, being non-zero, we can have a proportional vector Hm such 
that 

A, m = IA, (21 
Further, differentiating the egnation(9,9) Mirar Mehave 


; Aut rA, v" = 


..(2.10) 


behave 


..(2.11) 
..(2.12) 
obtain 


..(2.18) 


(2.14) 


..(2.19) 
"c in ( 


(2.16 


such 


m 


(217 


Making use of equation (2.11), we derive 
Ia E 


v, ==H, 0" (since X, #0). (2.18; 
This proves the theorem. 


Theorem 2.4 Under the decomposition (2.1), the vector ó,and the tensor Y 
satisfy the equation 

IAS Qian (2.19) 
where Wirm = Vj one 


Proof. Differentiating (2.1) covariantly w.r.t. x” and using equations (1.16), (2.1) 
and (2.12), we get the required result of the Theorem. 

Theorem 2.5 Under the decomposition (2.1), the curvature tensor and 
holomorphically projective curvature tensor are equal if 


Vim (05 -68/) + 0, (F/F! -F/F/))+20/4, ELE! =0. (2.20) 
Proof. The equation (1.11) may be written in the form 

Pi, = Rin + Dj, (2.21) 
where 

Dj, = ee = R8; 75 S $, S,EF/ E 2S, F"). ...(2.22) 
Contracting indices h and k in (2.1), we have 

R= v" OA ju. (2.23) 
In view of equation (2.23), we get 

(2.24) 


S; 5 Elo" ion > 
Making use of relations (2.23) and (2.24) in equation (2.22), we have 


A inh -08A (ER me) +20" Ov PER JO 25 


h 
ijk 


From (2.21), it is clear that BA = Ri if DM =0, 


which, in view of (2.25), becomes 
h un hy € 
v" Wan (0/51 - 0,88) +0, (RIF EM Ep) * 20" OW FF! = 0 p) 


Multiplying the above equation by A, and using relation (2.2), we obtain the 


required result (2.20). . 
Theorem 2.6 Unden.the.decompesition. A, khe scalar curvature R, satisfies the 


QM —ÀÀ—À 


es m O 
va 


i IR c iR = Ou 
| Proof. Contracting indices h and k in (2.1), we have 
R, =) ANa: «(212 
Multiplying (2.27) by g cU both sides, we have 
R= g'v QN 4: ...(2.28) 
Multiplying (2.28) by X, and using (2.2), we get 


A,R- BOM jy 
or, R,= GO y, [by using (1.18)] . 
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paN ABSTRACT ` 
The object of this paper is to derive integral involving modified H-function 


of several variables. These integrals are used to establish the Fourier series for 
Jim generlized function. By suitably specializing the coefficients and the parameters 
we can obtain many (new and known) interesting results involving multi-variable 
Cien. H-functions [8]. The results obtained by Srivastava and Panda [8], Kaul [1,2], Mac 
Robert [3] and Sneddon [7] follow as particular cases of our results. 
erian 2000 Mathematics Subject Classification : 33C70, 42A24 | 
Keywords and Phrases : H-function, Multi-variable H-function, Modified Multi- 
s ina variable H-function. 
1. Introduction. The modified multi-variable H-function employed as 
iban? ^ kernel of multi-dimentional transform defined by Prasad and Singh [6] on the 


lines of Srivastava and Panda [8], Prasad and Maurya [5] is as follows: 


| 
4 
| 
| 
i 
in a | 
| 
| 
| 
| 


rent | | 
Ji r XT ry fe 5 at ry irk i 

ensor m adf: Z (a;; Oa. Bo F :(e, HU Ej" sU u; ry Jr IR Gps 15), ee ee (G yy 13 n 
ens y n : ie Aa M, E 3 ee ba | E | 
ences, PyAr|R: pyqy ii p, My " (o ; AW sn BO L, A -U e Sun jf je : (2595). idi [ ls | 
Ids in age | | 
37069. = i 5 afo (E,)...0, (Sr WE sena 2-27 06:06, xls) i 
rent (2no) ;, i, à il 
where | 

ls in a: - i) 
1 Tr(a", ro cy 0E,) | 
Ye $ | 
$; (&)= a Jz | z (a m (i=1, De ey r) (1.2) it 
£ 5 1 
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is 


| o gorro de) iso) 
x Ir bi + Se) 11 ^ ji rfe - Soe) Tiro Sos) ! 


j=l 


LS) 
The multiple integral (1.1) converges absolutely if 
| 1 
' larg,[< Us, (il 2 m). 
where 
2 i i 1 i (i (7) A 
UNE y pRa- D zs Xr 
J=1 Janet) Jan m+ 
IR 
D ai dir >0 (i=1,2,...r). 
Jan +1 
We shall require the following known iritegrals. 
1. — ['cosP0(cos0/2)" (sing/2)” de 
-T(Pp*1/2)r(p, 1/2) y [py *1/2,-P,-P 1/2; 0 
2r(P+p+p, +1) 2) —P—p+1/2,1/2; 


provided that Re(2p+1)>0, Re(p, +1)>0 and P=0,1,2,.... 


2. f sin(2s+ 1)9(cos0)” (sino) de 


iN _T(2s+2)P(s+p+1/2)T (p, +1) pı +1,=s.(=s +1/2); 
h T(s+p+p,+3/2)P (2s5+1) “ *|(-s-p+1/2),3/2; 
la provided that Re(2p+1)>0, Re(p, +1)>0 and s= 0,1,2,.... 


2. Main Integrals. 


" 2p 29; 2h XS 2h. 2k, ( 
(D NC) (sing) 2l) (sing) C) (sns) =| 
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| : 
VE oca E (1/2=p5h,,....h,) (1/2=P,:k,,....k,), 
| | 2, (95h sh, e), 
(3) oM (r) . un o” tr) gr) Dela ") Lir)! j 
(a, Ae " E (e, j 458 jas; 8j I: 1 (1). Path ae Ux 
Mn (r) [T meg (r) pir) LEA RO sales 20241 
(o, Bi B7). (1:055 UP) A E le: 
provided that Re 2p +2), ha.) 20, Fel 2p, +2) kh, + i 20, 
t=1 i=] 
where A, = min Re(df [87)) (Uc Te RM oan 
Gi) O<@<n 
along with the conditions for the convergence of modified H-function. 
2p A20; 2h 2i 2h, 2h 1 
an fe cosPef eost) (ing) a (6) (sing) 22085) (sing) «| de 
: (-1)” (zB) (52 + 1/2), POR samy an, usum, um 
= 1 TET A : ALAN uy 
mau m (1/2), N! am x. 
z, LE (a Sor GP 
e RE AR AR Lp 
f (-P -p-p;;A, th, tk), (b; : Pr) 
LES" r) tr. ARTEL AMA 
(e, UL Ere 5 AG Wi Jee 
^g r) etr 1 S7 [g0 str) ...(2.2) 
(os UAR UD fi js d (4;,85),,, i (as Š; je 


with the conditions of (2.1) and P=0,1,2,... - 


jJ de (IIT) IN sin(2P +1)@(cos6) ° (sin 0)" H | [cose (sin®) 4 (cos0)"" (sin0)” z, y] de 


= (op +1) y dll EAM ar Hips 


<i 


€" 


ie 


De Qu . 
d (1/24 N -P prb (EN - piis ient is : 
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g j 2 E : 
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(in a), ALS, ED 


provided that Rel + Y ko; * | 20. 


izl 


where Q; = min Re(d\) 189),(j herent =152;....7) and P =0,1,2... 


3. Proof. In order to prove (2.1), we substitute for the modified multi- 
variable H-function in terms of its contour integral of Mellin-Barenes type and 
change the order of integration (which is permissible). We then evaluate the inner 
integral using Gamma function. On interpreting the resulting contour integral 
by means of (1.1), we obtain (2.1). Proceeding in the similar way, using (1.4) and 
(1.5) we obtain (2.2) and (2.3). 

4. Fourier Series. We shall now establish following Fourier series fot 
mdified multi-variable H-functions : 


(cos8/2)” (sin8/2)” HI (cose / 2)" (sin®/ 2)" Z,,...,(cos0/ gy (sin®/ 2)" zA | 
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1 o AO... QC; : 
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z, (-p— pi; hy + pst, +R) NCI Pia obi Au 
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y (b; B. B le s (6, Oo A js (28). mala 1 8^ de «(4.1 
with the conditions as in (2.1) and P=0,1,2.... 
5. Proof. In order to prove (4.1), let us suppose | 
E019) f (8) = (cos0/2)" (sing; 2)" (5.1) | 
| 
H| [(coser2y^ (sin@/2)"" z,,...(cos0/2)'" (sin0/2)”" Z jj | 
1 Z | 
- 5C + 2.6, cos PO. 4912) | 


Integrating (5.1) between the limits 0 to 1 and using the reult (2.3), we 


multi obtain the value of C,. Again multiplying both the sides of (5.1) by cosPe@ and 


e an 
: integrating from 0 to 7, w.r.t. 0, using (2.2) we get 


> inner 
1tegral 2 
4) and C, = =» (1/2), S! p*2,4 VIR: py yi P, Hd, | 
: | 
ies for z | 
. 1/24 s— P— ph, hi, (12- s piss): | 


| 
| 
| 
| 
| 
eo (-1), (-P), (-P 3r 1/2). Eee My See My My i 
| 


(-P=p—pysh + hc +h): | 


ip 
| . ye | 


(a, s one (e; 5 WB) d (c5), ico 
(CEU UMEN ON 
of C, and C, in (5.2). We obtain 


...(5.83 
Ot hs id 
with the conditions as in (2.1). Putting the values 
(4.1). Fourier sine series may also be obtained similarly. 

6. Special Cases. 
1/2), Case I. If we put m = IR'= IR=0 
7 | H-function defined by Srivastava and Panda [8] as : 


in (2.1) we get result in terms of multi-variable 


n P 9 $ 2h 2h, B 2h, i 
f(eos0/2)"(sino/2)!* H[(cos0/2)"™ (sin@/2)"*,-.(e0s0/2)"(sin@/2)2,]a@ | 


edes irl : 
& [1/2- pibus A UP ihe (2 P le . 


= pgosenm, megy Gi E (da Haridwar. An)eGangotri Lb p Br) 
E 12 gy rJ? Pje 


p+2,q+L py Qi s- 


La 


214 (my e) 
(es: JT : | J Yi Lp, 


| (27,5, jo Emu 


4 


provided that 


1. Re 29+23 ha, nee, 


i=l 


9: Rel 2p #2Y kA, +1)> 0 where À; = min | Refa /3%)] (i= WP ip) D 

V i=l 1 sysm, Y T 
Case II. Similarly if we put m=/R'=IR=0 in (4.1) we get the result in terms of 
multivariable H-function defined by Srivastava and Panda [8] as: 


cos0/2)" (sin/2)" HI (cos0/2) " (sin®/2)"" z,,....(cos0/2)” (sin0/2)" z, 
1 


= ah onim PUT NT 
v b lap: pan ep, qu] n 
T p*tqvip Qiii P, A, 


Ze 


Zi 


eu. 
(1/2 — p; A, ,...,h,), (cree jet 
(-p Pp, 3, +R,,...,h, + hk.) (b. Br BO) 


lg 


D up ltt) lj 
T 


(2,8;)... (80,59) ge um 


1g, 


e 
H” 2n magma, | : 
PL pin pd, | ^ 


(1/2+s=P=p3h,..h,),(1/2=8=Py5 Py». ki) 5 
(-P—p-p,;h, 1. h, tk): 


5 (r) D cup ) [5 
à. QU... o ) de 1 ( 6 2 
(a; jov Lp jo Yj apn C; Yj B 


(brBr B) y (27,85). (010,80) | 


Lg, 


„(6.2 


with the conditions as in (6.1) and P=0,1,2.... 
Case III. Again in (1.4) if we put Pp, =0, then the hyper-geometric function 25 


reduces to , #c-I'keruraskdtantoisctwritkerwasAn eGangot Initiative 


WAS ee ee 


Li 


l'(c)T(c- a - b) 1 


INEA) == 27 0) \ 
F (a,b,c;1) Eeo O 


f 
| 
li 
| 


| On further simplification, we get 


| cos PO(cos0/2)” qe = — T(2p«1) - 

0 ` 2"T(px P--1)* 16.3) 
provided that Re(p)>-1/2 and P=0,1,2,... which is known result given by 
MacRobert [3]. 

Case IV. Also if we put p=0 in (1.4) and proceed similarly as above, we obtain 


P(2p, +1)P(1/2+P) 


rms of Sd 5 2% Ja — i 
f cos PO(sin@/2) de = gn P EP (6.4) ij 
at A | 


provided that Re(p,)>-1/2 and P=0,1,2,... which is the known result due to 
Sneddon [7]. 


Case V. In (2.2) using r=2 and putting p¡=k,=k3=.....=k,=0. we obtain a result in 

terms of H-function of two variables (defined by Mittal and Gupta [4]), given by | 
koul (1,21. 

Case VI. Again using p,=h,=h3=.....=h,=k,=k3=...=k,=0 and taking r=2 in 


(2.2) we get the result given by Koul [1.2]. | 
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ABSTRACT 
All the 65 results of this paper from (2.1) to (2.65) are false in view of the 
fact that in right hand side of each results k's run over summation signs of multiple 
hypergeometric series and therefore no such k should be involved in the left hand 
side (i.e. left hand side of each result should be free from %'s.). 
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‘hag ABSTRACT 
In this paper, we consider a potential function to study the phase shift | 

difference of s-wave Schródinger equation. |! 
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65D20. | 
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1. Introduction. The phase shifts have great importance in the | 
computation of scattering and experimental work of nuclear and atomic collision. | 
Several mathematical problems are solved to find out a potential function from a | 
observed phase shift difference. Many workers namely Mahajan and Varma [8], | 
Raghuwansi and Sharma [9], Kumar, Chandel and Agrawal [5], Agrawal and Kumar | 
(1), Chandel and Kumar [3] have determined the phase shifts from a given potential | 


function. Recently, Kumar and Singh [6] have obtained an approximation formula 


for phase-shifts of s-wave Schrodinger equation with the application of binomial 


potential function and study the phase-shifts variation with respe 


parameter involving in binomial potential function. l 
d r a potential function of Mittag- Leffler 


thod [12], we obtain phase shift difference 


ct to the 


Here, in our work, we conside 
function and then on applying Tietz me 
in series form involving hypergeometric 
on shift variation verses the parameter 

The Mittag- Leffler function (see Er 

| Manocha [11]) is defined by 


function [8] and then make some studies 
s involving in the potential function. 
délyi et al. [4] and Srivastava and 


ci 7 (1.1) 
E, (2) = DY cioe arukul Kangri Collection, Haridwar. An eGangotri Initiative 
mol (an+1 


We consider the potential function in the form 


U(r)- E, (e), a,beC. (1,2) 
Particularly, on setting a=0 in (1.2), we find the potential function of Bhattacharjie | 


and Sudarshan [2] in the form 


U'(r)=U(r)-1= Ye" . zy 


nzl 


Also, from (1.3), we get 


O y Ae" at a=0. (14) 
The potential function considered in (1.2) have more parameters to that of 
Bhattacharjie-Sudarshan [2] function and Chandel-Kumar [3] function. In this 
paper, we determine the phase shift difference formula for this function on making 
use of the Tietz' s techniques and make some studies on shift variation with respect 
to the parameters involving in the parameters given in the potential function 
(1.2). 

2. Formulae Used. In our investigation, we need the applications of 
following formulae : 

For the s-wave radial Schrédinger equation 


aao ln L(L+1 
dr? *|K^-U(r)- e Vu 70. (2.1) 
Using Tietz [12] method and the Luke [7] formulae, Chandel and Kumar [8] have 
given the phase shifts in the form 
nk et -du ; L+2: 
Vo Qu pb m , TK ; 
¿Mia = gar Nan dr! ana i ja 
provided that (2L+3)>0 ...(2.2) 
The second phase difference formula is given by 
A -dU L+1; 
Mia Mia = z FE SUMA 2 | Kr 3) 
SO Ort ASA dr! Ta 3/2,2L 2; s je e 


3. Phase Shift Difference for Mittag-Leffler Function . From (2.0). 
we find 


d U be Bie 
dr = Gi rae Gurukul Kangri Collection, Haridwar. An eGangotri Initiative mee! | 
n=l 


Fu 


pr 


Now, making an appeal to (3.1) and (2.2), we derive 


112) In K* T(L+2) < 
., pe lin = aan l 27 Lyc L+2,L+2;-4K* | 
rarjie | ab P(L+3/2) Ta) N ana | 
provided that (2iK | nby)| si 
| 3,2) 
-(13) Further, Making an appeal to (3.1) and (3.2), we get 
Me Ny -MR T(Le1) 1 L+1,L+1;-4K? 
ab” T(L41/2) 4 A 2s A 
(14) $e ý 
provided that (2i /nb)" e x 
at of ADD) 
| this 4. Particular Cases . For L=0 , (3.2), gives 
aki ; 
Es. Noh = Unk Sic 1 E 2,2; AK 
aan ^" ab (3/2) (anys as SES 
while for L=1, (3.2) gives 
ns of 3 E : 
| mon =e NK eee 
1 , F = 3 es D 
: wr) baba (4.2) 
From (4.2), we may write 
AO) Si MAH - 
ies JnK 1 52.2) (-4K* /n%b*) 
1 y == ZO NL i Nl p A Je 
have PNG Aan 3(4), s! À p 
Now, adding (4.1) and (4.3) we can write 
5 n sm Ink e il 2,2; ARE 
2 abr (3/2) £& (ny T (an)* 1| 4. nij | ...(4.4) 


(2.2) 
If we put L=1 in the equation (3.3), then we get directly above result (4.4). 
5. Analysis. To analize above results we start with 
| Table No. 1 
r (2.3) 
(21), 22.9226033 


11.40945107 
3.800312649 
0.950005558 


; Haridwar. An eGangotri Initiative 


Di A do ah ar 


| 


Variation of no-nz with respect to a 


25 
20 
1 15 on 
= —-—Series1 
o 
TO, 
5 
0 
2-1 a=2 a=3 a=4 
a>, 
Figure 1. 


The above graph shows that when potential function Figune vo.1 woar i 
reduced on making increment of a so that internal kinetic energy of the particle i 
increased and by this kinetic energy the particle of s-orbit gains frequency of d 
orbit and then phase of s-wave tends to phase of d-wave (n, — 1 ), thus differenc 
in phase in sharply reducing. 

Table No. 2 


22.9226033 


0.239139705 
0.085412647 
0.027040942 


Variation no-no of with respect to b. 


No-na> 


—@— Seriesl 
10 


b=1 b=2 b= b=4 


b> 
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(Ls 
of t 
dift 


I1] 


[2] 


[3] 


woar it 
rticle i 
y of d 
ferent 


By second graph it is shown that when we make 
potential energy of the particle is grown up, 


the increment in b, the 
: as when as the potential energy is 
changed, the internal energy of the particle is also changed so that the particle is 


' disturbed due to slightly scattering and thus phase shifts is sharply decreased but 


when potential energy becomes large the difference in phase shifts becomes constant 
and the particle remains is s-orbit. 
Table No. 3 


3458.01349 


61372.55374 
465550.4809 
2233824.794 


Variation of n, -n, with respect to b. 
2500000 - 


2000000 ^ 


& 4500000 @b=1/2 
= 8b-1/3 
1000000 e 
500000 85-1/5 

0 ^C EP 


b-1/2 b=1/3 b-1/4 b=1/5 


b> 
Figure. 3 
From third graph, we have that when we reduce the value of bin the formula 
(1.2), the potential energy of the particle isreduced and thus internal kinetic energy 
of the particle is increased and by this effect the particle is scattered so that the 
difference in phase is sharply increasing. 
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ABSTRACT 

In this paper, authors derived some generating functions (partly bilateral 
and partly unilateral) involving exponential and Mittag-Leffer's functions in view 
of Exton's result [3]. 
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1. Introduction and Definition. The function 


m 


DR SELL c. (1.1) 
«(2) > ar) am 
was introduced by Mittag-Leffler's [5] and was investigated systimatically by 
several other authors (for detail, see (2, Chapter XVIII). 


The function 
k: 


3 2 
E.(z)= > = 72 w o, p » 0 (1.2) 
«el ) 2 rohs) 


has properties very similar to those of Mittag-Leffler's function E, (z) (See Wiman 


[9], Agarwal [1]). 
In 1971, Prabhakar [7] introduced the function E/,(z) in the form 


(Y), 


ES z 
Y No a,B,y>0 (1.3) 
Eis (2) 2 T (ok B) k! 


k50 


. is the Pochhammer symbol (Rainville (81) 


k 


(Y), = L(v), = y(y+1)(y +2)--(v+k-1). oe | 
The function E! y (2) is most natural generalization of the exponential 


function exp(z), Mittag-Leffler function E, 


Where (y) 


(z) and Wiman's function E, y (2). 


We note that.o Gurukul Kangri Collection, Haridwar. An eGangotri Initiative 


M 1 E 
E! (2)= Bap (2) Bun (2)=E.(2)-Es(2)= c ALP) 


By, (2)= En la) El) e E,(2)- coshz 
An interesting (partly bilateral and partly unilateral) generating function 
for F" (x), due to Exton (3, p.147(3)] is recalled here in the following (modified) 


from [see [6]): 


ex[s +t X Y, SAN ee), (1.5) 


mz- nam? 


where E" (x)= ,F,(-njm+1;x)/m!n! = L; (x)/(m * n)!, (1.6) 
and L” (x) denotes the classical Laguerre polynomials, (see [8] and in what folllows 


m* = max(0,—m),(meZ = 0,1,2,...) (1.7) 
so that all factorials in equation (1.5) have meaning. 

2. Generating Reltions. 
Result-1. If p,g,leN, then 


lurah Mr 
exp| s" +t’ — (E ) mn (7x) 
1 F 2 à" P r C E zt) (2.1) 


p q 


Special Cases 
(i) For p=q=l, equation (2.1) roe to 


D t a X st" -n/p a 5 
exp| s” +t? — (8) = a N ? Y 
| s DOCE mlp-*l ; e 22 | 


(ii) When p=2 in (2.2) or p=q=/=2 in equation (2.1), we have 


2 ct 2 te m ` ti —n / 2 . 
exp| s tE) j IS P 9 
| S Aa l m/2+1 F (2.3) 


D 


(ii) When p=q=1 in equation (2.1), we get 


i —n+j-1 
no L2 m L- 1 > 
exp s+t—(xt/s) ER yp 1 
| JE Eus 2 (m*1)*j-1, = (2.4) 
J=1 l 
For /= 1, egaxinotul Radgrrerdentes, tisriawg) Anteca Aei Gbtained from (2.1) bY | 


taking p=qg=l=1. 7 


14) Result-2. If p,g,le N and E/ » is defined by (1. 3), then | 
| i | 
1 pu P Tz ‘i Y i | 
Lion El (s' oU (4 E; Be (-xt/s) | 
led) 5 5 s" [" gi Jr mir (rs 2 ai «yr 3 Ya). (—x)’ 
= = La 
io det ALT BIB EI 70 mu (2.5) 
: 1 p } “G3 r | "tar 
L5) Special Cases. i 
(i) For y, = Y, 2 Y, =1, equation (2.5) reduces to 
1.6) SEN 
Ba Bsa) (2) 
lows s Y 
(1.7) 2i y S snp s p (x) 1 
z NEXT. (2.6) i 
mul «PB, TB, TB, r=0 (B, eun (B, ), {ct | (B. hi r | 
p Ad [| 
(ii) When p=q=! in equation (2.5), we get | 
) 4 —xt E | 
E, (s! )E, s, QE I 
(2.1) i ll 
sg" [n/p] (Ex)? | 


ne TB TBT B; 2 (8. Palen) lar : (2.7) 


For a, =a, =0, = B, =P, = D, =1 equation (2.7) reduces to (2.2) and also to 


(1.5) for p= 1. 
(2.2) © ii) If ¥,=Y. =Y3=1 and p=q=l=1 in equation (2.5), we get known 


generating function of Kamarujjama and Khursheed [4] 
—xt 


Ens, (s) Eu, (t)E., Bs ES 


(2.3) 


mz-en-m? 


Site (- -x) 
eS (2.8) 
= > mme Tp, TB, Bs DB, Jasmani (Bo Ie aUr-r) (Ba Jc » 
(iv) When p=2 in equation (2.7) or p- =q=l=2 in equation (2.6), we get 


A t 
E, s (s ) E, 8; (P). s. E 


ES st" (212) (=x x) : 29) 
hangotri Initiative 


2 X erna apaan ngan Maha 


a — 


For o, = 0, = 0 =P, =P, =P, =1 equation (2.9) reduces to (2.3). 


: : f di 
| Now DOE ay =O, 20, = 2,B, 2 B, - B, 71 in equation (2.9) and using a 
| relation (1.4), we nm a generating function of hypergeometric function gl. in 
| terms of hyperbolic cosine functions. ' 
ee e -n,-n-41/2 £ zt 1 
Pre cia cae 7 x 14 c 
coshs cosht cosh(xt/s) 2 2 am)! (an) Wan? Z a Punah . (2.10) 
(v) For p=q=1, equation (2.6) reduces to 
l 
i Sei) 
} E, 5, (s) E, s. (0). 9, Ga 
* Sur |n] (5x) 
-Y Y) : ; (2.11) 
Ea UE TB, FB.FD, LB. 2, (P, Ip TUE vin (Ba Me mE 
For œ, =a, = 0, = B, = B, = B, = 1, equation (2.11) reuces to (2.4) and also to 
(1.5) for l=1. 
Wi) For o,=0,=0,=1 and /=1, equation (2.11) reduces to a following 
generating relation ` 
A elen a 
Mz nm? »(B, ),, (B. se En tm ¿Ba > H p 
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ABSTRACT 

In the present paper, we introduce the subclass S, (n,p.4,g,0) of functions in 

(212 with negative coefficients, which are analytic and multivalent in the unit disc | 


| 
U= fz: lz|< 1}. The fractional calculus of functions associated with ingegral | 


236 | 


operator J., in the class S;(n,p,h,q,0.) as applications of the Saigo fractional | 
York, calculus operator J!" are established here. | 
5 : : y kown results for 
iting Corresponding to our main theorems some known and un 


the multivalent functions are also shown to be deduced as the special cases. 
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' ; i f the form | 
d'une 1. Introduction. Let T(j,p) be the class of function of the fo | 
| 

alized A ; | 2 * | 
f(z)=2? = Y az (a, 20;pjeN- (more | 

ction k=j+p | 


Jsea Which are analytic and p-valent in the open disc U =fz:|z|<1}. A function 
F 2 . H LA CE 
f (z)e Tj p) is said to be p-valently close to convex of order a in U if it satisfies 
fatha 


the i : P. 
| he Inequality CC-0. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative 


UO ——————————————————— — 


1 Rel z ^f (z)) » o. (zeU;0sa<p;peN). 9 

! A function f(z)e T (j, p) is said to be p-valently starlike of order o in Uif 

| 

| it satisfies the inequality | 
pl > (zeU;0sa«p;peN). 8 A 


Further more a function f(z)e Gs p) is said to p-valently convex of order a in U 


if it satisfies the inequality 


Re SUME) 20 TES peN w 
10) (zeU;0sa<p;peN). (4) 
W 
For each f(2)e T(j,p) we have [2] s 
jn = p 1 ¿14 ` k! ¿1-9 5 ie 
en Ai @eNe=NUOP>g). 6 
Sal ; I 
For a function f(z) in class T(j,p), we define the following differential operator: 
DY yp (z) = pa (z), E; 
D! f" ) pf )- zE [2 UC ) a 
z)- z) =- 
ia (p+A= q) Sc 
| | p! 5 e (k+A-q) R! = ? 
(2 = z^ OF — az" ] ; 
| (p -q)! np (D+ -q) (k-q)! D 
1-h 
Di le Je D Di, eu = Z fo 
[ (z )] cgi 2 »x (z 2)] 
— OU NS ee | k+à-q J i TI 
= BPS = EO aye 
(p-a)! en (R- q)! ptr-q 
and so on ` 


o EA o = aa] 


ptà-q 
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= | p! 2! k! ktk-q a 

1 . NUS ix pcm [AA ah 

EU (p - q): à Ten (R-q acd ‘ ! 
(P.J e N;qe Nu p>q;h20) (6) 


(3) A function f" (z) e T(j, p) is said to be in S, (n, p,A,q,«) if and only if 


j | Dif" (2) 
gni e [ral 20; (7 
l Dena) 1) 
(4) where z€U,A20,peN,q,neN,,0<u<p-qg,p>q and D’, is defined in (6). 
We note that for g-0 the operator D’, f(z) in view of multiplier transformation | 
is defined and studied recently by Agharaly et al. [1] and Singh et al. [13] for positive | 
coefficients in the following form : || 
(5) 5. ir 
SI (kh i | 
I, (n,X)f (z) 2 z^ + Ee az (2 0,n € Z). n 
ator: 2. p* || 


Earlier the operator I, (n,A) was investigated by Cho and Srivastava [6] and Cho 


and Kim [7]. Whereas the operator 1,(n,l) was studied by Uralegaddi and 
Somanatha [16]. 7,(n,0) is the well known Salagean derivative operator [11] 


defined as | 


D'f(z)- 2+5 k'a,z*,n e Ny = NULO). 


152 
2. Coefficient Estimate. 


Theorem 1. Let the function /( 


z) defined by (1) be in the class Tp). 


| 
Lg) i ly if 
Then the function f(z) belongs to the class S; (n, p.? ,q,0) if and only i | 


> ( kt+A~-q E E -a)é(p.q) 
A o O ENDE 
[E239 aaa 


h=j+ p 


(O<a< p-q;p,jeN;qne Noi p> A20) 
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Hn na I A AN Ba EIUS Naa Naga PA ES -~ 
—. Seas IE na 


: 2 - [p(p-1)..(p-a-*1).a #0 
Bag 1 ,q=0 (9) 


Proof. Let us suppose the inequality (8) holds true. Then in view of (7), we have ' 


k-X-q >|? 
aD aayi (z )} x mgl res 


| p^ 
| n (gl 
Da (2) R+A— -q i-p 
| Pr. ame | Es “a(n 
(p,q)= È E a Al 


-(p-q) 


5 E (k- p)5(%,a)a, 
h=j+p 


£ Ne 
AA PHA- POSU 


k+r- =) 
oup. Sala 
E din El q ( q) h 


Therefore the values of function 


z| D? f (z)] 


MO ge 


(10) 


lie in a circle which is centered at w=(p-q) and whose radius is (p -q — c). Hence 
the function f(z) satisfies the condition given in (7) 


Now conversely, assume that the function f(z) is in the class S j (n, p,A,q ,o). Then 


we have 
Ò 5S ktÀ-q | kg 
z| Dp. f" (z)] (p-a) (p,q)- ` (R=q) — 5(k,q) a,2 
Beng Re ee o um o >0 
IDEE 1 (z) xs peu 80 " EP (1) 
h=j+p ptr q ye 


for some 0.(0<0< p-g;p,jeN;q;jne Ny; p» q;2 0) and z e U . Choose value o 


z on the real axis so that ¢(z) given by (10) is real. Upon clearing the denominat? 


in (11) and letting z _, 1- through the real values we can see that 
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(p-q)8(p,a) SW ae) 812) 


ASE Do 


b ae 


.~(k+h-g)' il 
(9) 2018(p,g)- | ved SENE 1 
a | 2, DG jas, (12) j 


which leads to inequality (8). It completes the proof of Theorem 1. | 


Corollary 1. Let the function fz) defined by (1) be in the class S,(n,p,2,q,0) 


then the following inequality hold true 


e (p-q-0)ô(p,q) | 
k+r-q) 

EE ms (13 

(xs) e q ct) 5(z,q) ) 


(k2 j* p;p,je N;qjne N;A20; p» q). | 
The result is sharp for the function f(z) given by i| 


F(z) ==" - : AA ,9) z! | 
Lad (14) | 
— —Á k-g-a)s(k, 
(10) (aji q-0)5(%,g) 
Hence (k2 j+p;p,jeN;qneN,;k20;p>g). 
3. Applicatrions of Fractional Calculus. In our present investigation, 
Then We shall make use of the familiar integral operator J. „ defined by [5,p.676, eq. 
| (1.8) 
(TA) E t f(t) (f eT(j.p)ic» —p:peN). (15) 
Definition 1. The Reimann-Liowville fractional integral of order is defined, for 
>a 


(Ü. afunction f(z), by [15,p.224,eq.(3.D] 
De e CUBE ae (16) 
TO (ec m 


jnato! Where 150 f(z) is an analytic function in a simply-connected region of the z- 
inal , 


E k-1 - x m 
Plane containing the origin, and hte multip'icity of (z—¢) isremoved by requiring 


] 
98 (2 = C) to be reabh enis EI E Haridwar. An eGangotri Initiative 
a a : rivative of order A is defined 


eee e En — 


for a function /(z), by [15,p.224, eq. (3.2)] 


D'f(z)- >, (17) 


r(1 TN dz 
where 0€ X « 1,/ (z) is an analytic function in a simply-connected region of the z- 


plane containing the origin, and the multiplicity of (gear is removed as in 
Definition 1 above. 

' Definition 3. Under the hypotheses of Definition 2, the fractional derivative of 
order (n X) is defined by [15,p.225, eq. (3.3)] 


d" 


n) 
Dipa) 


—-D'f (z), (18) 


where 0<2<1 and n e N, = NU(0] . 


Srivastava, Saigo and Owa defined the following fractional integral operator 
involving Gauss's hypergeometric function: 


Definition-4. For real number «>0,8 and n, the Saigo fractional integral 


operator Ij!" is defined by [12, p.112, Eq. (8)] (See also [15]): 


7 

IN = cp NE AN a 1 : 9 | 

s Hg, z-C) oc B,-5;0;1- 2 |f (O)dCG. Gt) | 
where /(z) is an analytic function in a simply-connected region of the z-plane 
containing the origin with the order 

m olla), as 250, 

where €>max{0,B-n}-1 and the many-valuedness of (e-f is removed by 
requiring /og(z —&) to be real when (2-E)>0. : 
From Definition (1) and Definition (4), it is easy to see that 
D;"f (z) = Is!" (2). | 
Lemma. If œ > 0 and k>PB=n-1, then [12] F 


LP F (E 1)r (k- B+n+1) 


0,2 Lp 
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ive of 


(18) 


rator 


;egral 


(19) 


plane 


red by 


(20 


For a function f(z) € T(], p) defined by (1), we obtain easily hef , * 


llowing results 
in view of (15) and (20) respectively : 


c+k (c» -p;p.je N) (21) 

and 
giri)" ubt DG emen S T(Re)r(b-yemel) a 
F(p-y*l)F (p+P+n+1) e F(h- +1 \r(k+B+n+1)- ye 22) 


Now for the function (J. Fl) defined in (21) we obtain the following fractional 


integral in view of (22) 


MPAA 
gean Alf z) = A AAA 
E Il 3 ) Tip-y+1D)Fíp+B+n+1)" 


c (=) r(k+1)r(k-y+n+1) m 
eth JU(k-y+1)P(k+B+n+1) * 


Siam 
(B»0,c» -p;p,je N). 

Theorem 2. Let Qu and 1 satisfy inequalities B>0,y< ptl,y-n« p- 1l, 

B+n> -(p +1). Choose a positive integer such that n 2 y(B-n)/B- p-1. 


If f(z)e S, (n, p,À,q,c) then 


Dan r( #1)F(p-y#nt1) x (c+ p) 
e Ce Ik on ctptj 


Ppt j*1)T(p*j-Y*n*D(p- -g-o)8(p.a) la? Hal" 


: jtpt^-q ej p, 
l'(p* j— y -1)r (p+ jabensa [I T Jose q-@)5(j+ p,q) 


(28U,:0<0.< Pp=q;c> —p:p, je Niqn€ Ny p» 9:20) PA 


and 

r(p+1)r(p-y+n+1) | (c+ p) 
Ta Eu 
È- (gard Ires du: tie A A ER Initiative 


meester 


[ 
Ü 
D 
| 


meno mos DLL. 


Pp j*1)F(p*j-y*n*1)(p-a-0)89(p.q) 


i ' j*p*k-qY 
Mp+j=y+DT(p+j+B+n+1) : 


ptr-q 


(zeU,;0€a«p-q;c»-p;pjeN;q,ane Nu; p»q;A2 0), 
where 


U, y<p 
U, spas ws je e 


These results are sharp for function f(z) given by 


zZ 


Mp en oa. a 
eps) (J* p-a-o)8(J* p,q) 


(J, ee: Ke DP-c-eWbg  — Hp 


Proof. To prove the Theorem 2, we start from eq. (23), i.e. 


T A 
me erstelle Natan 


c (e) F(R- 3)F (k — y 4-1) E 
== A AA 
eel) O (R=y+1T(k+P+n0+1) ^ 


Now on setting 


pap ae M 
pete F(p*1)T(p-y-n41) aps 


TP) 


the above result given in (27) takes the following form 


Lo 


H(2)=2"- Y a,w(k)2" 


k=j+p 
where ; ; 
(1) = T(k+1)r(k-y+n+1) Mo-Y+1)r(p+P+n+1) (c+ p) 
Y P(k-y+1)r(k+B+n+1) I(p+1)l(p-y+n+1) (c+k) ` 


Since y(%) is aGteexensin Kafanevietpo Haide ARNI Initiative 


| (i+ p-a-a)8(i+ p.a) 


(25) 


(29) 


co 


(25) 


(27) 


(28) 


(29) 


xi »P,4,0,r. To Ie we give few cias 


y(k) € w(pt J) , 


F(*p*ll(-*p-y*n«1) l'(p-y«1)F 


ixi SE (p*B*n*l) (e+p) Hu 
Fü*p-Y*UTl(j*psBenel) T F(p1)F (p-y*n41) (c+ p j) e il 
Now in view of Theorem 1, we have 
a (P=4 -c 
2 SE = Pe) 5 
he j+p ^ Jtptk-q 4 O14 
9(J + p,q = +p-q- (31) 
(Fe " EU. Ju »-g-40) 


Hence for the function H(z) obtained in (29), we have 


IH (z) |z|” - || "^w p) Y a. 


h=j+p 


Therefore, making an appeal to (30) and (31), we obtain 


TERRENUM. M(j+p+1)P(j+p-y+n+1) Mo-Y+1)F(p+PB+n+1) 
"(c+ p+ iT rpo ra (Gep+Ben+1) Torir (p yni) In 


(p-q-a)é(p.q) |" 


> B i 


Now in view of (28), we atonce arrive at the desired result in (24). 
The second result of Theorem 2 is proved similarly in view of (30) and (31) for 


inequality 


I) s |o! e yo) Y a | 
h=J: 4 


p= jrp 


4. Special Cases. (1) If in differential operator Dj af " (2) defined in (6) we 


Consider ) =( then it reduces to the operator studied by Aouf [2]. Therefore, at 


A=0 all the results established in the Sections-2, 3 will provide the known results 


of Aouf [2]. 


(D The results aa in class S,(1,p,19, «) here in turn provide many ! 


known results studied in various subclasses by specializing the parameters 
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(i) S, (n, p,0,g, 0) 7 S (n, p,q) (Aouf [2]) , I 
Gi) S,(0,p,0,q,a)=8,(p.9,0) and S,(1, p,0,q,%) = C, (p,q,%) (Chen etal.[4)), | 


(iti) S,(n,1,0,0,a)=P(j,0,n)(je Nine Ny;0<0.<1) (Aouf and Srivastava 


13), 
(iv) S, (n,1,0,0, 0) = T(n,o)(n e N,:0 Sa « 1) (Hur and Oh (8)), 
à ( 
T’ (pa) (Owa[10]) 
S,(0,p,0,0,0) 2 4,” : 
(v) (0, p,0,0, 0) ee (Yamakawa[17])’ : 
Cip) (Owa[0] à 
1) , 
(vi) 5; (1, p,0,0,0) = 3 $ 
CT, (p,j) (Yamakawa[17)) [i 
(vii) S,(0,p,0,0,0)=T*(p,0) and S, (1, p,0,0,a)=C(p,a)(peN;0sa< p) 
(Owa [9] and Salagean et al. [11]), à 5 
(viii) S,(0,1,0,0,0)=T, (j) and $,(1,1,0,0,0)=C, (J)(n e N,:0 < œ < 1) [ 
(Srivastava et al. [14]), : Ik 
(ix)  S;(n,p,0,0,a)=S,(n,p,o)(p,j e N,ne N,¿0<0.< p) (Aouf [2)), 
lb 
(III) Ifin (24) and (25) we take y--P then these inequalities reduce to the 
following inequalities involving R-L fractional integral operator defined in (16) l 
Corollary 2. Let the function f(z) defined by (1) be in the class S; (n, p,,9,9):] 
Then l 
: (p+) E 
DANI. 2) 24 - 
not i 
[1 
(c + p)r(p + J ap 1)(p OE «)(8(p,q)) le ir (32) 
; 5 . AS n à 
(e+ ps DT (rejet 1523279 | Gs p-q-ejsG pa) |n 
ptX-q | 
(zeU;0sa« p-g;ß>0;c>-p;p,jeN;neNp>qià>0), [1 


and CC-0. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative 


T(p+P+1) | 
T 
14]), | it 
il 
tava AY ; : ; i 
(c+ p) (p+ j *1)(p-q- «)(8p,q)) E | 
; . j+pti-q)' Fi ne (33) 
C+ pe J) (poe ta Ser E IN Vela ms 
( )'(p*j«p A gue | (J+ p-q-o)8(j* p,q) | 
(z € U;0 << 1g — q;p > 0;c > —P;P,J € N;n € Nsp 2 q;^ 2 0) k 
Each of the assertion (32) and (33) is sharp for the function fiz) givey by (26). 
The results in (32) and (33) in turn at 1-0 give the known results due to Aouf [2, | 
p.32, eq. (6.6) and eq. (6.7)]. i 
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ABSTRACT 

In this paper we give a t-generalization of the new zeta function duc to 
Goyal and Laddha. This function is defined by series and the corresponding integral 
representations are established. We have also investigated the properties of the 
new type of generating functions such as integral representations and generating 
functions. Several interesting results obtained earlier by Goyal and Laddha [3], 
Katsurada [4] and Bin-Saad [1] follow special cases of our main findings. 
2010 Mathematics Subject Classification : 11M35, 11523. 
Keywords: t-generalized zeta function, integral representation, bionomial 
theorem, and gamma function. 

1. Introduction. An interesting definition of the zeta functions, due to 
Goyal and Laddha [3], is as follows: 


y e (1), Y 
ly z,a)= 5 —— ——,y|« 1, (1.1) 
| ) Zea] n! b 


T(A+n) 
Rea > 0,u 21, where (A), = FO 


This is expressed as the integral form 


» 1 Ea AN (1.2) 
ly za) ——[ £e (1-ye“) dt, 

) neh ( | 

This function is continued to a meromorphic function over the vise Sabang (cf. 
Erdelyi (21). Obviously, when p=1, (1.1) reduces to the zeta function studied by 
Erdélyi [2]. Katsurada [4] introduced two hypergeometric type generating 


functions of the Riemann zeta function as follows: 


so m 1.3 
«Ge Y tem) lle WP 


, 
m=0 meca. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative 


and 


m 


fivx)s > lv), Ket m) |x| «1 (1.4) 


m=0 
where v and z are arbitrary fixed complex parameters. 
Recently, Bin-Saad [1] has defined two new type of generating functions suggested 
by (1.3) and (1.4) as: 


m 


, h 

Way a)- 26st mar «1 (L5) 
4 , x" 

6 (55a) e (0), o(%2+ ma). «tls < lal (1.6) 


m=0 


where ọ is the generalized zeta function. , 

Its various properties are investigated including the integral 
representations, generating functions, partial sums and N-fractional calculus. In 
a recent paper the authors [5] have defined a new generalized zeta function and 
derived its hypergeometric types of generating functions. In this paper we aim at 
giving t-generalizations of the generalized zeta function.and at deriving their 
various properties and formulas including their integral representations, series 
and genrating functions. 

2. Definition and Integral Representations. A t-generalization of 
generalized zeta function is defined in terms of series representations as : 


d (91,2,a)= y A (2.1 
molattn) n! 
where || « 1jRea > 0, 21,te R,t» 0. 


It is interesting to note that for q= 1 , (2.1) reduces to the generalized zeta functio! 
studied by Goyal and Laddha [3]. 
Here we prove the following integral representation 


eg . = 1 SA Su yA 
0, EM e (1— ye" )" dt, (23 


where |y|<1,Rea > 0,1 21,te R,t » 0. 
Proof of (2.2). Let 


1 re -u 
=== ¿gal | BurukarK gridgollection, Haridwar. An eGangotri Initiative 
T(z) f, { y 


Or 


| r(z) n=0 n! 


este; On interchanging the order of i i 
a Integration and summation, we see that 


n 


SU 
Tea Y EA Ee ieee 
2 raj E "p !dp 
(1.5) 
(1.6) n=0 (a E mn) n! 
2, (552,2). 


egral This completes the proof of (2.2). 
us. In. We now ENG the t-generalization of two hypergeometric type generating functions 
nand Of the Riemann zeta function in the form: 


im al 2 E 
their @&(%3t) = Y C (zem), 
J m=0 m! (2.3) 
series ; 
| lx|« e, te R,t>0, 
ion of and 
1 > É os : x" 
a) Fuzz SONGA (2.4) 


I| « 1, e R,t » 0, 


Where v and z are arbitrary fixed complex parameters. 


m - s 
nctio 3. Integral Involving e (x; 1) and f (v;,x) . In this section we evaluate 


definite integrals involving the function es(x;1) and f? (v;t,x) in terms of other 


kind of zeta and hypergeometric functions. First, we recall the Eulerian integral 


(2 formula of first kind (cf. Srivastava and Manocha [6)), 


Pap at TEC pola) Oe) =O (3.1) 


| From the term by term integration, we can derive the following 
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: orem 1. Let Re(c-e das, n An eGangotri Initiative 


T (c) 1 h-1 DEL 1 yi ^ 
ES go pen e! dt = G!(b,¢;1,x) c 
re=o)r TOL (1-0) ef (atit) 3j 
and 
T (c) LON xul ART TAS 
I" ü-t fi (vit, xt" )dt = G; (v,b,c;t,%) | (3: 
PSI OA 33 
where 
> O) OS ut | 
Bae!) No) 1 (c+ tm) S Gen (3.4) 
and 
, dk (P vomer xt : 
SO M0) £ (c+ mo (o) mi (3.5) 


For t=1, (3.2) and (3.3) reduce to the known results given earlier by Katsurada 


[4]. 
Proof of (3.2). Denote for convenience the left-hand side of relation (3.2) by 1. 
Then in view of (2.3), it is easily seen that 


r(c) 1 
= 2a z+ um) nro br (^ 


Using (3.1), we obtain 


formel (1 x: 2i =b=1 di . 


E Pe x" x" T'(c)F (b wm)T (c - b) 
E PEN mi Pe -5)F DIT (e+ m) 
re) E F(b-- m) x" 


Pj 22r BERG. (zr un). 

After some simplification, we get the right-hand side of (3.2). 

Proof of (3.3). Denote for convenience the left-hand side of relation (3.3) bY J 
Then in view of (2.3), it is easily seen that 


T (e pel c-b-1 
J= Dm megah (1-1) dt À 


Upon using the Eulerian integral formula (3.1) and the definition (3.4), We a 
finally led to right- and Ride of formula (3,3 eGangotri Initiative i e | 


an 


l| 


rem 2. Let R A 
(33) Theo Rez » 0, Reu » 0 and Rek<1,teR,t>0 Then 


-l5-u 


1 = paie n-l,z -auyt P 
| mme p E C [nue "yi ajdudu 


(3.3) E " x , n 
É A EM Y 
= 2,0, =31,2,a ; (3.6) 


Tan (a * 1n) (a+tn) 
and 
(3.4) I r3 - 
re te 'G xe" ,y;t,z,a)dt 
n" = x y" r : 
7249 ya] v (3.7) 
E^ E € tn) (a+tm) 


urada, Where 


by I. C (x, y; t, z,a) - Y 9 (osz tama) >, (3.8) 
m=0 O 
lyl<1,Te Rt» 0,Re(z)>1,a + 0,-1,-2,... 
and 
vt co : x" 
6 (x, 937,2,@) = Y. (v), 6 (ps2 wma), (3.9) 


m=0 


|< 1,1 € R,T > 0,Re(z)> Ha #0,-1,—2,... , 


while $* is the t-generalized zeta function. 
For 1=1, (3.8) and (3.9) reduce to the hypergeometric type generating function 


Studied by Bin-Saad [1]. i i ) 
Proof of (3.6). Denote for convenience the left-hand side of equation (3.6) by L. 


by 4 Then, in view of (3.8), we have 


TAM sa x" y" 1 al “du 1 feta do 
muc ee | 


m.nz0 M Ya + ge) r'(u) » 


i : : ion, we get 
Applying the Eulerian integral for gamma functio g 
& 2. ul 
= n 2 x/(a T an)] 
L= NS : TCollection, Haridwar. An eGangotri Initiative 


n=0(Q + TN i m=0 mila ar Ue 


vo EE — 


After some simplification, we obtain the right-hand side of (3.6). 


Ji 


This completes the proof of (3.6). Proof of (3.7) can be developed on the 


same line. 
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ABSTRACT 

Our aim is to solve the quadratic Diophantine equations 
1161146329226x?-y2= 41. Two starting least positive integer solutions are 
obtained. 
2000 Mathematics Subject Classification : Primary 11GXX, 14GXX; 
Secondary 11D09. 
Keywords: Quadratic Diophantine equation, Quadratic irrational number, Period 
of Continued fraction, Convergent of continued fraction. 

1. Introduction. In Diophantine equations only positive integer solutions 
are calculated. Solutions of quadratic Diophantine y?+19=7" and y? 411-3" 


are discussed by Devi [1]. Unique integer solution of the equation x?-13=3" is 
discussed by Sharma, Singh and Harikishan [3]. Ternary cubic Diophantine 
equation is solved. Different patterns of non-zero positive integer solutions are 
obtained by Gopalan arid Somnath [2]. Our aim is to find the positive integer 
Solution of quadratic (in two variables) Diophantine equations. S ad 
2. Formation of the problem. Consider the quadratie Diophantine | 
equations Nx? — y? = x1, where N isa positive integer. Let /N = x, be the quadratic 


irrational number. The algorithm for JN generates the infinite simple continued 


fraction of the form [5.402.450] as follows | 
1 
1 

q Rl | 


= i 
In rebana ollection, Haridwar. An eGangotri Initiative i 


We do not have to calculate infinitely many q;s. Since the quadrati, 


irrational number is always periodic. So it is of the form 8 


KEREN he E It is also found that q,.,, = 2q,. Number of terms m, 


from q,4, to q, ,,, 15 called period of continued fraction. If the period is odd both 


[ 
the equations have positive integer solution. 


We denote the n!? convergent C, = P, /Q, , therefore [: 
y= == = o, [3 
Q 1 
C, o Ge 
Sae d (23) 
C = 2 = Cres + DES , for n> 9 
Q, q,Q,. ar Qi 


The solutions of the equations Nx? - y? =+1 are 
x2Q,yz-P,fori-0,2,4,.. and 
x-Q,y-P for 1=1,3,5,.... respectively. 

3. Results and Discussion. 


V1161146329226 = [1077565,2155130 |. 


In this case, period of continued fraction is one, which is odd. Therefore both the 
equations can be solved. 


From (2.1) and (2.2), we get 
j 
F, = 1077565,Q, = 1,P, = 2322292658451,Q, = 2155130, P, = 50048425770085811% 


Q, = 4644585316901, P, = 10786086382990825883438801,Q, = 10009688515401500726! 
Solutions of the equations 1161146329226x2-y?= 41 are 


x y 

1 1077565 

4644585316901 5004842577008581195 
and 

x y. 

2155130 2322292658451 


100096885154015007260 10786086382990825883438801 
Only two starting least. positive int 


úru ROSY ol teger solutions ara ntalculated. 


a 


= —— ”— 


h the 


11961 
7260 


of other Diophantine equation of the fo 


[1] 


[2] 


[3] 


rm Nx ~ y? = n, where n isa positive integer. 
REFERENCES y 
Usha Devi, On the Diophantine eugation x2 
Education, 39, 2 (2005), 84-86. 

M.A. Gopalan, 


+19=7" and x2+11=3". The Mathematics 
Me 1 * : 

anju Somnath, Ternary Cubic Diophantine equation 
2a-1 (2 2 al : a 
D(a? + y )=2", Acta Ciencia Indica, 34, 3 (2008), 1135-1137. 


PK. Sharma, M.P Singh and Harikishan, On the Diophantine eugation x^ -132 3" . 
Acta Ciencia Indica, 34, 3 (2008), 1043-1044. 


CC-0. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative | 


CC-0. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative 


to 


q-OPERATIONAL FORMULAE 
LNG T GEN ros 
q-LAGUERRE 
POLYNOMIALS 
By 
R.K. Kumbhat, R.K. Gupta and Madhu Chauhan 
Department of Mathematics and Statistics 
J.N.V. University, Jodhpur-342005, India 
E-Mail : drrkkumbhat@yahoo.com, drrkgupta@yahoo.co.in 


(Received : October 15, 2010: Revised : March 20, 2011) 


ABSTRACT 
In the present paper, certain q-operational formulae for the generalized q- 


polynomials J!4 (x;r; p,l) are developed and make an attempt to unify the various 


results. The results given earlier by Gould and Hopper [6], Singh and Srivastava 
(91, Al-Salam [1], Das [4], Carlitz [3], Joshi and Singhal [10] follow as special cases. 
2000 Mathematics Subject Classification : Primary 33D05: Secondary 33D45, 
33D50. : 
Keywords: g-operational formulae, g-Hermite and g-Laguerre polynomials, 
Fractional g-derivative, q-Gamma function. 
1. Introduction. Burchnall [2] made use of the opertional formula 

E OT da) 
to prove the well-known relation 


minfn,n) [m n 
H pan (x)= 2, (-2)' | jas. (x)H,., (x) : 0.2) 


Gould and Hopper [6] studied operational formulae associated with classical 
Polynomials and established that 


a A (xr, p) D", du 


Where the symbol Y is defined by 
Y =D- pra” +0/x% 
and satisfies the relation 


Y A 

| | A j i i i tri Initiative 
x" D" — (xD ct prx GPR Sul Kangri Collection, Haridwar. An eGangotri Initia' 

J=0 
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EZ KK 
Ta - 


PM 


1% 

and 

H’ (x,0,p)=(-1) x“exp(px )D" di exp(=px"); 4 
defines the elegant generalization of the Hermite polynomials to which it reduces i 


when @&-0,p-1,r- 2. 
The relation (1.3) provides a generalization of the formula of Burchnall [9] 
quoted above as well as of Carlitz's formula [3] 


n n ak 
A " QW x (ath) fa k 
[penser awe (DA aa 
for the Laguerre polynomials. 
Joshi and Singhal [10] also introduced certain operational formulae associated 


with a class of polynomials unifying the generalized Hermite and Laguerre 


polynomials. 
The Rodrigues formula (c.f. Joshi and Singhal [10]) is 

Ne, ey) =e atu an exp( px") D" jeu exp(-px")), ...(1.6) 
where 


ENE 
SE 


ni(2-1) 


c(l,n)= 


l being a non-negative integer. 

In a recent paper, the authors [8] established certain operational formulae 
for the generalized basic Laguerre polynomials with the help of known results. 
The object of the present paper is to define g-analogue of the generalized poly nomial 


Ji (x;r,p,l) and develop certain g-operational formulae for the generalizeli 


polynomials Jl" (x;r,p,1,q) and make an attempt to unify the various results. 


For real or complex 0,0 « |g| « 1, the g-shifted factorial is defined as 
( ) ( a ) 1 n=0 
aq), = (aa), = Ane M 11 
(1=q \(1-¢ (190) ne N > vol 


In terms of the q-gamma function, (1.7) can be expressed as 
P, (e n)(1- q) à 
FP, (a) 5 


. CC:0. Gurukul Kangri Collection, Haridwar. i Initiati 
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(eq), = so 08 


tł 


s the q-gamma function (c.f. Gasper iven DY 5 
.£. Gasper and Rahman [5]) given 


(aza) | : 


959). 


r, (a) a a a- li 
NM — (asa) (ays | 
luces] | 
B Indeed, it is easy to verify that | 
io) MEM __ (arsa) 
limT, (0) =P (a l T= 
Bar, (e) (9) gd T (1-gy Ky. (1,10) 
(1.5) where 
(a), =0(0+1)..[(0+n-1),n>1. LAL) 
iated : bees pre 
N. The fractional q-derivative of arbitrary order 450 for a function 
f (x) 2 x'*' , is given by 
son c 
1.6) Dan aS E ] 
( 1 ( ) IF (u-A) , ...(1.12) ! 
where 1 #0,-1,-2..... | 
For 2 =], the equation (1.12) reduces to i 
H 
r xi UE Ne i 
n, (e) UIS (ce 210213) il 
nulae a MI 1 i} 
sults. Further, we shall denote the infinite product | 
mials 


„lized peee ln Fe “gl 1.14) 


1-0 (1-4,9’)...(a-b.q') eccles 


į 
Its. in thi 
In what follows the other notations and symbols employed in this paper 
have their usual meanings. 


2. g-Extension of J^? (xir. pal). In this section, we define a q-extension of the 


SERIO aan 


Anne 


ou generalized porel. ya (xr, p.l) due to Joshi and Singhal [10], by means of 


1 


the following relation 


18 2 n(I-1X1-2) te (px”) NIS p 
3 JG pl) E l: E n; (o "ej (ps )- 271) 


D 
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AK 


SS ob 


where r,p,l and constants assume integral values, and the g-Leibnitz rule is 


n I 
n LA r Tu j al 

D (UV )- "| IDE (U) D .(I ), (22) ; 
where 

Y E. (q: q), 

r] (aa), (ara), (23) | 
By virtue of (1.12) and (1.10), we observe I 
ALVA ple (e. Oe 

lim( (sd) = 1 (xir. pal) - (2.4) 


3. The g-Operational Formulae. In this section, we prove the various q- 
operational formulae by working use of the q-differential operator ô= xD, , which Al 


possesses the following interesting properties: 


(i) F(6)[x" f (x:4)] -x*F(8- a) f (x:q), E ssl 

Gi) F(8)[e, (3). (4) = e, (g(x)) F(8-- xg") f (x:q), (30); TI 
Op 

Gi) x" F(6)F(6+a)..F(6+(n-1)a) =[x*F(8)]". (83) 


We now consider the expression 


| e, (o aes Di [x"e, (=px") y] 


i -e, ( DN [ae (- px" ) Y ] | 
4 
= ¢, (px) "6 (6—-1)(5=2)...(5—n+1) [xe (- p^) Y] [By (3.3)] Si 
a ja (d+a-+(I-1)n— pra + j) tay making an appeal to (3.1) and (3.21. 
We thus have 3 
prp (ov )r] =x Co Ine Eo 
yal : | wh 


CG30. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative : i 
On the other hand, replacing n by n-k, o by a+ Ik in (2.1), we find the ma 


n 
q 


p [ee (=px")Y | = a (Epa) 
12.2) | ie 


1] pena | n 
(-1) Horr 


5 í (ll) 
gui | Joi (xir, p,l,g). (35) 


(23) Comparision of (3.4) and (3.5), gives the q-operational formula | 


n 


Us +0+ (1 e" 1)n - pra" + J) £ EE iH | 


Jel h=0 k 


3 E 3 (a+th) if 
us q (56). a Ji (3r, plg) (3.6) Ig 
vhich Also ; 

REA ( px") D? [sta (=px")Y | f 
(3.1) 


= e, (px")8(8—1)(6 - 2)..(8— n1) a" "te (=pa")Y ]. 


32). Therefore, making an appeal to (3.1), (3.2), (3.3) and (3.5), we obtain the another q- 
operational formula 


.(8.3) x(6 Se 1) [nt (-px")¥ | 


ana 


D 
hs nift- n n-kyt-yt-2y2 | Y | y. x! 
SU Mer (-px") (-1) l; A z 2 


k=0 


I 
i 
SI anna ba) ary (3.7) | 
Next we observe that the g-analogue of the recurrence relation of Joshi and 

Singhal [10]: | 


c(l,n) tud f 
(x'D, tax = pa "It (zr, pla) = coal) rm (xir. p,lq) 


Suggests the g-operational formula 


c(l,n) {a-nl) f 4. l 
meet =e queues pg) ...(3.8) 
Bi (xr, p, q) c(lm n) ( 


where 
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[4 c 
Y =D, rax 


Ko 


Ta 
which corresponds to the q-analogue of the formula of Gould and Hopper [6] to 


which it reduces when /=0. 
4. Particular Cases. In this section we shall deduce some interesting (i 

particular cases of the q-operational formulae. | 

(i) If we let g—1 and make use of the limit formula (1.10), in (3.8), we obtain 
known results due to Gould and Hopper [6]. 

(ii) Ifwe put /=0 and make use of (1.10), our formula (3.6) reduces to (1.3) due to 
Gould and Hopper [6]. 

(iii) Again, if we set p=l=r= 1 in (3.6) and make use of the limit formula (1.10), we 
arrive at the known results due to Carlitz [3]. pi 


(iv) On the other hand, when /=1, q > 1 and use of (1.10), (3.7) yields the known 


formula due to Joshi [7]. 
(v) “If we put 121, p=r=1 k=0, Y=1 in (3.7), we get the q-operational formula of 
Das [4]: 


(20,0) {x"e, (7) = 2e (-2)(a54), 12 (50). 
(vi) Ifwe take /=k=0, Y=1 in (3.7), we obtain 

{x(xDq + yr eta (-px' ) = (-1)’ er H" (xr, p,q). 

5. Applications. Setting Y=1 in (3.6), we have 


(i 


n (1-1)n 


Q 6+¢at(l—1)n—-— - Pa. A ae ja : 3 
(i) IT (1—-1)n— prx J) e(ln) " (sr, p,1,9), x 
so that 

x conn) 5 | 

min eon na) 

- HG s(- m en)- pra’ + j) wi 

jal 
(v. 


li 


E (8o (1-3) (m n)- prx" + j+n) 


j=1 


| (6 *ac (1-1) (m n)- pre +j) 


j=1 


xn m 
= (0 |+ a +n + (11) - prë" * yer C eerte 
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Therefore in views of (3.6) we finally derive 


fo: 


T RT 
c(Lm)e(t,n) jj i 


tin Gi) c(l,m+n) men (257, p,lg) 
tam A 1 Y (n m a 1 
(ig) an k20 k 5). s 
ue to (pelan ares E ras(I-1)m 
Jae” (BP, Bolg) Di J," (xr, pia). (5.2) 
), We If we reverse the order of the operator on the left hand side in (3.6) and 
proceed as above, we obtain an alternative q-operational formula 
10wn c(l, m)c(1, n) 
P: > , J” 3 l P 
(iii) c(l, m * n) men (xr, p, ,9) 
ula of 
b 1 Sapa m x! i 
Qoam h=0 k diaz ze 
Jim (aire plg) Di JC (ar, p,q). (5.3) 
Comparision of (5.2) and (5.3) gives the identity 
1 k 
, c eya) | M fra, x 
(iv) 2 (-1) E [eee | 
k=0 
(5.1) 


Jeruk) (xir, plg) DEJE 0" (rpg) 


m-k 


m 


i k 
«(1-1y2-2)/2 | mM x 
, = XC) 1)(2 aeaa gem] 


h=0 
JEM Gcr. p,a) Dida” (x57, p.a); 154) 
which in the special case /=0, making use of (4.6), gives 
A A 
+J) k=0 
=) ai r (seen na) Di H; (o " na) 55) 
k=0 a 


for the generalized -Hermite polynomials. 
- Next in oa eye Kenori Colgction, HadwasithshyimgiOmidhoth the sides by 


UU —_— —É—_É ———————————— — 


p > 
TX and summing from m=0 to m=, we arrive at 
(14) n 
> (Lg 
(vi) ( lx, acm (x;r, p,l,q) 
EE 
2 GG tan T 
z (Lia) Jy” (x Y Atx",r,p,l,q) 2 TA | Qr p,hq) 
mzü , m 
where 
A, zc (A (oy te 6.8) 
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ABSTRACT : 
5.6) In this paper discrete Hermite transformation of generalized functions 


belonging to a certain testing function space has been defined and an inversion 
formula has been established. 


tin 2010 Mathematics Subject Classification : 44A20, 46F10 
zo Keywords: Hermite Transformation, Generalized Functions. 
1. Introduction. The aim of the present work is to extend the Hermite 

transform 1 
ials, SEA 

F(n)= H {f (x)}= [e H, (x)f (x)dx (1.1) I 
^1, and its inversion formula il 
v fi) YH, (lees) a» — || 
cad. nao Misia: ; eee : || 
275 Mis died by Debnath [1], to a class of generalized functions. Zemanian [3,4] has also I 
ress, Studied (1.1), but here the present work is treated differently. 1 


2. The Notation and Terminology. Throughout this work the interval i 


ions; (-co,co) is denoted by 1. DKI) denotes the space of infinitely differentiable functions t 
! onI which have compact support in I. The topology of D(I) is that which makes its | | 

jals, | dual DD the space of Schwartz distributions. E(D denotes the space of all infinitely | 

differentiable functions on J. Its dual E'(I) is the space of distributions with compact - | 

ted | 
i Support. The symobls Os denotes the operator | 

peat Sio (De De" ) - (D? «22D, +2) , D, e d/dx (2.1) | 
3. Testing Function Space H(I) and Its Dual HI). We define H(I) as ! 

ert? the space of all complex valued infinitely differentiable functions (x) on 7 such | 

982) that | 

3.1 | 

P, (9) = sup loe(x)| ke , HOLM gi ¢ (3.1) | | 

ALI K0 1 I Y " * ; | 

| It is seen that H(D is ee space E be topology óf N a aehd "iterated by the | 


—- MEE — 


countable multinorm [p,], 4 = 0,1,2,3,.... - 


BD is the dual of HU) and is equipped with the usual weak topology. The members M 
of H(I) are called generalized functions. It can be seen than H(I) is a complete | 
if 


countably multinormed space. 


; : : f 
We now give below some properties of the space H(I) and it dual A’(/). 9 


(i) D(1)< H(I) and topology of D(D is stronger than the topology induced onit | i.e 
by H(I). Hence, the restriction of any member of H'(I) to D(Z) isin D'(I), 

(ii) D(I)c H(I)c E(I). As D(I) is densein E(I). H (1) is also a dense subspace 
of E(I). Consequently, E'(I) can be identified as a subspace of H'(). 


(iii) For each feHI), there exists a non-negative integar r and a positive We 


pr 
constant C such that for all y e H (I). T 
I f,$»| € C maxp, (0) 83. ol 
y Osksr OI: 
The proof of this statement follows from the boundedness property of 
generalized functions. 
(iv) If f(x) be a function of x defined in the interval (-e,ee) such that 
Le 
Ma (x)|dx <, then f(x) generates a regular generalized function of HI) 
Le 
defined by > 
(-R 
<f9> |" f(x)0(x)dx,o e HI), (33 
) 
This result can be easily established as follows | 
ES for 
ln f,0>1< p, (6) | If (a) < co, Th 
[ 


(v) For each positive integer n, the function 


La H, (+) (= <X< 00) is a membe defi 


of H(I). This can be proved as follows : 
9. [e7n, (x)|= -2n{eH, (x)} 
Hence, | 
vH = sup [2 fe le" H, (x (-2n) LH i 
p. le (x |= cen Gurukùl Kenek (9 gU ar. Moet Ini (ls vz, a= Le 1, 2 25 | 


sitive 


mbe 


4. The Discrete Hermite Tran 


sform of G i Tang 
Members of H'(I) are called Discrete Hermi Pm 


te transformable generalized functions. 


The Discrete Hermite transform of fe H'(I) is defined as an application 


of f e H'(I) to the kernel e y (x)e H(1), 


n 


ie. F(n)2 H(f)s (F(x),e" 8, (x)),n SO (4.1) 


e I H,(x)H, (t) 


As we require time to write the expression Y PRA 


Te 2" n! : 
therefore again, we denote it by T,(t,x);N being any positive integer and 
TOLAK, o<t<o, 


We now state below several lemmas whose proofs are based upon the similar lemmas 
proved in Dube (2). 


Lemma 4.1. Let fe H'(I), then for any positive integer N and for any arbitrary 
(x) € D(Z) 
de J -U =<" 4 = -ü ^ q Aly 
PA), (52) e(x)dx -(ro.f e Ty (t,x))e $(x)dx (4.2) 
Lemma 4.2. lim [^ Un (t,x)e dx =i, 
Lemma 4.3. Let (x) be an arbitrary member of D(I) with compact support in 


(-R,R) el, then 
R dici 
Lus (t,x)[ (x) - e(£)]e dx>0as N30. 


5. Inversion Theorem. We now prove the following inversion theorem? 
for our Discrete Hermite transform: 


Theorem 5.1. If F(n) denotes the Discrete Hermite transform of f(t)H'(I) as 


defined in (4.1), then in the sense of convergence in D'(I J^ 


> 1 H, (x) t),o(t)>as No. (5.1) 
(Sa 2" n! a ael) Kasigri 30529 El An eGangotri Initiative 


vr y 
Now $(x)e D(I) es e ola) e D(I), hence (5.1) is equivalent prove 
4 y ) 
€ (x t),e* $(n)» as N40. (5.2 
mo 2'n! dl )> «(e vo) | 
We have 
O) 63) 
“avn 2n! 
R x 1 H, (x) =a” va 
d ex 2^n! F(o)e deos go * su 
ac 
[As (x) € D(I), the support of (x) is contained in (-R,R)) be 
20 
R N 1 H (x) -x -x* K 
= eat t,e* H, (t))e x)dx (5.5) 
TX ga V e" B, Oe“ (x) E 
R pus zu) D 
E DEBER e ax) dx (5.6) De 
A ) e 2 9! n! ( ) T} 
5 | lin 
= J (FET, (63) (Jas EN sa 
9n. 
= (F(e)e" TE (t.x)e"6(x)dx) (58. fol 
y cor 
> (F(t),e°6(t)) as >, Qu an 
» cor 
H, (2), ES 
(5.3) equals to (5.4) as the function xc E (n) is locally integrable over t 
nz xn 2 n! 5 Th 
interval I and $(x) is in D(I) with support in (-R,R). E 
; ; y 
(5.6) is obvious because of linearity property of functionals. From Lemma 4-1 y 
obtain (5.8). Finally, Lemma (4.3) helps us to derive (5.9). in { 
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ABSTRÁCT 
In order to improve the efficiency of the 2N-Ray algorithm; we propose a 
variant of the Dps-Triangulation. A nice property of this triangulation is that it 


*subdivides all the subsets, on which the 2N -Ray algorithm works, into simplices 


according to the D Rs”Triangulation. Numerical tests shows that 2N -Ray algorithm 
based on D; ,-Triangulation is much more efficient. 
2010 Mathematics Subject Classification : 55M20, 54H25; Secondary 57915 
Keywords: Artificial algorithm, triangulation, admissible lebelling, nodes, 
barycentre, extra layer. 

1. Introduction. Various triangulations have been situdied by Dang [1], 
Dang and Talman [2], Scarf [4], Todd [6,7], Vander and Talman [8], Vertgeim [9]. 
The 2”-ray algorithm was proposed by Wright in [10] to compute solutions of non- 
linear equations. The 2”-ray algorithm partitions R” into 2n cones which have the 
same vertex. Then a triangulation of R” subdivides each cone into simplices. The 
2"-ray algorithm starts at the vertex and leaves it along an edge of some cone. It 
follows a sequence of adjacent simplices with varying dimension. Under some mild 
conditions, the 2”-ray algorithm terminates at an n-dimensional simplex that yields 
an approximate solutions to the system of non-linear equations. Since heg 2n 
cones have 2” edges, each of which is a ray, the 2"-ray algorithm has 2” possible 
tay to leave the vertex. ' , 

Motivated by this work, we have tried to cadar a new Dps triangulation: 
The D rs-triangulation based on the 2” ray algorithm is much Hone pEr el 
2" ray al gorithm proposed by Wright [10] as well as K, triangulation introduce 


by Kuhn in [3] and J, triangulation given by Todd [5]. i 
sane ed Definitions. The following notations have been used 


in this paper : 

R: Set of real numbers, 

Z: Set of all integers, 

N: Set of positive integers (1,2,...7), 


NO. Set of all non FepetivsuegerSoNctak Yaridwar. An eGangotri Initiative 


- 


R": n dimensionally space, having co-ordinates indexed 1 through n, 
? > i 
R^*!: n+1 dimensional space, with coordinates indexed 0 through n, 
Tt: Group of permutation on (1,2,...,n) and 1+1 group of permutation on 
(0,1,2,...,72), à ! 
With unit vector in R™ jeN and u= a up 


R?: Non negative orthant of R"i.e(xeR";x 2 0). 


Now we consider some standard definitions and explanations which will ho 
used in this paper. t 
2.1 Standard Simplex. The standard n dimensional closed simplex Sr is 


the convex hull of vv! ,... v" i.e. S" = {xeR!""! RL s” denotes the face of y 
opposite v'i.e.S” = { xes” DES 0) and boundary of s" is denoted by ds = UNST 
Again a j-dimensional simplex or [j-simplex] is the relative interior of the 


convex hull of j+ 1 affinely independent points y, y!, y? y?..... y’, called its vertices, cc 


We write o= yen . A simplex t is a face of c if its vertices are a subset 


vertices of the c. It is convemient to call the closure of a (j-1) dimensional face of 
the j simplex o as a facet of c. Two į simplices are said to be adjacent if they share 
a common facet. 

2.2 Triangulation. A trangulation G of S" is a collection of n simplices 
and satisfies the following two conditions : cedi 
ill The simplices in G together with all their faces form a partition of S" and 
2. Each point of S” has a neighbourhood meeting only a finite number o 
simplices. 


2(a) Pivot Rule. For a given simplices G and a vertex y of o the rules for obtaininé 


the simplex of G whose vertices include all vertices of o except y, are called the 
' pivot rules of G. 


ill 


2(b) Mesh. The mesh of a tringulation G is sup... diamo. We shall use tH 
Euclidian norm through out this paper. 


2.3 Definition. For each sign vector seg", let 
E(s)={ xeR" : six; = xl whenever s, #0) 2.3 


= cone { eR" :¢is a sign vector, and s, # 0 > s, m. 


' 
In case s has k non-zero components for k>0 than E(s) is a polyhedral cone 1 


dimension n — hac. Alsewveda eiO) Harmmanyig FEE Wah», each E(SINA 


f the 
tices, | 
ubset 


ice of 
share 


olices 


1 and 
er oli 
ning 
] the 


> the 


D 
f 


re | 


nr 


isa polyhedral of a cubical subdivisions of B"- where p^ donote the unit ball in 


1^ norm. 
Wright [10] has also defined another subdivision of R^ 


into closed convex 
cone as n-dimensional geometric form given as follows : 


x, -0 ifs, =0 


LetC(s)- xeR" : e 2 0ifs, # 0 


=cone ís;u; :s0, for each sign vector s). If s has k non-zero components 
then C(s) is an orthant of a k-dimensional coordinate subspace of R^. 


Wright [10] has proposed two type of T-triangulations of R” with the property 
that E(t) is a subcomplex for every sign vector seR" for (40. The first 
triangulation is called a K! triangulation. This triangulation is obtained by taking 
the triangulation K, due to Kuhn [3] in the first orhant and the reflecting through 
coordinate hyperplane to triangulate the other orhant. A vector v! of a n-simplex 
<u!,....u"tl> of K! specified by choosing a sign vector s with all the non-zero 
components, is a member of C(s) when all its components are integrals and x is a 
permutation of (1,2,...,n), then v!eK! is defined recursively as : 


y> =y’ NU for i-1,2,...n. 


For n=2 triangulation K! can be illustrated by following diagram. 


The second triangulation J, is defined by Todd [5] and whic.. can be 
illustrated for n=2 by the following figure : 


X2 (central vertices are heavy dots) 


1 
25, The D ns Triangulation. Define x 
Ww [xe :%, = maxx,,i = 23 a} 


taking a vectoreY SAGs ganar GpliFotiony etait. An eGangotri Initiative 


IEE 


yi HE aa | 


[x, |+1 otherwise | 2 
where [a] is the greatest integer less than or equal to a. Let D is the set of all | 
) 
Yew" where Y, is defined above. If YeD, we define | 
T | 
I(y)= fieN;y, = 9.) and J(y)- [jeN : y, > 9,) $ I f 
If 
Let son, s) be a sign vector such that ; 
( 
1. For ¡eN , if y, =0 then s, =1, and if y, #0 then s, =-1. 
Let K(y,s)= fiel (y) gom 1} ; 
Let / denote the number of element in Ily) and A the number of elements in 
K(y,s), we take integer p such that 
1. when h=0, if ( =n then p=0 or 2, 
2. when 70, if h=n then p=0 and if «n then 0€ pEn-1. 
Let m= [n(1),n(2),..,n(n)] be permutation of N. 
When h=0, for ¡=1,2,...,n, | n 
1. If 7=1, define TI 
-an |—Lifiel(y) 
(ns (1) 

&) 6 otherwise : ; 6 
of 
: for i=1,2,...,n. th: 

2. If 7 #1, we define y D 
a |S ifi=j 

§i\J) = i (2 for 

l G) l otherwise 
] | for 1=1,2,...,n. E 
| When h>0, for j=1,2....,n. 17 
1. If x(j)eK(y,8), define Th 
Be) ai EUR a 
0 otherwise | (a) 


for ¿=1,2,...,n, 
i 


CC-0. Gurukul Kangri Collection, Haridwar. An eGangotri Initiative 


f al] 


———  — —— 


ts in 


2. If n(j) € K (y,s), define 


s, (J) ifi 2 n(j 
g,(n(J))= ge!) 
0 otherwise (4) 
for i= 1,2,...,N. 
If y, x, s and p be as above, then vectors Y^, ,..., y” are defined as follows : 


for p=0, we have 

y =y 

y = y* &(n(k)),k -1,2,..,n 
and for p21, we define 


y =y+s, 


y! = y?" - sU ORR =1,2,...,p=1. 


y =y+8(m(k)),k= p,n. = 
The y”, y*,..., y” vectors obtained from the above definition are affinely independent. 


Thus their convex hull is a simplex. Let us denote this simplex by Dp, (y,m,s, p) or 


yt. y") . Let Ds be the set of all such simplices. Then Dpg is a triangulation 


of W”. Note that simplices of the Dpg-triangulation can be represented in more 
than one way. Moreover, the triangulation of a whole cube in W” is the same as the 
D -triangulation. z 

To be more expatiate let us illustrate Dp,s-triangulation of W” 1 »=2 and 


for x<4, Obviously, we have that for y, <4. 

D={(0,0,0)" (2,2,0)" (4, 4,0)" (40,4) (444) |. AD) 
“Let y= (0,0,0) Then, Z(y] Hie 9, = 9) = (128) end es 

Then s must be (1,1,1)7. Thus 

K(y,s)- lig (y) :s, =1}=(1,2,3) and h-3. We have p=0. 


9) Let m= (2,81). Then 77 =(3,12) and by applying (8) 


a(x (Uje 2)= PIC Le y a (2 An eGangotri Initiative 


) Be 2), 83 


Therefore, y’=y= (0,0,0)" , 


y? 2y*a(n(2) - (0,1), (b) 
y =y+8(m(3) - (0,0) , Let o' =(y°,y',9",9°)- 
(b) Let x — (3,2,1). Then x7 =(3,2,1) 
: g(n(1))=a(3)=(L41), | 
g(n(2))=8(2)=(1,1,0)", pe 
g(x(8)) = g(1)=(1,0,0), 


Therefore, yy =y= (0,0,0)" , 


+8(m(3))=(1,0,0)”, Let o! =(y",y',y, 9°). “Wet 


2: Let y-(2,2,0)' . Since I(y)={ieN : y, = y] - {1,2} and ¢ =2. 


Lt the d 


So s must be (-1,-1,1)7. Thus h(y,s) = fiel (y);s, = 1} = and h=0 while p can bea) 
one of 0,1,2,. | ES 
N ow by applying (1) and (2), we have 

g(1)=(-1,-1,0) 

£(2)=(0,-1,0) — | 


g(3) = (0:83) Gurukul Kangri Collection, Haridwar. An eGangotri Initiative A 


MEMO RA —— 


— 


| Now considering different values of p and applying (6) and (7), we obtain the 


following simplices: 


| (a) For p=0. Let 1 -(1,2,3). Therefore, 
) T 
| y” =y = (2,2,0) , 

y =y+8(n(1))= 


(1)) =(1,1,0)", 
| y! = y*&(n(2)) - (211,0), 
(3)) = (2 1)". Let o* AP ap). 


(b For p=1. Let x -(1,2,3). Therefore 


y =y+8(n 


| Let o^ = (^y^ y^). 


| 3. Let y= (4,4,0)" . Therefore, I (y) ={ieN;y, =y,}={1,2) and (22 


| We have that s rust be (-1,-1,1)7. Thus K (y,s)- fiel (y): 5, =1)=0 andh=0. 


| We have that p can be any one of 0,1,2. We also have 
&(1) = (-1,-1,0), 
g(2) = (0,-1,0), 
any | 8(3) = (0,0,1). 


| (a) For P=0, and = (1,2,3), we have 
Y = y=(4,4,0)" > 


Y s y*g(n( (3,8,0)" 


| (1))= 
| P sell (2)) =(4,3,0)" 
| (3)) 


y =y+ 8 (m( ES È qae AP Prey grin eGangotri Initiative 


5. 00 E 
(b  Forp-l, let 7=(1,2,3). Therefore | 
y = y+s=(3,3,1) , 
y = y+8(n(1))=(3,3,0)", 
y = y+ g(n(2))=(4,3,0)", | 
* = y+g(n(3))=(4,4,1) . | 


J Leto’ =(y",9', 979"). 
(c) For p=2 and r = (1,2,3), we have 

y! = y+s=(3,3,1)', 

y 2y Su z (441), 

x = y+8(m(2)) - (43,0), 


y = y*a(n(3) - (44,1) Let o° (9, y? 9°). 


) (a) 1 
4. Let y=(4,0,4)' . Since I(y) = {ieN; y, = 3] * (,3) and ¢ = 2. so that s mu 


be (-1,-1,1)7. Thus K(y,s) = fiel (y): s; = 1} = and h=0. we have that p can be an 
one of 0,1,2. We also have | 

g(1) = (-1,0,-1), 

g(2) = (0,1,0), 

g(3) = (0,0,-1). © 
(a) For p=0, and n= (1,2,3), we have 


y'2y-(40,4) 


— 


3). Leto? =(9,9!,9, 9) se 
(b  Forp=1. Let 7 =(1,2,3). Therefore 
y =y+s=(8,1,3), 
y =y+8(n(1))=(3,0,8)”, 
y+ a(n(2))=(4,1,4)", 
Y 2y*&(n(3))- (4,0,3)" Tet o (yl y! yt : 
(e) For p=2 and x 


| CC-0. d = a25 8) Geha Rewar. An eGangotri Initiative 


| 5. 


E su E (AMS ES 


y 
y'=y 

y 2y*&(n(2)- (41,4), 

y" = y+ &(x(3)) = (40,3) . Let o" =(y", yt y?, y?) 


Let y =(4,4,4)’ . Since I(y) = {ieN; y, = y,) =(1,2,8) and ( = 3. We have that 


a must be (-1,-1,-1)7. Thus K(y,s) - and h=0. we have that p.can be an A 


| 0,1,2. We also have 


g(1) = (-1,-1,-1), 
g(2) = (0,-1,0), 
g(3) 2 (0,0,-1). 


| (a) For p=0. Let x = (1,2,3). Therefore 


S mus 


be ani 


= 


Q 
- 


y =y=(4,4,4), 

y = y+ g(n(1)) = (8,3,3) 

Y = y+ g(n(2)).=(4,3,4)" 

y =y+8(1(3)) (44,3) Let o? = (y^, y y^ y?). 


For p=1. Let n= (1,2,3). Therefore 


For p=2. Let x= (1,2,3). We have, 

Y = y+s=(3,3,3)', 

Y = 99 -sgu = (43,3) 

y= y+ e(n(2))=(4,3,4)", 

y = y+8(m(3)) =(4,4,3)". Let o" =(y°,9',9,9”) 


It can be seen that le 2p = Ao 414] form a triangulation of W? for x « 4 


| 8$ sho 
wo by the following. Bigume kangri Collection, Haridwar. An eGangotri Initiative 
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ABSTRACT 

In this paper, making an appeal to computer we try to provide the measure 
of efficiency for triangulations through graphical represenation. hi 
| 2010 Mathematics Subject Classification : 55M20, 54H25; Secondary 57915 4 
Keywords: A.D.D., Triangulations, Fixed Points, Complementary pivoting 
techniques, Leveling of vertices of the triangulation, Homotopy, Simplex, Sandwich, 
Finer meshes, Directional density, Theoretical measures of efficiency. 

1. Introduction. There are number of algorithms for comuting fixed points 
using triangulations and complementary pivoting techniques. Each algorithm 
requires a triangulation, a lebeling of the vertices of the trinagulation, and a 5 
particular starting point called as pivoting point. | 

The successful algorithms use one of the two methods to obtain better and " 
V better approximations of fixed points, without wasting previous approximations. 
The first one is the restart method of Merril [4], independently developed 
in|  byKuhnand Mackinnon [3]. called "sandwich" method. If a fixed point of a mapping 


f on R” is sought, one triangulates. RR” x[0,1] rather than R”. All vertices lie in 


` either R” x(0) or R" x{1}. The latter are labelled according to the given mapping 


| 

f the former according to the constant map taking each x into CeR". Tue point c 
760) — can be so chosen as to be closed to a previous fixed point and provides a starting Ñ 
Point for the algorithm. | 
The second method was introduced by Eaves [1] and extended by Eaves and | 

| 


os| Saigal [2]. Instead of triangulating S" or R”, one triangulates S'x[L=) 


m ER X[1,»). Roughly, triangulations of finer and finer meshes are placed on the 


top of the another so that boundary simplices match. Then one generates an infinite 
Path of simplices following fixed point of piecewise linear maps that approximates 

better and better as the artificial coordinate tends to c. This process of 
deformation correspdttds $t E "no hobtdpy Heótueept. “Singgii weialife not concerned 


with the maping f we shall say such triangulations have continuous refinement of 
grid gize. f 

The efficiency of fixed point algorithm is very sensitive to the triangulation 
used. In this paper we are concerned with comparison of the theoretical measures > 
of efficiency for different types of triangulations. The technique of the comparison 
is based on the concept of directional density of the triangulation in the direction R 
das defined by Todd [8]. The only other measures known are rather crude measure iL 
based on the number of simplices in the unit cube and the diameter, introduced by 2 
Saigal, Solow and Wolsey see [7]. 

In this paper, we have also tried to provide the measure of efficiency for 
such triangulations through graphical representations by using the computer. 


Nn 
a Ls 


T 
2. Notations and Definitions. The following notations have been used if. 
in this paper : 
R Set of real numbers 
Z Set of all integers | ol 
E: Set of all even integers. ca 
O: Set of all odd integers 
N: Set of all positive integers {1,2....}. sh 
N°: Set of all non-negative inegers N U{0} tri 
| ri 
R”: n dimensional space, having co-ordinates indexed 1 through n. | K. 
R"#1: n+1 dimensional space, with coordinates indexed 0 through n. ste 
T: Group of permutation on (1,2,...,n) and 7,41 group of permutation on Ja 
(0,1,2,...,n) 
ui: ith unit vector in R^, ieN and u= Div’ is a vector of ones in R?*!. K' 
: : : : y ZEN G 
V: — jf unit vector in R^*1, ieNy and v=}, „p is a vector of ones in R^*t. 
RP: Non negative orthant of pm i.e. [xeR" 1x2 o). 1 (1) 
A: Set of {~1,+1} | (2) 
B:  Setof {0,1}. | . wh 


Now we consider some standard definations and explanations which will 
be used in this paper. 


2.1 Simplex. An (open) j-dimensional simplex o = o, ae Oe ) is the | 


relative interior of the convex hull of ¡+1 affinely independent points y”, y' Hay | 


called its vertices. ASiPaplexu Soa Ue UF ridere epengot e vertices sits 


Jr 


the o. Again the closure of a (7-1)-dimensional face of the j-simplex o is called a 
facet of c. Two j-simplices are said to be adjacent if they share a common facet. 


The diameter of c is sup Tl - all: x,ze0) = max ll - Weis i} 


2.2. Triangulation. A triangulation G of an m-dimensional subset F of 
R" is a collection of m-simplices that satisfies the following two conditions : 


1. The simplices in G, together with all their faces form a partition of F; and 
9. Each point of F has a neighbourhood meeting only a finite number of 
simplices. 


We denote G/ the set of j-dimensional faces of simplices in G for OS jsm. 
The members of G° is called the vertices of G-and the closures of members of G- 
are called facets of G. 

2.3. Pivot Rule. For a given simplex geG and a vertex y of c, the rules for 
obtaining the simplex of G whose vertices include all vertices of y except y are 
called the pivot rules of G. 

2.4. Mesh. The mesh of a triangulation G is sup(diameter of o:0€G ). We 
shall use the euclidean norm throughout this paper. 

2.5 Some Standard Regular Traingulations. We define some regular 
triangulations of R”. We take the basic grid size of all triangulations to be unity, 
KH and J; are triangulations of R^. K (some times called Kj) is based on the 
standard subdivision of the cube. His very closely related to K and J , is the "Union 
Jack" triangulation of R^. R 

All these triangulations use integer vectors as vertices, let 
K'=- H°- J; = Z" . Also denote the set of central vertices of J, by JW ={ yeJ? : y, 


is odd for each jen ). 
Now we give definitions of some standard trianglations given as 'elow : 


(1) K triangulation is the set of o = (Pras) where y! = y;y =y +u, 1<i<n 


(2) H triangulation is the set of o= (9,....y") where y’ = yy =y eq! 1<i<n 


. Where @ is the ]-th column of the n xn matrix 


= —— 


| The simplex o is the set all xeR' satisfying : nm 
| a | a 
x(1) a(n) | (6) 
1> Y tam Ian) ? - -Dhe Yam)? 0. | 
f = 
(3) J, trinagulation is the set of o= o", song D where ASS Ern a(t) | LO 
1<i<n and y”, 1 isa permutation of N and se" is a sign vector, where s, =+, | pel 
The simplex c is the-set of all seg" satisfying : 
1> Sx) (ene Ya) > > San) ega -Y, m 20. 
The triangulations by J, and by K of any unit cube with integral vertices 
are isomorphic upto orientation of the edges of the cube. 
K,(m) and J,(m) are the triangulation. of S” having mesh ¿K,(m)=m" | Th 
and mesh.J,(m)- 2m ^! . Each Ky(m) and Jy(m) can be expressed as follows : do 
3 (7) 
(4) K,(m)= lo:c- k, (9,7), ifoc s”) is the collection of all such o where 
for 
=(9%,...,9") satisfying y = y" +m"g"" for each ¡eN and where q! is the jth 
column of the (n+1)xn matrix denoted by Q and expressed as 
ma 
-l. 0 
Jeu mo 
O = 
e Q = | 
0 B +] | 
i and its vertices are given by K; (m) e: { yeS” : my, eZ) is: cado E 


(5) J,(m)= (9:0= Ja (97,5), foc Ss" is the collection of all such o whe 
SIKA n : ; i i- 3 í - 
=(y EoI ) satisfying y! = y^! +m epu for each ¡eN and q'? is as define! The 


in (4) and its vertices are given by Jj(m)- K? (m) and Ponte vertice Y 2 
0c <a 0 z ; i 
J, (m) - (yeds (m): my, is even for 1<i<n and my, is odd} for each ieNy | 29 


We now describe two triangulation with continuous refinement of grid a (8) 
The first Ja, basgdconwlkumascitttmodicedbinyangaa ayo oder is due to Baves ET 


e MM 


| Saigal [2]. 


(6) K; triangulates (0,1] X R". Let K" denotes set of vertices 1yc(0,1] xR":y, 
\ 
( =2* for O<keZ and y,/y,ez for ¡eN }. A mapping t: K" — A"! is defined by 
9 | [eO or if y;/yg is odd and 1 otherwise. (Thus / o0) always equals 0). Let z be a 
m | permutation of Ng and K be the set of o= k(y, T) zy EXT y") where 
agis | = : 
| y =y 
| y =y yo ,0si« jan (0) 
| y -y CS aon. ue D ua 
ices | $ an 
| y =y + 2y,0" = jshsn. 
ai | The simplex c is the set of all x e R"! satisfying 
»2 T Apo) 7 Yro) + £540) (o = Jie F ES 1) Es =a) 2 Yo ~ Ko: 
m. (7) J; triangulates.(0,1] X R^. Let J? denotes set of vertices ( yeR"^! jy tes ZN 
for 0< KeZ and y,/yoez for 1<i<n}. 
jth | 
Let Jet = {ved :Y,/yy20 for1<i< n] be the set of central vertices. A 
mapping @:Jy + A" is defined by w,(y)--1 or +1 according as y,/y, is 1 or 3 
mod 4. Let x is a permutation of Nọ and sgA". Then J, is the set of simplex 
| o= ji su ¡ESA where 
| y'=y, 
| y = "HOS UNDE <i< j= 1 (0), 
| yi = yn = » Y o, (x), 
ser? fe cfe q x(k) 
| MD COR UR: SISSE 
ine 


| The simplex o is the set of all 4 e R"*! satisfying : 

jg yx» ` 

| | E (0)(x, (D) $25] 2.25409 Ios Layan 2 Xo Ya 
net) (ua > San es 2 © gm) (Fam Han) 0-20: 


ie | 
and S K; triangulates (0,1]X £^. The set of vertices of K;, denoted by K;” is 
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EL ————————<&x«< — 


E yt = y" 228 0 


ut, (y) 8 L- y, / y, mod 2. | 
Let n is a permutation of Ng with n(j)=0 and i, (y) 2 O if x^ (sj then 


K; is the set of o = k(y,n)- Oey ad") where 


yay yw sic j 
j-l 3 n 5 
A + yp? Y yu 
izü k=j+l 
y=y 429 0 j « kb & n. 
Let z=x-y, n2 (y) and w-z-zgp. Then simplex o is the set of all xeR"i! 
satisfying : 


yo 2 US 2. 2 Ut, 2 Wo: 


nfn) — 

| 
(9) Let A be a non-singular n Xn matrix and G be any triangulation of R” (here 
we consider G=K or J,). Let AG= (Ao:oeG} with Ao={Ax:xeo}. Then trinagulation| 
AK of R” is the collection of simplices o=(y®,...,y") such that all components o| 
Aly? are integers and y'=y'!+a()) where a") is the jth column of A. When A ha 
T's on the diagonal, -1's on the upper diagonal, and zeroes elsewhere, we have the H 
triangulation. 
(10) J'isa triangulation of R” defined as follows : The set of vertices of J isthe 
set of vector veR” with each component an integer, such that there is not precise!) 


one even component or precisely one odd component. Each simplex o = / (yl 


ES (e US = | 
= (y 9 Y ) , Where y is the vector each of the whose components is an even integet 


m=(n(1),..."(n)) is a permutation of (1,2,....2) and s is a sign vector (each $; zi 


| 
| 


With Sy) Sm 71, | 


For convenience s, e" is denoted by & , where e^) is the ith unit vect?" | 
R^. For n24 the vertices of Į are defined by 

yay 

y! Ex y +28, 

yay BP 

yl = yt +6',8 SjsSn-2, 
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by 


then | 
\ 


Ret) 


(here | 
‘ation | 
its of] 
A has | 


the H! 


is the| 
cisely] 


n 


For n=3, each simplices of J' is of the form 6 = j'(y,x 5) 2 (y^, 77,9") where 


V5 Y, 
a(l) 
y' — y -+ 2e ; k 
Des ERE UI 5 212) r(3J 
C Spo are 
y= «A Fe 4) 


Now below we give three triangulations of K. pJ. and J., of R" » [0,1] given 
by Todd [8]. 


Let R" x[0,1]=[xeR"" MERO Removal of the bar from a vector in 


R" x[0,1] denotes its projection into R”. Let g! " be the unit vector in R”*+!: 
thus u!,...,u" are the unit vectors in R^. 
(11) Suppose yeZ” x {0} and 1 is a permutation of (1,2,..,n-- 1). Then k (y,n) 
denotes the (closed) simplex greca where 
y zy, 
y = te qe. 
The triangulation K', is the set of all such k, ( y,n)'s. 


(12) Suppose yeZ" x(1) has all y''s odd. Let x be a permutation of (1,2,...,n.-- 1) 


and let seg" x{-1} be a sign vector, each 5, is +1. Then j, (7,7,5) denotes the 


simplex pinoa] , where 
y? ee 7 
Quy s SS 
The un J! is the set of all such 7, (y,1,5)'s 
13) y i} is the triangulation derived from J',. Let yeZ" x{1} has all y''s odd. 
Let x be a permutation of (1,2,...,n-- 1) with m(j)=n+1. Let ser" x{-1} be a sign 


vector such that, for j<k<n+15n(k) is -1 or +1 according as ¥*" is 1 or 3 mod 


4. Priced 
Then J, (9,1, 5) denotes the simplex [x^ gn |, where 
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y' 2 y" «2sn(k)u j «E n1. | 
um | 


The triangulation J . is the set of all such J,(¥,7,5)'s. | 


3. Measure of Efficiency for Triangulation. Todd [8] has indicated tha | 
the theoretical measure of the efficiency of different types of es ior used 
to be compared by counting the number of simplices into which unit cube is divided 
However, Todd [8] has also observed that many triangulations have differeni 
efficiency but yield the same number of simplices e.g. K, H or J 1 triangulations) 


] vs G MS 1 ‘al | 
having n! simplices and K' or Jy in [xeR BST S x, < 1; having 


(2^ *1—1) n! simplices. | 
Saigal Solow and Wolsey introduced the concept of diameter, for, those | 
triangulations that subdivided unit cube. This is the maximum over all peius of 
facets of the triangulation lying in the facets of the cube as well as of the maximum 
number of simplices that form a path of simplices linking the two facets. They 
computed the diameter of K triangulation and an obvious extension to K; 2 
found them compareble, Saigal [5] calculated the diameters of triangulations A | 
and H and obtained results and suggesting that the number of iterations using į 
increases with n? while that using H increases with n°. It is easy to see that Ur 
diameters of J, and K are equal. | 

The diameter of triangulation is a "worst best case" measure similar to tlt : 
diameter of a polytope which is of interest in linear programming. 

4. Main Result. First we shall give the concept of average direction! 
density, then formula for obtaining different triangulations and numerical result 
for obtaining a.d.d. for different values of n. At last we shall give the computers” 
graph of average directional density of different triangulations. | 

4.1 Average Directional Density. Todd [8] proposed the concept % 
directional density as an alternative measure, for the 


nl 
comparison of efficiecy | 
the triangulations and is global in nature. 


n 
For x,deR"and 150, let [x,x+Ad] denote {x+pd:pe[0,A]}. Let G 2 


| 
T 
triangulation Sham cimesh Jy), Let N(G,x,d,A) denote the numb? | 

5 li of G free! urukul Kangri Collection, Haridwar. An eGangotri Initiative 
simplices ol 


the 
wrung [53 Ad] divided by 2. Let N(G,x,d,A) denote 


limit as p > 2, (if it exists) of the average of N{G.x,d,7.) for x uniformly distributed 


in B(0,p) = [xe |x| < p}. Let N(G,d) be the limit as 4 — (if it exists) of 
N(G,d,X). Finally, N(G) is defined to the average of N(G,d) for d uniformly 


distributed on dB" = {deR" Jali = 1). 


We call N(G,d) the directional density of G in direction d and NIG) the P 


| average directional density. It should be noted that N(G,x,d,1) is the number of 
that | 


| simplices met per unit step size. To eliminate the effect of starting point, we average 
Used | over x lying in a large ball. To eliminate any effect on the ending point, we let || 
ided, | : E Bins 
kb. Finally to get a measure independent of the direction d, we average over i 
reni if 
tions) AB" - | 
aa For any triangulation G of R”, G"-2 is a countable collection of sets of & 
ving | 1 
dimension n-2. Thus for almost all x, [x 4 Ad :AeR] meets no simplices of G^ of $ 


hose) dimension less than n-1. If x and x+Ad lie in n-simplices of G.[x,x+2d] meet one 
rsof| more simplex of G than (n-1)-simplex of G. Thus the computation of N (Gx,2.d) 
nun! lan be made by counting the number of points of [x,x+d] lying in facet of G, for | 
hey! almost x. | 
and! From these observations the following results are obtained by Todd [8] and | 
ag K| Vender Laan and Talman [19] : 


ng f (i) N(&,d)- Y |d,|+>., Ia, -a | | 


the 
| © Maa) Ya], a aaa; | 
, the | Í 
j (iii) NE, d) X Ehl 1 
ona| à i 
ae) 0 NU-S Hla raf) | 
sed 3 i | 
| M Let g, = 2F(n/2)/(n—1) dnt ((n-1)/2), then " 
t 0 i | 
| if 
y’ (a) N(K,)=N(J,)=(n+v/2(n/2)8,) | 


aN (b) N(H,)= Y, (n+1-d (Dg, i 


O N(K )=1/2 m+} [n is, | 
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the (d) NI) (n VE 


R 


O 


IMD, (n+ 1k) (4) + eS 


(e) N(J,)=(1/4 
(1/2) SY Ds j+l- (i+ j2)/ 02)" |e, 


g, {n(n +1)/8)' Ee +2)) ifn is even 
g, {nr (n+1 yer | 
Table 1 gives N(G)/g, for various values of a n and G equal to K, Jp H "n 


N(A* K)= 
(f) ( n4 1) if n is odd 


1 

A*K and Y. 2 
TABLE 1 3 

The a.d.d. of the K or J}, K , H, A*K, J, and J’ triangulations for various 4 
values of n(mesh equal to Jn ) à 


TRIANGULATIONS G 


y 
Kor J, K H A*K J al 
1 0.7 1 1 0 3 
3.41421 2.44949 3.4142 | 2.44949 1414. | | 
1.24264 5.4641 7.5604 | 4.89898 4.2426 JE 
12.4853 9.94936 | 13.7047 | 7.74597 8.4852 | 
19.1421 16.0797 | 22.089 | 11.619 143421. | rh 
59.9117 60.0249 | 82.0105 | 33.541 50.9116 || for 
163.492 197.748 | 269.433 | 87.6356 148.4924 | bec 
288.701 392.453 | 534.189 | 151.987 268.7006 | inc 
645.183 1044.16 1419.95 | 334.066 615.1829 | ver 
1782.41 3636.58 4942.12 | 910.357 1732.412| slo 
7100.36 20108.5 | 27316.5 | 3588.52 7000.351} J, 
4.2. Asymptotic Behaviour of A.D.D. of Different Triangulations !! 
^ G is any triangulation then Asymptotic behaviour of N (G)/gn for G=K Jul | a 
x 


J,,H and AK are given as follows : 


i ami 
(i) For K or J, triangulation : N(K Ve, =~ n2/V2 = N(J J= N(J"), 

(ii) For K* triangulation : N(K *)/ g, =~ nn?" /16, | diff 
(iii) For H triangulation: N (H)/g, =~ 4n??15 dey 
i ; n , b t 
(iv) For A*H triangulation : N ( As K) le, =~ n? [Af8 ) m 

Ui > | 
v) For J; triangulation : N(J s, yg, == (n nl), | fast 


n 
In table SS worn baret Ayi oki ¢ behaviour of different tr rango ju 
— o o S SSSS o TP EN 


| 


for different values of n. 
TABLE 2 


The asymptotic behaviour of a.d.d. of the K or J, K ,H, A*K and J’ 


| 
| triangulations for various values of p (mesh equal to fn) 
| d 


Tringulations G 


S [1 : 0.70711 0.19643 0.26667 0.35355 1 

| {2 2.82843 1.11117 1.50849 1.41421 5.65685 | 
3 | 6.36396 3.06202 4.15692 | 3.18198 15.5885 l 

1s 4 113137 6.28571 8.53333 5.65685 32 
5 17.6777 10.9807 14.9071 8.83883 55.9017 | 
9 57.2756 47.7321 64.8 28.6378 243 | 
cee 15 159.099 171.172 232.379 79.5495 871.421 | 
20 282.843 351.382 477.028 141.421 1788.85 | 
30 636.396 968.295 1314.53 318.198 4929.5 | 
| 50 1767.77 3472.4 4714.05 883.883 17677.7 | 
| 7071.07 19642.9 26666.7 3535.53 100000 
| 6. Conclusions. In table 1 it has been observed that acL.d of H and J l | 
1 trinagulations increases very fast as n increases and both remains nearly same | 
6 | for 2730. However, their values remain parallel for 30<n<50. But values for J | | 
24) become faster than H for 50<n<100. For J' triangulation, in beginning its value | 
06 | increases very slow rate but as n increases its value pick up very fast and become | 
29 | Very close to a.d.d. of K or J,. Rate of increase of a.d.d. of A*K triangulation is the Ñ 
12 || slowest of all the triangulations. In the same way rate for K is faster than K or | 


a J, but less than H and J T 

| From these observations one can conclude that those triangulations whose 
K | a.d.d. increases at faster rate are considered to be inferior than the other. On this | 
| basis H can be considered inferior among all the triangulations and A*K is superior | 
among all the triangulations. 

Again if we take a graphical representation of tabulated values of a.d.d. for 
different triangulations, one can easily observed that for n<15 there is not much 
deviation in the different graphs of a.d.d. for most of the triangulations. However, 

etween n=15 and 20 the deviation between the value of H and K became faster, 
but that of H and J , remain the same, that is H and J , begins to deviate at the 
faster rate than K' . Deviation of A*K is not much from its previous values before 
jo, "15. In addition to this enaicanaalscaeaily, abasuved tbatadt (eap) and J remain 


Se 
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nearly same for values of n<20 and their deviation is also less than H, J ; and Jn 
K , but greater than A*K. However, from n=30 and onwards the values of a.d.q, 
for triangulations H and J , as well as K. steeply deviated from all othe, ' 
triangulations. Whereas the deviation between A*K and that of K(or J,) and J jg 
not much upto n=30 but deviation between them pick-up between 30 and 50 thereby 
the deviation of K (or J,) and J' became steeper than A*K. Though A*K is superior 
of all the triangulations. But in practice K (for J j is found to be much convenient 
to determine the fixed point of given mapping. | 
Table-2 compares the asymptotic a.d.d. for different triangulations and | 
chart-II gives their graphical representation of asymptotic a.d.d. based on tabulated | 
value of table-2. Both tabulated values as well as graphical representation again 
highlight, A*K as superior among all the triangulations. | frac 
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ABSTRACT 

The purpose of the present paper is to derive a number of key formulas for 
fractional derivatives of generalized multiple hypergeometric functions of several 
variables, multivariable H-function and genralized multivariable polynomials. 
Each of these formulas can be shown to yield interesting new results for various 
classes of generalized hypergeometric functions of several variables. Some of the 
applications of the new formulas provide potentially useful generalizations of 
known results in the theory of fractional calculus. Our results include all the 
results of Chandel-Kumar [12] as special cases and all recent results of Ram- 
Chandk [32] as special cases for the Fox-Wright generalized hypergeometric 
function of Wright [54]. 
2010 Mathematics Subject Classification : Primary 33C99; Secondary 26A33. 
Keywords: Generalized Fractional derivatives, generalized multiple 
hypergeometric functions, Multivariable H-function Generalized multivariable 
polynomials, Fox-Wright generalized hypergeometric function, Bessel-Maitland 
function, Wright generalized Bessel function, Mittag-Leffler function. 

1. Introduction. The theory and applications of fractional calculus are 


based largely upon the familiar differential operator ,D7 defined by (cf., e.g., (31, 
P.49); see also [47, p.356]) | 


: y Fat (Re(«) « 0) 


194 T(-a) : 
1D pDE{f(x)}= a 2 0< Rela) m.m e N,) 
dx" p ET if (x) j = : 


Where N, = NU {0}, N ={1,2,3,...}- 


For B=0, equation (1.1) defines the classical Riemann Liouville fractional 
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derivative (or integral) of order a(or-a). On the other Tand, when B— œ, the 


LIU i 


* 


equation (1.1) may be identified with the definition of the familiar Wey] fractional 


derivative (or integral) of order a(or — 0) (see for details, (16, Chap. 13] and 134p, 
For the sake of simplicity, the special case of the fractional calculus operator ,D' 
when B «0 is written as D; . Thus we have 


(12) D'=,D*(aeC). 


For 0<a<1;8,n,~eR;meN, the generalized modified fractional derivative 


operator due to Saigo is defined in Samko, kilbas and Marichev [34] as 


sa d( ¿ón 
(1.3) D (x)= =| ( 


| da T (i-a) x" 0") F(8-0)1-nd-o;1-t" /z" | (gap. 


The Multiplicity of y” _ 4” in above equation is removed by requiring log( ag” =t") 
as real for y" —¿" > 0 and is assumed to be well defined in the unit disk. 

It is remarkable that | 
(4). Dir f (x)= Def (x), 


where D“ is the familiar Riemann-Liouville fractional derivative operator defined 
by Miller and Ross [28] 


For 0s a<1meN;8,n,xeR; > max(0,8— m), the refined form due to 
Bhatt and Raina [2] is given by 


apr —ljm D IP — 
(1.5) Dee E) 1) } = oe 


Making an appeal to the Saigo modified fractional derivative operator Dj 
Miller and Rasso [28] investigated fractional derivative formulas. Kilbas [21] 
established the fractional integral formulas for the Wright ([54]; see also Erdely 
[14]) introduced the generalized hypergeometric function s Fox-Wrigh! 
generalized hypergeometric function, defined by 


(a,, A;) 


p 
Lp 2 TIT (s; + AK) z! 


izl : 


(1.6) pY (2) nu AN y E e. 
: jaj . hz0 a 
Cc . Burukuik ngii Collect ikrar eho Initiative 
jal 


suit 


and 


defi: 


(1.7. 


whe 


ed 


1 q p 
where a;b; € C; A, > 0, B, > 0L Y -V A, 20; (A.B, #0)(i= jm Dana 1,...g), for 


J=1 i=l 


suitable bounded values of |z|. 


For details of conditions of its existence and the H-function due to Mathai 
and Saxena [26], see Kilbas [21]. 


Wright [55] also introduced the special case of (1.6) (called Wright function) 
defined in the form: 


RE A 2 
(1.7) o (o. ; z) = QW mi = Loro) e > 


where œ,zeC and Be R. 


Kiryakova [23] introduced a function J5 (z) called Bessel-Maitland function 


or Wright generalized Bessel function defined by 


: MM th 
8) Ji (z)= 1,52) =N a eR 
o Vo) reve] kl 
Other special case of (1.4) but generalizing the classical Mittag-Leffler 
function (Erdélyi [16]) is given by Kilbas et al. [22]. 
Srivastava and Garg [38] introduced a general class of multivariable 
Polynomials defined by 
Ind teluk sb 


(1.9) Si (at ay X,) = y Por O a 


REPE 


Where h,,...,h, are arbitrary positive integers and the coefficients A(L:h, yk), 


0» 


ON. MA. 5-1, the polynomials (1.9) reduce to the polynomials of 
Stivsataya [37] defined by 


(L;k, EN i= Des) are arbitrary constants real or complex, N, = NU(0J. 


Wal 3 
(10) S^). X Us eer ES D 
k=0 5 


Where h is arbitrary positive integer and the coefficients A,, (1; e N,) are arbitrary 


Con 
Stants, real or complex. ku Kangri 


j gwar. ie Ab i ls) of 
The computation of fractional derivatives CO Aes iiid : 


special functions of one and more variables is important from the point of view of 
the usefulness of these results in (for example) the evaluation of series and integrals 


(ef., e.g., [29] and [52]), the derivative of generating functions [46, Chap. 5] and | 
solutions of differential and integral equation (cf. [29] and [39], Chap. 3. see also f 


(27, 30 and 51]. Making an appeal to the operator (1.4), Chandel and Vishwakarma 
((61,17)) have obtained fractional derivatives of confluent forms due to Chandel. 


Vishwakarma [5]) of Karlsson's multiple hypergeometric function ear 


and 


other multiple hypergeometric functions of Lauricella [24], Chandel [3], Chandg 
and Gupta [4] including their confluent forms. Srivastava and Goyal [43] derived 
fractional derivatives of the multivariable H-function of Srivastava and Panda 
(1481-1501). Srivastava, Chandel and Vishwakarma [40] obtained several new 
fractional derivative formulas involving the multivariable H-function defined by 
Srivastava and Panda (see (50, p.271, eq. [4.1] et Seq.]) and studied systematically 
by them (see [48] [50] also [44]). 


Further for special interest, Chandel and Vishwakarma [8] and Chandel- ` 


Sharma [11] derived fractional derivatives involving hypergeometric functions of 
four variables defined by Exton [18] and Sharma-Parihar [35], while Chandel- 
Sharma [11] established fractional derivative formulas for their own hypergeometric 
functions of four variables ( [9],[10]). Recently, Chandel and Kumar [12] derived 
generalizations and unifications of various key formulas of Srivatava Chandel 
and Vishwakarma [40]. Very recently employing the operator (1.3), Ram and 
Chandak [32] derived a generalized derivative formula involving the product of 


Fox-Wright generalized hypergeometric function » V, defined by (1.6) anda general 


class of multivariable polynomials defined by (1.9), Some special cases are also | 


discussed. 


> ; f 
In the present paper, employing the operator (1.3), with the motivation of 


Chandel-Kumar [12] and Ram-Chandak [32], we derive generalizations and 
unifications of various key formulas of Chandel-Kumar [12] and Ram-Chandak 
[32]. Each of these formulas can be shown to yield interesting new results [or 
various classes of generalized hypergeometric functions of several variables an! 
genearlized hypergeometric functions of one variable. 
key formulas provide potentially useful genearlizatio 
theory of fractional calculus. Some special cases are 

2. Key Formulas, In this section, making an 
derive the following key formulas on generalized fr 
multiple hypergeometric function of Srivastava 


Some of applications of the 
ns of known results in th 
also discussed. 


actional derivatives involv 
-Daoust ([41],[42]; also $ 
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appeal to the result (1.5), V^ ¢ 
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Srivsatava-Panda ([48]-[50]; see also Srivastava, Gupta-Goyal [44]) 


: ; "EET ap ali ; 
multivariable polynomials due to Srivastava-Garg [38]: pu 
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Where ps; 


cp; p» 1$ generalized multiple hypergeometric function of Srivastava and 


Daoust ([41],[42)), 0< a, < Lm, e N;B,,n,x; e R;p, > max(0,B, -n,), 
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function due to Srivastava and Panda ([48]-[50]). 
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F LE an is generalized multiple hypergeometric function of Srivastava and Daoust 


([41],[42]] while S^*-^ ' (x,,...x,) are generalized multivariable polynomials due to 


Srivastava and Garg [88], defined by (1.9). 
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3 multivariable polynomials due to Srivastava and Garg [38] defined by (1.9). 0% | 
result (2.5) among others also includes (2.4) as special case.’ ji 
3. Proofs of the Key Formulas. In this Section, we prove the results | 
Section 2. 
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Therefore, left hand side of (2.2) 
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Proof. of (2.3). For brevity, if we denote 
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Which can be adjusted in the form of right hand side of (2.3). 
Proof of (2.4). For brevity considering 
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which can be reformed as right hand side of (2.4). 
Proof of (2.5). We can write 
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[By making an application of (1.5)], 
which can be expressed in the form of right hand side of (2.5). 
3. Special Cases. In this section, we mention the special cases of our 
results. 
Special Cases of (2.1) 


(a) Form,-1,- Os, = L...,r; (2.1) reduces to Chandel and Kumar [12,p.107, (2.1)]. 


D) For r- Lin, - LB, a, (2.1) gives Srivastava, Chandel and Vishwakarma 
[40, p. 567 (3.1)]. 

(c) For p= 2,m, =m, - Lp, =0,,B,=0,, (2.1) reduces to Srivastava, Chandel and 
Vishwakarma (40, p. 568 (3.3)]. 

Special Cases of (2.2) 

(a) For m,=1,B;=a,;i=1,...,r; (2.2) reduces to Chandel and Kumar (12, p.108 
(2.3)], which also generalizes the earlier results due to Srivastava, Chandel 
and Vishwakarma [40, p. 563 (2.1), p. 564 (2.3)]. 

Special Cases of (2.3). 

(a) For m,=1,8,=a,,i=1,...,r; (2.3) reduces to the result due to Chandel and 
Kumar [12, p.107, (2.2)]. 

b) For, 1,m, =1,B, = o, (2.3) reduces to the improved version of the result due 
to Srivastava, Chandel and Vishwakrama [40, p. 567 (3.2)] 


() Due to general nature of the functions involved in (2.3), specializing the number 
of variables and other parameters in (2.3), we can get several interesting known 


E and unknown results. 
“pecial Cases of (2.4). 
() For y - LA, =0,n LX 0oSuulg! chee ii Amebadethemesult due to Ram 


o 
144 So 


and Chandak [32, p.53 (10)] involving Fox-Wright generalized hypergeometri, 


h Yel | 
" of Srivastaya | 


function ,V, and generalized multi-variable polynomials Sj" 
and Garg [38]. 
(b For r=1,1,=0,n=1,0' =0 and f, =o, (2.4) includes an interesting result 


due to Ram-Chandak [32, p. 53 (11)] for the Riemann-Liouville derivative prc 
operator-defined by Miller and Ross [28]. 


(c) For s=1, the polynomials S^ (x,,...,x,) reduce to the Srivastava polynomials 14 


S! defined by (1.10), therefore, for r= LA, =0,n=1,0', =0,8, =a, and s, -1, 


(2.4) gives an interesing result due to Ram-Chandak (32, p. 54(12)] involving 
Srivastava polynomials. 


(d) Similarly all the following results due to Ram and Chandak [32, p.54 (18), 
(14), (15)] are included in our result (2.4). Since for /20,A,, = LS OS 
therefore our result (2.4), also gives the result due to Kilbas [21] as special vali 
case of Ram and Chandak [32 (15)] - whe 
(e) Since Wright function pla, B; z) defined by (1.7) and Wright generalized Bessel ae 


function J°(z) defined by (1.8) are special cases of Fox-Wright generalized 


hypergeometric function ,,w,(z) defined by (1.6), therefore, all the results due 
to Ram and Chandak (32, pp. 55-57 (16), (17), (18),(19)(20),(21)] are included in 
our result (2.4) as special cases. FER 


Special Cases of (2.5). Our result (2.5) includes (2.4) along with all results al (4.4 
Ram and Chandak [32] as special cases. j 


4. Other Special Cases as Application of one Fractional Derivative | 
In this section, making an appeal to (1.5) and one fractional derivative we der!" 


| 0... 0 9 Negi]. | fp) . ale) Whey 
P OE E UR E iu yx. oo IE ds 


Wa v^ |: ta») Volks ;[ (a): | i 


— 
Q 
— 


9" | fu 135A, | [pa +n=B: Nis ag heal : 
c) 2 Wi Dt -B 5 he lft Ena: UV : 


= xt y Arz sup [(a) :9 pr) 


C32:D 5. D K valid 


(4.6) 
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wore] s 
(We (a): ge je YX | 


provided that 0< d <1,meN;B,n,A, >0,x € R,u > max(1,B—n) and 
Gee and $ A p" 

1e, Vj + 2,87 - 307 - $97 20 IL... 
Jel j=l J=1 J=1 

(4.2) Dy, a fx ae Ey (u —B,6, yb, USC ao oe” )} 


alt B= > Pu) r(u: Ln- B) 
= El" Ez 3 
F(u - B)T (u- n- o) (u+n- B,D, 5-50, s+ T= oz yx" joer y, x ) 


valid if Os o «1, me N;B,n,xe Rh» max(0,B — 1), hnal « L.....|y, x] cd, 


where p is well known Lauricella function [24]. 


(4.3) Dern E FY) (u +9-0,b,,..,0, N= Bro xu" sy, x" ) 

Ap) (a+n-B) pin 
l(u -B)T (u-- n- o) 

Provided that all conditions of (4.2) are satisfied. 


EO Dinda Piu en ab, by 1930"... 9,2") 


(mb iC 0 B; yx Ela a): 


aan 


Se ewe (ome B) zc 
F(u-B)r(um-a) E 


Where all conditions of (4.2) hold true. 


(4.5) Dyin ommo) (u -Bbo bp ;u t n- Broyx" s y, x ) 


(u+n=B,b,,.. > lo gn B; yx pa AE 


A $ 
= u+n-B) (n) : _R- a 
=> HE A ub ob HEN Biyax dt X 
(i -B)F(u € n- o) D ( 1 ¡Y p 
Vah . 

alid if al] conditions of (4.2) are satisfied. 

(4.6) Dem E = 


(u-1)m Re | o 
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au 


({u-B-1)m Z3 * 
EEUU B) (unen-Biu-Bien-o;z , 
T(u-B)F(u+n=a) 


provided that 0<a «1,me N,B,n,x e R,u>max(1,P- n). 


Am Denat e (u- Bie n- oci e n-Bia")) 


NE MEO D) (on 
Ee cms) 2) 


where all conditions of (4.6) are satisfied. 


co Dj" mu Baz") 


T(u-B)r(un- o) 


which holds true if all conditions of (4.6) are satisfied. 


Fia n-Brucen- oz"), 


(4.9) pada aa F (u +n- aua” ) 
$ at r(ur(un- g) &( F 8 B ^) : 
A —B;u-B;x 
T(u-B)r(u+n-«) 1^1 u T u , 
where all conditions of (4.6) hold true. ` 
(4.10) Deen Ea ¡A (u E Bu -n- o; x" )} 
AMA Eme B) 
NN en a Na J, 
provided that all conditions of (4.6) hold time. 


Gan Data On urnas en) 


_2 (ur wenp) 
“Ta Drama) Abu Bio"), 


which is true if all conditions of (4.6) are setisfied. 


apo), Ln 
(4.12) DEA a ) 
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val 


whi 


(4.1 


pro 


vali 


(4.1 


whi 


4 ME (wrp +n- B) (1 E -(p*n-f) 
F(u - B) (& 9 n- o) ) 


valid if all conditions of (4.6) are satisfied. 


(4.13) DI F (u-B;u+n-o,u+n-B; 2" )} 


A T ee s 
r(u-B)r(u+n-a) i 


where all conditions of (4.6) are satisfied. 


(4.14) D Jaen ¿F, (u — Bu -T7- Qu x" ) 


S x" (u)r (u+n-p) E E) 
F(u -B)T(u * n- o) Pa 


provided that all conditions of (4.6) hold true. 


(4.15) Dyin {x" (UB e n=") 


x "T(u)r(u* n-B) _ myn) 
r(u-8)r(u+n-a) (1 z) 


, 


valid if all conditions of (4.6) are true. 


(416) Deen Tea F, (u +n-œ;p p+n=P;x" J 


tecum ae (ERE el) s su 
F(u-B)r(u*n-e) ° nebat, 


Which holds true if all conditions of (4.6) are satisfied. 
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ABSTRACT 


The purpose of this paper is prove some strong convergence theorems of 


the modified Ishikawa iterative sequences with errors for unformly equi-continuous 
and asymptotically quasi-nonexpansive mapping in the setup of unifromly convex 
Banach spaces by using condition (A) instead of comletely continuous or 
demicompact condition. Our results improve and generalize the corresponding 
results of Rhoades [6], Schu [7,8]. Tan and Xu [10,11] Xu and Noor [12] and many 
others. 
2000 Mathematics Subject Classification : 47H09, 47H10, 47325. 
Keywords: Asymptotically nonexpansive and Asymptotically quasi nonexpansive 
mappings, Fixed points, The modified Ishikawa iterative sequence with errors, 
Strong convergence, Uniformly convex Banach space, Uniformly equi-continuous 
mapping, Uniformly Holder continuous mapping. 

1. Introduction and Preliminaries. Let E be a real normed linear space, 
K bea nonempty subset of E. Throughout the paper, N denotes the set of positive 
integers and F(T)=(x:Tx=x) the set of fixed points of a mapping T. Let T:KK be 
à given mapping. 
(1) T is said to be asymptotically nonexpansive [2] if there exists a sequence 


{k} 1,00) with limk, =1 such that 
(reg »y|s ^, li — 3, (1.1) 


for all X,ye K and neN. 


D Tis said to be asymptotically quasi-nonexpansive if F(T)#o and there 


exist; ; E 
Sa sequence {hi édp Ebro), wath dim. A Such An a Initiative 


t 


| 
[Tx - p| Sk, |x - p]; (19 | 
for all xe K,pe F(T) and ne N. 


(3) T is said to be uniformly L-Lipschitzian if there exists a positive constant 
L such that 


|T" x- Tys L|x-yI, (13 


for all x, ye K and ne N. 
(4) Tis said to be uniformly Holder continuous [5] if there exist positive 
constants L and « such that 


[Tx - T"y| € L|x- yl, (14) 


for all x, ye K and nec N. 
(5) T is said to be uniformly equi-continuous [5] if, for any €» 0, there exists | 
$»0 such that 


Ira -T" j| <e (13) 
whenever |x - y|« 8 for all x,y € K and n 2 1 or, equivalently, T is uniformly equi 


continuous if and only if |T"x, — Ty, | > 0 whenever | 50 as ne. 


2, =I, 
Remark 1.1 (i) It is easy to see that, if T is asymptotically nonexpansive, then il 
is uniformly L-Lipschitzian. 
(i) If T is uniformly L-Lipschitzian, then it is uniformly Holder continuous with 
constants L > 0 and a = 1. 
(iii) If T is uniformly Holder continuous, then it is uniformly equi-continuoUs f 
This class of asymptotically nonexpansive mappings was introduced b 
Goebel and Kirk [2] in 1972. They proved that, if K is a nonempty boundel 
asymptotically nonexpansive self-mapping of K has a fixed point. Moreover the 
set F (T) of fixed points of T is closed and convex. Since 1972, many authors haw 
studied weak and strong convergence problem of the Mann and Ishikawa iterati"| 
Sequences (with errors) for asymptotically nonexpansive mappings in Hilbe" 


Spaces and Banach spaces (see (2,6,7,8,10,11,12] and references therein). 
Recenyly, Liu [3] studied modifi 


: o 15 
fi : : ed Ishikawa iterative sequences with D 
ion asymptotically quasi-nonexpansive and unifornly Holder continuous mapp 
= mu Ka Banach spaces arid established some strong converg” | 
eorems which extended : | 
Th e purp SO parut Kang Comesponding Alan and Xu [11]. (3) ftf 


15 paper is to extend and improve some resuits of 


for ; 


(19 


tant 


13) 


tive 


i 
| 
| 


uniformly equi-continuous and asymptotically quasi-nonexpansive mappings. Also 
our results improve and generalize the corresponding results of [6,7,8,11,12.] and 


many others. 


In order to prcve the main results in this paper, we need the following 
lemmas : 


Lemma 1.1. (Tan and Xu (101). Let[o, Y" |, (B, J^. and [r, ),., be sequences of 
nonnegtive numbers satisfying the inequality 
Onn € (14 p, Jo, +7,,V, 2 1. 


n? 


If 3B, « and 2,5 <, then Lima, exists. In particular, {œ}; has a 


n=1 n=1 


subsequence which converges to zero, then Lima, =0. 
no 


Lemma 1.2. ([1]) Let X be a uniformly convex Banach space and B (0) be a closed 
ball of X, Then there exists a continuous increasing convex function 


8:[0,00) > (0,ee) with g(0)=0 such that 
[Ax + uy + vell Sel ul vla Aue (lx ul) 


for all x,y € B, (0) and 2,1, y e[0,1] with X*uy-1. 

2. Main Results. Now, we give the main results of this paper. 
Lemma 2.1. Let E be a normed linear space and K be a nonempty convex subset of 
E. Let T:K>K be an asymptotically quasi-nonexpansive mapping with F(T)s6. 


and a segeunce (i, Jc fo) with Limh, =1 and Y (12 -1)<=. Let tr} be a 
Sequence in K defined by 


2 smpn y 
Yn = qx, te bT x, ar C,U,; 


=a,x, * b, T" y, +c,U,, NZI, (2.1) 


nl non 


Where {un}, (v,) are bounded sequences in E and (a, ), 05, ), GG, 3,027 3, (5 ttc, ) 


x 


are sequences in [0,1] and a, +b, c, ai +b, +c, =1 with the restrictions 


2 e. <eo and $^ b.c, «c». Then we have the following: 
E Lima, — p|| exists for any p eF (T), 


(ii) Limd(z,, F(T)) udis where d(x,F(T)) denotes the distance from x to the 


Hh=>0 
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Proof of (i). Let pe F(T). Since fu, and {v,,} are bounded sequences in K, Sp v 


| 


can set TI 
M = max (sup, - p|,sup|u, = al. 
n2l nz 
Then it follows from (2.1) that s 
ly, - pl “asa, € T", +c,- p| 
(T's, =p) +6; (o, =p) s 
=la’ (x, - p)t 6; (T"x, - p)+c; (v, - p 
bn : 
«a; lx, - p|+b,|7"x, - pes lv, - pl sul 
€ a; lx, - p|-- bk, læ, — p] * c; llo, = pl no) 
S[a, + b; ], t s, - p] ci len — pl | an 
AASA E 
Sele, - p] c; M. eaj 9 
Again from (2.1) and (2.2), we have : (ii) 
(Eta zs pl = lanx, + b, T^y, » CU, m p| Thi 
= la, (x, = p) 5, (T^y, - p)+ c, (u, — p)| Pr 
s a, tn pl + b, IT", z pl + C, lu, = p| 
< a, en - p|*5,5, ly, — p|*c, lu, - pl | E 
<a, |x, - p+ by, [k,en = pl tes Mo, lu, — p| f an 
1 <a, +b, ]%? |x, - p|- 5, RM c, M - 
| . AS 
i 5 [i =C, |k; Ix, -p| ap b kc M ar c, M 
| r= 
< kilen = p|+b,k,c M + c,M 
< 2 (23 lu, 
X [2 x (E s 1)]ls, PN pl + (b.c; b, ar €, )M 


: UEM 4 re 
since = - o Y 
Lama js o s» and 2. <œ, it follows from Lemma 14" 


know that lim|s, - pl 


| exists. This completes the proof of part (i). 
| Proof of (ii) From conclusiogietiparb. (dridwer.kavEangori Initiative 


) We 


2.2) 


lx... — Pils [1 + (R2 — 1) ls, - pl (6,cik, +c,)M. 


This gives that 


d(x, ,,F(T)) < [1+ (k = 1) |d(x,,F(T))+ (b,c;k, +c,)M. 
Since Ns UR -1) « ae «co and È bic, «ce, it follows from Lemma 1.1, 


we know that lim d(x,,F(T)) exists. This completes the proof of part (ii). 


Theorem 2.1. Let £ be a uniformly convex Banach space, K be a nonempty convex 
subset of E and T:K>K be a uniformly equi-continuous and asymptotically quasi- 


nonexpansive mapping with F(T)«$ and a sequence {k,}c[1,<) with limk, =1 


and Y" (k? -1)<oo. Let (x, )bea sequence in K defined by (2.1) with the following 
nzlN 1 n 


restrictions: 
(i) 0<b, <b<l and b, 


nl 


€ b, for all n>1. 


(ii) 15,79, (i) ime, =0, (iy) Y" c 2e» and Y^ bc «e. 


nz] 7 
Then liminfl|x, — Tx, | 2 0. 
Proof. Since T:K—K is asymptotically guasi-nonexpansive, we have 


Ir" Y, = p| Sk, = pl 


»,-p|s E 


forany pe F(T). By Lemma 2.16), we know that lim||x, — p|| exists. Hence (x, - p} 
and [T^y — p} are bounded sequences in E. Set 7 -sup(lx, - pl: 21}, 
T, 5 sup (zu, — p|:n> 1), 1, = supl[]u, - p|:n2 1), 7, =supllo, -p]:n2 1) and 
"=max(r, :i=1,2,3,4) for any fixed pe F(T). Then we have fx, - PT", - p]; 
lu, = D).[v, — ple B,(0) for all n2 1. By using Lemma 1.2 and (2.1), we have 


a’ (x, - p)+ (T^, - p); (v, -»f 


ly, ES pl S | 


Tx, = pl +c/ llo, - pl - a;b;&(|s, -T"x, |) 


x, — pi xb, 


Sa, 
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Xa dux — pl + x, Pl tC 


non 


<a; «Ies - f ec 
= [1-01 ]A? px, - pf etr? 


<k? |x, = pl rcr. (24) 
Again, using Lemma 1.2, (2.1) and (2.4) we obtain 


len- pl = le, (x, —p)+b,(T"x, — p)+ c, (u, — p) 
<a, lx, — pi tb, Ir», - pl +C, loe, — pl - ab, (lx, SD, |) 
Sa, |x, -p| t bk |y,- pl tc, Ju, -p[ -a,b, gll, -T"y, I) 


E a, |x, = pl t bk; Es ll, - à pl iis er + e,r? T a,b,g (læ n T |) 


) 


s [a, + b, ]h; lx, ex pl + b,c, kir” ar um 3 a,b, 8 (|x MASS T"y,| 


=[1=c,]k*llx, - pl) + ber *e,r^ ^apb,g(|x, - T^y, |) 
Shi |x, - pi + * (b.c; ki +c, y - a,0,8 (ll: y T^y, |) 
sn«(u -1]ls, - 


Note that Ya =1)< 0 is equivalent a 1)< 0 and so, setting 


" «(b.c ki +e, r= a,b, (lx 27 T^y, I) (i) 


=r* (ki -1), then DEED Furthermore, since g:[0,0)—[0,-) is ? 


continuous increasi i : i 
creasing convex function and ls, p yn} is a bounded sequence! , 


E, we assert that 8 (lx, -T" 


Yn ) is bounded. Set 6, = ¢,8([x, -T"y, I), we have | 
| D 9, <>. Since (%,) is bounded, and by sE yb eo 80 


| Pra TA & oo . Now, Set 


8, =P, +0, +b ek? jue 


nnn 


Then io. «oo, 
(2.5) that 


ka -p[ sk, -p| - 


By the assumption (i), we have (1-b,)>(1-b). It follows f°" | 


CC- 0. pa ein b Si) ak ee aby Dial tc - 


exis 


1 +3, (2.6) 


which leads to 


(1 x: b)b,g (lz, Sa T'Y, I) < Ix, i pl Sa F pl + ô, , (2.7) 
and 
(1- )b, se (fena = T°" Ele p ea Toral REO 


Adding on both sides of (2.7) and (2.8) and using the condition 5 4 Sb, for all 


n>1, we have 
(ue b) Z bin [ellen Tai) elle 7", )]<=. (2.9) 


Since D a0 =œ by the assuinption (ii), we have 


tmint [e (fe, T7". al. -n= 210 
By virtue of the continity and monotonicity of function g, we assert that there 
exists a subsequence (X, } of {x,} such that 

lo, LS E O. -Ty | 0 as joo. (2.11) 
By the assumption (iii), we see that 

ly, =x, |s 5 ls. -T"x, [+ cl [>0, as nb. (2.12) 


It follows from the uniform equi-continuity of T that 


[Ty, - Tx, | >0 as n=. (2.13) 
Now we observe that 
lx, - 7^2 [|< a, - T7», | ry, -";]. (2.14) 
It follows from (2.11) and (2.13) that 
le, = Ip | — 0 and | LTA | —0 as joe. 
Since [r^ Xn, = X. | —0 as je, we have 
(2.15) 


X . 
| DE Ika |r" = de |y TASA as Jo. 
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E. - un 


a uM 


It follows from the uniform equi-continuity of T that 


pre 
[ieee a T X, | x 0 as J TREO (2.16) 
Again, from above inequalities, we observe that ' 
è Sin 
pra, ma. | e] ms rn nl | 
It follows from the uniform equi-continuity of T' that 
[Ce -Ts, |^ 0 as j>, (2.18 
It f 
and 
ny nj i fi 
fn Paks E ns] as| 
Therefore, it follows from (2.17), (2.18) and the above inequality that |: 
x, -Tx, |—0 as joo, 
Is... Ta, | J 9 (220 | The 
This completes the proof. 
Let {z,} be a given sequence in K. Recall that a mapping T:K—K with 
F (T) 0 is said to satisfy condition (A) [9] if there exists a nondecreasing function 
with f :[0,00) > [0,co) with f(0)=0, Ar)20 for all r e (0,0) such that 
lz, -T2, | 2 fla (z,,F(T))) for all n 2 1, 
where d(z,,F(T))=inf lle, - pl: pe F(T}. 
| = By using Theorem 2.1, we have the following : 
Wee cM E be a uniformly convex Banach space, K be a nonempty conve 
i (0) . e D . i 
| and T:K—K be a unifromly equi-continuous and asymptotically qua I 
| none i i i wap = ° 
| xpansive mapping with F(T) +4 and a sequence fk, ) =[1,00) with lim! p lo 
i d * > 2 3 à | J*1 
and $7 (Me 1)« e. Let {æn} be a sequence in K defined by (2.1) with the followin 
| restrictions: 
| W 0<b,<b<1 and b, «b, forall n>1 
(ii) " b => (iii) limb? 20, — EROS | 
£u n , (iii) ae) 3 (iv) , uat os and xo AK co. | Hen, 
satisfi iti x JESS 
es condition (A), then the sequence {x converges strongly ® | "he 


es ive Ren 
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19) 


20) 


il 


proof. It follows from Theorem 2.1 that 
liminf |, e | =0. 

Since T satisfies condition (A), we have 
lim inf f (d(x,,FT))- 0. 

From the property of f it follows that 
lim inf d(x,,F(T))- 0. 


It follows from Lemma 2.1 that d(x,,F(T)). 0 as noo. Now, we can take an 
infinite subsequence (X, }of{x,}and a sequence {p,}c F(T) such that 


ls, -p;| S2”. Set M= expl >, (k -1)} and write n,+1=n,,, for some />1, 


Then we have 


2 
s ha ai= [pez = Pj 


< [1 + fem. - 1] 


< exp|(F sn ls. T p| 
I 


con if ep 


m=0 


iod -= Dj =e, a 


Xn sit 7 Dj | 


<M : (2.21) 
2 ` 
It follows from (2.21) that 


+ Xn jt - pl 


DL 
2M +1 (2.22) 

Saja 
Hence {p;} is a Cauchy sequence. Assume that P; 
p F(T) is closed, which implies that x; > P as j > œ~. This completes the proof. 


Smark 2.1. We note that if [af 1S completely tertmtlidusithen it must be 


P ~ pj s 


— p as j 2 «. Then peF(T) 


ee 


demicompact [8], and if T is continuous and demicompact, it must satisfy condition ifm 
(A) [4,9]. In view of this observation, Theorem 2.2 improves the corresponding 
result of Liu [3] in the following aspects: | 
(i) K may be not necessarily compact or bounded, r D) 
(ii) T may be not uniformly Holder continous. 
Remark 2.2. Our results improve and generalize the corresponding results or 
[6,7,8,10,11,12] and many others from the existing literature. 
Remark 2.3. Our results also extend the corresponding results of Cho et al. [1] t; 
t h e G est Qecontinuous asymptotically quasi-nonexpansive 
mappings. 
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ABSTRACT 

In this paper, we develop a model to study the spread of HIV infection, which 
can cause Acquired Immunodificiency Syndrome (AIDS), through Vertical and 
horizontal transmissions and introduce the concept of sexually immature and 
mature individuals by considering maturation rates in heterosexual community 
and partially analyzed. We obtain equilibrium points of the system at two states 
(Disease-free and Endemic). We investigate the criteria for existence of endemic 
Steady state of the system. We determine local and global dynamics of these steady 
States of the system and conclude. 

2010 Mathematics Subject Classification : Primary 92B05; Secondary sana 
Keywords: Vertical and Horizontal Transmissions, Heterosexual Community, 
Asymptotic Stability, Epidemiological Parameters. 

1. n NU Immunodeficiency Syndrome (AIDS), 
has Shown a very high degree of prevalence in populations all over the world, 
Which is caused by Human Immunodeficiency Virus (HIV). The nature of human 
interactions, the uncertainties in the current estimates of epidemiological 
Parameters and the lack of enough reliable data make it extremely difficult 3 
understand the dynamics of the virus transmission without the frame "BE 
Provided by mathematical modelling. The study of mathematical modelling is 


= Wee v 
helpful in determining the. demographic end ecenamiojmpachof the ep idemic whic 


ee E A 


is turn help us to develop reasonable scientifically and socially sound investigatio, | 
plans in order to reduce the spread of the infection. | 


| — 


In recent decades, several mathematical modelling studies have bom | ( 
conducted to describe the transmission dynamics of HIV infection for homogeneous r d 
and heterogeneous populations, Anderson et al. ([1][2]), Bailey [5], Knox UM 7 


Pickering et al. [21], May and Anderson ([17][20]), Grant et al. [12], Hethcote (13) 
Anderson et al. [3], May [18] Castillo-Chavez, Cooke, Huang and Levin (171,181.19), 
Sun [23], Chen [10] and Hethcote [13]. Si 

In particular, Anderson et al. [1] described some preliminary attempts to 
use mathematical models for transmission of HIV in a homosexual community 
May and Anderson [17] presented simple HIV transmission models to help clarify 
the effects of various factors on the overall pattern of AIDS epidemic. Blythe and E 


Anderson (1988) considered HIV transmission models with four forms for the an 
distribution of incubation period by assuming that the infectious period is equal 2 
to the incubation period. Castillo-Chavez et al. [7] analyzed a model where the ^ 
mean rate of acquisition of new partners depends on the size of the sexually active j | 
population. Most of the above mentioned models consider only one population but =" 
HIV transmission takes place in the population that are heterogeneous in a variey | ; 
of ways and this aspect should be taken in modelling HIV, Knox [16], Colgete et al = 
[11], Jacquez et al. [14], Koopman et al. [15]. i 

2. Model Formulation. Let us consider a heterosexual community of size x 
P with uniform promiscuous behaviour and taking only heterosexual encounters bir 
and assume that the birth and death rates are same, making community size to be ie 
a constant. We have assumed that infection passes in the considering population an 
through the member of one male or female class to the other female or male clas | Th 
respectivley. The infection can also be transmitted vertically to the offspring o | De 
infected mother, | 

Let any instant of time ¢ this considering community be subdivided into $ 


classe of Sı (t) mature male susceptibles, I 1& mature male infective having Hn Mc 
infegtion, S(t) mature female susceptibles, I(t) mature female infective havi" 
HIV infection, X(t) immature susceptibles and Y(t) immature infectives havi"! set 
HIV infection. The susceptibles become infected with transmssion efficiency 42” | 


immature susceptibles and infectives being sexually matured at rates m andil i 
respectively. kS 

This leads to the following system of ordinary differential equations z : 
dS : DI : | 
cu c X. bs, q il, omY —b'1,, ks 
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DS, . i1 +(1-0)mX -58, DA LAS, +( 
t 


1-a)m'y — OM 


E 08 + pb'I-omX -(1-0)mX, & = qj'1 con Y -(1=0)m0 Y, (2.1) 
with initial data 
S (0) Sio > 0,5, (0) = Sy, > 0,1, (0) = Z > 0,1, (0) - 1,, > 0, X (0) - X, >0, 

Y (0)=Y,>0,P=S+1+X+Y,S=S,+8,,]= 1,+1,,1= p+q,0<0.<1 (2.2) 


where 

a and (1-a)= The proportions of male 
and female hosts, who converts from 
immature class to mature class 
respectively. 

band b' are Birth rtes of immature 
susceptible and infectives respectively. 

P= The fraction of newborn offspring of 
infective parents, who are susceptible at 
birth. 

q= The fraction of newborn offsprings 
of infective parents, who are Bue at 
birth. 

b and b' are also Death rate of mature 
susceptibles and infectives respectively. 
The above mathematical modelling can 
be understood by Fig.1. 


3. Equilibrium Points of 
Model. The equilibrium points of the 
Model can be derived from the following 
Set of equations 


kS I, + omX° -bS =0, 
AST e om'Y' —b'I' =0, 
TES, «(1-o)mX' —bS; =0, 


EST +(1-0)m'Y' -bL =0, 
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/ 


| 
u- H de 


& 


qb 


i j 
LSTA 


pb 


b(S; +S;)+ pb'(1 + I,)—amX" —(1-0)mX" =0, 


qb'(1, + I,)-om'Y -(l-a)m'Y =0, whi 
SOL tl, bX +Y =p. (31) 
Disease Free Equilibrium Point E, (S',1,,S,,1,,X ,Y ). We 
oy =o dX 
EL o NU WP. P Qr 
S, = „S, = = =) H 
m+b m+b mtb en 
y e 2 Also 
Endemic Eguilibrium Point E, (SAA 
yA 
S. ta? +4a,a, 
i 2a, 4 
o o : km 5 whic 
S; »0; if ka(mX +m'Y )>bb'+ Ly : 
1 cal 
Ss aut a; + 4a,a, By I 
ts, -e cond 
Hen 
S; »0; if h(1-0)(mX° +m'Y")> bo 7 y" Furt 
q 
«_ Oa, -bS ae 
I, =—— 
2a, By D 
unde 
Ta A J : Henc 
120, if a(mX tm'Y 6S, We c 
y -0)a, - 5S; 
b' 2 
1,20, if (1-0)(mX” +m'Y") busi. | dy" 


Now we investigate a criteria for endemic steady state E, to exist. Here V 


need that the system of equatio 
ns (3.1) s iti i 
From the Ser OF edi Gud Rr tais hould: haveacpesitine solution. 


X -b-bY, 


which implies 


(3.2) 


X b, when Y” 30 and Y + 6,/b,when x’ 30 
.1) 


d o oa 

We can also obtain A =-b,, e 

a 0, if (m'+qb)b> pb'm' . 

Hence, we find that X" is decreasing function of Y 

Also, form the set of equations (3.1) we get 

aX * [2ka,, —(a,, *a4)]Y^ + 2a,, = | [enc * (0a, ha) Y ~a,,} +B'ay 5X" | 

«t -aJay X + (a -a)a, - kap) Y“ 5 as) ES b'awa X |, (3.3) 
which implies 
EVI a13 312 ^ 2 

o) (X ) sting (X ) SE LY (X ) +0)X +a,,=0, when Y =(, (3.4) 

By Discarte's rule of sign the equation (3.4) has at least one positive root under 

condition 2bb'> 1, Let us denote that positive root by @,. 

Hence, x* _, Q, when Y^ 50. | 

Further, when X 50, we get | 
\ 3 4 LJ i e. e | 
| n(Y ) + yg (Y ) + ds (Y y + Gay +a, =0. (3.5) 

By Discarte's rules of sign, the equation (3.5) also has at least one positive root 

Under condition 2bb'> 1. Let us denote that positive root by Qo. 

M Y s Q, when X' 50. 

* can also obtain 
a s, (1 a) ay, @33 035 k + 2ka, -( +a ) 
dx qu Us 1 1| Sap * uz (ras, 2 Ona t Qs 
PS ME DUE A AAN 
oa E l| an, Ge |y: 
ye Qio | Ds of ol qur gia J^ Ca 
l 1/2 1/2 Dan? a3 
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LU 
a. 
ke Ody, (1-0) ay, m dos y Dus 26 
<0 if Np 12 i: Gs a; 
dy” Y a. Q3 2 " 


E Het Jia > (an tan) 


1/2 
32 E 
and 
c lO) Gen Lil a Are |p: 
Ge Ady ( 3 4 =} LAS |b'|«1 
1/2 1/2 21 a.. a 
A», As, 32 34 
or if 


Gag, Gan 
E tg Jta, > (a, +a) 


inn, gy 


Oley LO) l(ay ag yl. 
ai? a 2 al? qi? í a 


Hence, we find that X^ is 
decreasing function of Y*. Form 
above, we see that the two isoclines 
given by (3.2) and (3.3) will intersect 


BES sD t 


b, 
provided that Q, >b, and DE Q.. 


Hence, with the above (tb) 
considerations, equation (3.2) 
represents X^ as decreasing from bi 
and equation (3.3) represents X' as 
decreasing from Qi. 

Therefore, two isoclines must 
intersect provided that intersection 


CE a SANA A EA 


stea 
anal 
suffi 


a NN 


exists in the positive Y-X' plane, shown in Fig. 2. 


. Note : The values of constants a; (1=1,2,....35) and b; (t=1,2) are given in the 
| Appendix. 

4. @ualitative Analysis. Now to determine local and global dynamics of 

steady states E, and E, of system (2.1) we use, Lyapunov's Second Method. We 


analyze local dynamic of above steady states (Ey and E,) and therefore find 
sufficient conditions under which E, and E, are locally asymptotically stable in 
the form of the following Theorems 4.1 and 4.2 respectively: 

Theorem 4.1. Let the following inegualities hold 

b'» (kS, + om')/2 (4.1a) 
b>{kS, * (1- o))m] (4.1b) 
b» (AS, «(1-o)m')/2 (4.10) 
m»(b-pb) ` (4.1d) 
m'»qb'. Jj (4.1e) 


Then E, is locally asymptotically stable. 
Proof. Using the following positive definite function in the linearized form of 
System (2.1), 


V = [n + An? + An? + An? + An; + A;n; |/2 , (4.2) 


Where, A; are arbitrary positive constants (i=1,2,...,6), a it can be checked ipa the 
derivative of V, with respect to ¢ under the conditions (4.1) is negative definite. 
| Hence in view of theory of stability, E, is locally asympotically stable. 


Theorem 4.2. Steady state E, of the system (2.1) is locally asymptotically stable 
if 


0 (kS? + kI; tam/)/2, (4.3a) 

(E x0)3 Pass +(1=0)m]/2, (4.3b) 

9» (kS; + hI: +(1-a)m')/2 , (4.30) 

12 (d+ pb), (4.3d) | 
E| usano (4.3e) 
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L40 


Proof. Similar to Theorem 4.1. d 
Now to show that steady states Ej and E, of system (2.1) are globally 


asymptotically stable, we coste a region of attraction for the system (2.1) ; in in tl 
the form of following - du, 
Lemma 1. Consider the set di 
R - ((S,1,S, 1, X, Y):0« Sim SS, < P,0< 1, <1, SP,0« So, SS, <P, m 

dt 


O«L, SI, <P,0<x,0< Y 


2m 
Which is a region of attraction for all solutions initially in the positive orthant | dus 
= Som = I 
Now we obtain here the conditions for asymptotic stabiity of positive steady | du, 

states Ey and E, of system (2.1) in non-linear (global) case in the form of following | dt 
Theorem 4.3 . The steady state E, of system (2.1) is non-linearly (globally) 


where S = [,, are positive constants. 


ln lm 


asmptotically stable in a region R, given by Lemma 1 where the following dus 
: d 
conditions are satisfied : i 
0<B,<2, : (45) | du, 
b'>(kP+am')/2 ` (48b) ge 
b>{kS,,, +(1=0. os 2 (4.5) | Lem: 
dV, 
b» ARP +(1-a)m'}/2 4501| —2 
{ i 3/2, dt 
m>(b+ pb'), b; (4.50) 
m'>qb', 5 (4.50 
| where B, is arbitrary positive constant. 
Proof. Consider the following positive definite function about Eo 
V, = [uz + Bu; + Bou; + Bau; + B,u? + By; |/2. e (4.6) Mani 
| Here 3 - a, 
B; are arbitrary positive constants (i—1,2,...,6). dt 
Take perturbations in E,(S;,1\,8,,1;,X 24 ) as u, (t) u (t), us (t),ua (t),us(¢) j hec 
TN 


us (t) respectively and putting 
S, =S +u (t) =I tu, (t), 


S, = S, ds Us (t) I, ce Guay (Elkangri Collection, Haridwar. An eGangotri Initiative 


> ÓN 


eX eu (Y =Y eu), 


lly in the system (2.1), we get 
| in A 

dy by, —h(S; +u, Ju, tam 

sa 1 ( 1 3 4 Us, (4.7a) 

d ? D i 1 

us u, * k(S, +u, Ju, + om Us, (4.7b) 
nt, dis eS, +u Ju — bu; +(1-a)mu 

di 2 Us Ju, 3 BG (4.7c) 
ady du id 1 ! 
3 E k(S; + Uy Ju, -b'u, *(1-a)m'u,, (4.7d) 
ly) 
| du; 1 
Ing A = bu, + pb'u, + bu, + pb'u, -mu,, (4.7e) 
5a) Ug 1 

meo U, + qb'u, -m'us. (4.7£) 
5b) | 

D i f f e Varwith respect to £ along the solution of (2.1) and using 


$)| Lemma 1, we get 


dV, 
5d) uo —bu? -b' Bu? - bB,u? —b' Bui] —mB,u? - m' Byug —kS, uu, + (omn. bB, usu; 


5e) 
i +kPB,u,u, + om! Buuyu, - kSyp, Byusu; + | (1-01) mB, + BB, Jusu; + kPB,U,t, 
5 
+(1- a) m'B,u,U, + pb' B,uzu; + pb' Byusu, +qb'B,u,u,+qb'B,uu,;. (4.8) 
Applying the inequality tab x (a? +0*)/2 and making algebraic | 
p manipulation, we derive | 
Tes Bu? a a 2 2 + Bus | (4.9) | 
d ` [ u^ + Bu; + B,guz + Bu + Byus + Pigs |, | 
nd | Where f 


Baa 
u=b—(kS,,, +0m-+bB,)/2, 


Bi 3 b'B, -(kPB, * kS,, B, +oun' B, + kPB, + pb' B, #gb'B,)/2, 


Bu 
13 bB, = LN B, + Ga, hills. Haridwar. An eGangotri Initiative 


— 


B, =6'B, -[RS,,, + RPB, +kPB, +(1-a)m'B, +qb' B, + pb'B, /2, 
B, <mB,—{am-+(1-0)mB, +2(6+ pb)B 2 


B,, -m'B, -|am'B, +(1-0)m'B, + 2qb' B, |/2. 


i, ae : is o 
From (4.9), it can be shown that, T is negative definite under following condition, 


b» (kS; +0m+bB,)/2, (4.109) 
b' B, » (kPB, +kS,,,B, + kPB, + om' B, + qb'B, + pb' B], (4.10b 
bB, » [kS,, B, * (1-o)mB, + BB, ]/2, (4.100 
b B, >[RS,,, + kPB, +kPB, « (1- o) m' B, * qb' B; + pb' B, |/2, (4.10d | 
mB, >[om+2bB, + 2pb'B, +(1-0)mB, |/2, (4.10 
m'B, » {am'B, +(1-a)m' B, +2qb'B,)/2. (4.100 
V; 


The above sufficient conditions for to be negative dfinite, may be further 


dt 
manipulated to get the following simplified conditions: 

In (4.102) choosing B, as 
0« B, «2, „(4.110 
and in (4.10b) choosing B, as 

Bs RS,,,B, + kPB, + pb' B, #gb' B; 

b'- (kP + om')/2 j 

the condition (4.10b) reduces to 


y 
b'>(kP+om')/2 . 2 Gy 


In (4.10c) choosing By as 


bB,)/2 
B ( 4 

E b-{RS,,, +(1-a)m}/2 ? 
the condition (4.10c) reduces to 


b> (ks + (1 = okeh ral Kangri Collection, Haridwar. An eGangotri Initiative 
2m d . 


E IES 


In (4.10d) choosing Bj as 


(kS,,, + kPB, + pb' B, + qb' B;) 
(B7 (RP + (La m2 > 


the condition (4.10d) reduces to 


b'>(kP+(1-a)m')/2 . (4.110) 
E^ In (4.10e) choosing B, as 
0a) (om+(1-0)mB,)/2 
1g m-(b+ pb") É 
0b 
the condition (4.10e) reduces to 
Oc 


m>(b+ pb"). " (4.11e) 
0d In (4.10f) choosing B; as 


(| y „ (om'B, #(1-0)B,)/2 

> 2m'- 2qb' : 

10) the condition (4.10f) reduces to 
m'»qb'. (4.11f) 

s Hence, disease free steady state E, of system (2.1) is non-linearly (Globally) 
asymptotically stable in the region R, under the conditions given by (4.11), provig 
Theorem 4.3. 

Ja Similarly we also determine that steady state E, of system (2.1) is also globally 
asymptotically stable in a region R, given by Lemma 1 under the following 
Conditions : 
kb» (om +kPD, +hS; +0D,)/2, (4.122) 

tb) (kr + b) > (RS,,,, +(1 = o)m)/2 : (4.12c) | 
b» + kp +(1-a)m'}/2, 4.120) | 

A m> (b + pb") (4.12e) f 
(4.12f) 


m'»gbr. E 


Where D, and D itrary positive constants. 
: 5. i and D, are eher po el, disease tree aná eddëémie'steady states, have 


A PPM 


been obtained, which are shown to be both linearly (locally) and non-linear], 
(globally) asymptotically stable under the conditions involving disease related dy 
parameters. 

From the qualitative analysis of the disease free steady state Ep, it may he » 4 
concluded that the infection will die out eventually in the underlying Population 
and only mature and immature population of susceptible males and females wil] 
exist. From the qualitative analysis of the endemic steady state Æ}, it may be 
concluded that the infection will remain always in the population provided the b, = 
following disease related parameters satisfy 


almx +m'Y”)>bS;,(1-0)(mX" +m'Y")>bS;,2bb' > 1,m > (b+ pb") and m» qb, 


APPENDIX Sci 
a, = kb, a, = oa X * (o4, -ka,)Y -=a d, = alax ta,Y), and 
a= (1- aja, X + [a 3 0.)a,, = hay, ]Y Sa =(1- 0) (a, X "i aY’) = 
' 1 Š Š E 5 
a, =mb',a,=m',a,=mX +m'Y a =km,a,, = 4ba,a,, = —,a,, = 4kpb',a,,km', i 
q [4] 
a; = bb',a,, = 4b(1—-0),a,, =a,, + 0,5 — 2ka,,,a,, =(1=0.)a,, — ka, yy = 00,4 — hy, b 
App = 0,0, = 2(1- 1) a,,0,, + 2k,,0,5, ny = Quya = 202,,0,, + 2ka,,0,;, | [7] 
an “(fo «01-29 (1-9)! Ja -2(0*+(1-0%)]aj ah), 
— 3 2 ‘ 
| da =[0 +(1-a) (0(1=0)+5/2)]4a%,0,, +[0*+(1-0) ba, +80, a5, | 
i - 3 


an = [o «(1-0)']aj, - 4|2{o? +(1-a)}- 5 laha 


155 


? 


= 2 2 5 
Qs; = 2fa t (1 = 0) lahat, = 6da,,as,, Qs = da? (1 = 4a.) 


Bee SAS pa 2 2 2 
05; = Ay, FG + yy + 201,07 — 20; (ai, + oh) 


—-942 2 e j 
Qs = 24,500 t 2a, —- 8g? e 2 2 6 : 9 2 | 
9721 1% 2 Aiga + Aig + ay, (a5, + a) = 40,5055 (ais + aj y 


n 
y 


A 4 Mp 
Ang = Aja FA. +2 2 Qa? a? yp 2 
29 18 19 Aas + 20,,.a.. + 2 y? 2 4 1 t 
1815 19 “15 24aia;. -2 apan A a 2 4 al, (ais + 9 
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Zee 272 3 2 
= 5 S — 2 2 x z 
a, = 2,40 + 205,0, - 320,4; 2[ anat, + Op dis + Lay (a, a.) - 8a,,4;, | 
a i m , 


2 2 AP ! > z 
|» Ap 4a; (1 ES Lara ) a; = Any t bana, X, Q5, = Oy X +(00,, Ra, )Y -a 


15? 


ly = a5, + O'o X as, = (1-e)a, X" + (a -0)a,, hays ve -a 
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Pb b (1+qb'/m'- pb'/b)m' 
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i ABSTRACT 

An exact value of the unit of length measurement used in Indus-Saraswati 
Civilization, has been determined from the precise scale discovered by Ernest 
Mackay in the 1930-31 season excavation at Mohenjodaro and further correlated 
with the present day units of measurement. It has been calculated and shown that 
the speed of light as given in the Vedic literture, when referenced with this erstwhile 
unit of length measurement (used in Indub-Saraswati Civilization), works out to 
be precisely equal to the speed of light as per modern measurements. The present 
paper is a step-wise process followed to unveil this equality. 

2010 Mathematics Subject Classification : Primary 01AXX Secondary 01432 
Keywords: Unit of length measurement, Indus-saraswati civilization, speed of 
light in vedic literature. 

1. The Precise Scale. In his 1930-31 season at Mohanjo-daro, Ernest 
Mackay discovered a broken piece of shell bearing 8 divisions of precisely 6.7056mm 
each, with a dot and circle five graduations apart, which suggests a decimal system. 
However, attempts by Mackay, to relate such a unit to dimensions in Mohanjo- 
daro, were not very successful (Michel Danino [2]) and thus were abandoned. 

2. Units of Length in Chanakya's Arthashastra. Chanakya was the 
Political of the legendary monarch Chandragupta Maurya of zn century BC. He 
Was aman learned in may disciplines and wrote the famous treatise on economics 
called the Arthashastra-meaning the books of money). In Arthashastra, Chanakya 
mentions two types of Dhanushas as units for measuring length and eee 
One is the ordinary Dhanusha, consisting of 96 Angulas, and the other D ge Me 
is mentioned as Garhpatya Dhanusha and consists of 108 Angulas. Chanakya also 


Mentions many other units including a Dhanurgraha, which consists a 


and a Yojana (In Sanskrit-English dictionary, by 
been defined as a mesure of distance=4 Kosas or 

Oniere William dictionary, several other books/sources 
Cgal to about 9 milSEP: Gu warsisdingl adi BOOM MR ARA otri Initiative 


Moniere William, a Yojana has 
about 9 miles. Apart from the 
also give the Yojana as 


AA 


gg 
3. Decoding the Mohan-jo-daro Scale. If we keep 10 divisions of yy, 


: spes oli 

Mohanjo-daro scale as equal to a Dhanurgraha or 4 Angulas, the precise length of T , 
an Angula works out to be 16.764mm. | 307 
A Dhanusha of 96 Angulas=96 x 16.764mm=1.609344m (1) 2 
and re 
A Dhanusha of 108 Angulas= 108x16.764mm =1.810512meters. (2) een 
1 Yojana=8000 Dhanushas (of 108 Angulas each) (3) MG 

Thus Sar: 
1 Yojana=8000 x1.810512m=14.484096km. (4) Ma: 
Further Jus 
14.484096km=9 miles, (exactly!). (b) s 

Also i Ang 
1000 Dhanushas of 96 Angulas each=1mile (6) Ten 
Interestingly, when we look into the history of mile, we find that the word at 1 

mile is derived from mile, which means a thousand. This points to universal mes 
adaptation of ancient units of length. part 


4. Corroboration from Other Sources à 

The Indus Inch: The Indus civilisation unit of length, widely known as Indus 
Inch was 1.32 Inches which is exactly equal to 2 Angulas of 16.76mm each. the 
Mohenjo-daro's Great Bath: The height of the corbelled drain forming the outlet 
of Mohenjo-daro's Great Bath [6, pp. 133-142] is about 1.8m, which is equal to a 
Dhanusha of 108 Angulas of 16.764mm each. 

Standard Street-Widths. Kalibangam, a city in the Indus-Saraswati Civilization divi: 
(in Rajasthan) had street widths [8] of 1.8m, 3.6m, 5.4m and 7.2m i.e. built to the 


scale 
standard dimensions being equal to 1 Dhanusha, 2 Dhanushas, 3 Dhanushas and 4 
Dhanushas respectively. Such widths are found at other sites also. Bigger streeló ) prec 
of Banawali [8] another town in Indus-Saraswati Civilization (in Haryana) measure Civil 
5.4m i.e. they were built with the unit of 3 Dhanushas. May 


Taj Mahal. A Persian manuscrip 


description of three principal bu 
Jamah Masjid. In the 


t "Shah Jahan Nama" contains a very particular 


L 


as well as in Gaz, are given by Hod 


0 
(80.52cm) from Table A and 31.66 inches=80.42cm from Table B. The avera£ i 


: : E > tha 
the two values in “198 Ingh esm kem iibetorperimepti liit en tion here 
Barraud [1, pp. 108-109, 258.259] in T j i 


Yi 


Columns of Agra, has taken a Gazas egual to 80.5cm which is very nearly equal to 
80.47cm as worked out above. Taking a Dhanusha of 96 Angulas to be equal to 2 
Gaz, the length of an Angula works out to 16.764mm. It shows that a Gaz of 48 
Indus-Saraswati Angulas was being used, even in the days of Shah Jahan's rule 
ie, 17th Century A.D. 

The Gudea's Rule. The Gidea's rule (2175 B.C.) preserved in the Louvre shows 
intervals in Sumerian Shusi of 0.66 inches, which is exactly equal to the Indus- 
Saraswati Angula of 16.764mm. 

Mayan Units of Measurement. Drewitt [3] and Drucker [4] made a study of the 
ancient city of Teotihuacan, belonging to Mayan Civilization, in Mexico, and 
hypothesised a unit of 80.5cm, which is very nearly equal to the Indian Guz of 48 
Angulas=80.47cm. 

Temple Wall-Engravings. Nearer home, two engravings on a wall of the temple 
at Tiruputtkali (12th Century A.D.) near Kanchipuram, show two scales [7] one 
measuring 7.24 metres in length, with marketing dividing the scale into 4 equal 
parts, and the second one measuring 5.69 metres in length and marking dividing 
the scale into 4 equal parts. It may be observed that each division of the first scale 
is precisely equal to a Dhanusha of 108 Angulas of 16.764mm each. Interestingly, 
the second scale is precisely equal to 1 times Dhanusha i.e. equal to the 
circumference of a circle with one Dhanusha as its Diameter. 

It is interesting to note here that Mackay reports [9] at Mohenjo-daro, a 
lane and a doorway having both a width of 1.42m, which is precisely equal to one 
division of the second scale at the Tiruputtkali Temple, indicating that both the 
scales were prevalent in Indus-Saraswati Civilization as well as in South India. 

It proves beyond doubt that the Units of measurement as derived from the 
Precise scale found at Mohenjo-daro were prevalent not only in the Indus-Saraswati 
Civilization, but also in South India and in the ancient Sumerian, Egyption and 
Mayan Civilizations. 

5. Speed of Light in Vedic Lit 
Mandal 1, Sukta 50, Mantra 4, which is in praise of the Sun god, Say 
(14th Century AD) writes: 

we a ER.. 

IVA wed gg Tag a aoi! 

Yor Ma maa AAR” All 


Meaning... "It is remembered that... cart 
Salutations to Thee (the Sun) who approacheth (at a speed of) 2202 yojana 


ina nimish $F "c | ae é 
e etes ant (or sunrays) that 18 meneioHed in the Shloka. 


erature. In the commentary on Rig-Veda. 
anacharya 


wy 


— re NAN — A 


This shloka is attributed [10, pp. 67] to the son of Kanva Maharshi (4000 B.C), | was 
Bhatta Bhaskara (10th Century) mentions [11] this shlok in his commentary pp | MEA 
the Taitreya Brahamana. "m 
To put it in mathematical terms, as stated by the above shlok, the speeq of P kno 

light would be : 1 
Speed of Light= 2202 Yojanas/Nimishardha (1) Spes 


We have already calculated the value of Yojana in modern unit of km in the | 
previous step (Equation 4), but we still do not know, what the Nimishardha was 
translates into. Let us dig into another Vedic text tLe Vishnu Puran to find how | Civi 


iz) 


the erstwhile units of time can be related to the modern units of time. | inch 
In the Vishnu Puran (Book 1, Chapter 3, Shloka 8,9), it is stated that: 
15 Nimishas = 1 Kashtha [1] 
30 Kashthas = ] Kala 
30 Kalas = 1 Mahurta [2] 
30 Mahurtas = 1 day and night (FERTAH) | ag 
Thus | 
one day and night = 405,000 Nimishas=810,000 Nimishardhas. (8) [5] 


(Literal meaning of Nimishardha being half of Nimisha). 
In Surya Sidhant (Chapter 1, Shloka 12), it is mentioned that 60 Nadis 6] 
constitute one Sidereal Day and Night («rera SRT). It is also well known that 
1 Mahurta=2 Ghatis or 2 Nadis. It is clear from this that in Astronomical 
calculations, the sidereal day was taken as the unit of time. A sidereal day is the 


time taken by the stellar constellations to complete one revolution around the 8] 
Earth. A sidereal day is equal to 23 hours, 56 minutes and 4.1 seconds or equivalent 
sec (Wikipedia (121). " 
Thus ; | ) 
Nimishardha=86164.1/810000 sec=0.1063754sec (9) a 
The speed of light as given in the Vedic Literature therefor comes out to be; i 
2202 * 14.484096 km 
= 75 
0.1063754 sec” 98 10 km/sec . (10) 


which is precisely equal to the speed of light as per modern measurement? 

6. Conclusions. The calculations and references int he article endeavor ” 
establish the vast reach of the Indus Valley scholars, not only academically, pr 
also geographically. Thousands of years before the modern scientists rediscove™ 
the speed of light; our ancestors knew of the exact same value and referred to it” 
an exalted property of the Sun-God to salute Him. It is also heartening to not 


thatthe modern-age concept of elahalizationand:shartmaar ives between culture 


aS 


| 


was also an established way of life in their era, as 

measurement found across far-flung cultures. 

civilization lost their influence is still a mystery, 
knowledge that : 

1. The Speed of Light as given in the Vedic literature, is precisely equal to ihe 

Speed of Light as per modern measurements. 

9. The basic unit of length measurement in the Indus-Saraswati Civilization 

was an Angul of 16.764mm. This unit was used not only in the Indus-Saraswati 

Civilization, but also in South India, and other ancient world Civilizations 

including Sumerian, Egyptian and Mayan Civilizations. 
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proven by the shared units of 
When and how these ancient 
but for now we can celebrate the 
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